Диссертация (1149979), страница 20
Текст из файла (страница 20)
.2);> er2 := int(-(a2*x+b2)*x*BesselJ(0, x*t), x = .2 .. .4);int(-(0.01979000000 x + 0.8645090000) x BesselJ(0, x t), x = 0.2 .. 0.4)> er3 := int(-(a3*x+b3)*x*BesselJ(0, x*t), x = .4 .. .6);int(-(0.03638500000 x + 0.8578710000) x BesselJ(0, x t), x = 0.4 .. 0.6)> er4 := int(-(a4*x+b4)*x*BesselJ(0, x*t), x = .6 .. .7);int(-(0.05321000000 x + 0.8477760000) x BesselJ(0, x t), x = 0.6 .. 0.7)> er5 := int(-(a5*x+b5)*x*BesselJ(0, x*t), x = .7 ..
.8);int(-(0.06953000000 x + 0.8363520000) x BesselJ(0, x t), x = 0.7 .. 0.8)> er6 := int(-(a6*x+b6)*x*BesselJ(0, x*t), x = .8 .. .9);int(-(0.09534000000 x + 0.8157040000) x BesselJ(0, x t), x = 0.8 .. 0.9)> er7 := int(-(a7*x+b7)*x*BesselJ(0, x*t), x = .9 .. 1);int(-(0.1672900000 x + 0.7509490000) x BesselJ(0, x t), x = 0.9 ..
1)> err := er1+er2+er3+er4+er5+er6+er7;> int(BesselJ(0, 1.2*t)*t^2*err/(t+.2), t = 0 .. 20);> evalf(%);> plot(int(BesselJ(0, r*t)*t^2*err/(t+.2), t = 0 .. 10), r = 1.5 .. 2);> restart;> ch1 := int(.2^2*BesselJ(0, .1*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.1/r)/(Pi*r);155(продолжение А3)> plot(ch1, r = 1.1 .. 2);> chlen1 := (.867518*.1)*ch1;> ch2 := int(.2^2*BesselJ(0, .2*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.2/r)/(Pi*r);> chlen2 := (.868467*.2)*ch2;> ch3 := int(.2^2*BesselJ(0, .3*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.3/r)/(Pi*r);> chlen3 := (.870083*.3)*ch3;> ch4 := int(.2^2*BesselJ(0, .4*t)*BesselJ(0, r*t)/(t+.2), t = 0 ..
20)-(.2*2)*EllipticK(.4/r)/(Pi*r);> chlen4 := (.872425*.4)*ch4;> ch5 := int(.2^2*BesselJ(0, .5*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.5/r)/(Pi*r);> chlen5 := (.875582*.5)*ch5;> ch6 := int(.2^2*BesselJ(0, .6*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.6/r)/(Pi*r);> chlen6 := (.879702*.6)*ch6;> ch7 := int(.2^2*BesselJ(0, .7*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.7/r)/(Pi*r);> chlen7 := (.885023*.7)*ch7;> ch8 := int(.2^2*BesselJ(0, .8*t)*BesselJ(0, r*t)/(t+.2), t = 0 ..
20)-(.2*2)*EllipticK(.8/r)/(Pi*r);> chlen8 := (.891976*.8)*ch8;> ch9 := int(.2^2*BesselJ(0, .9*t)*BesselJ(0, r*t)/(t+.2), t = 0 .. 20)-(.2*2)*EllipticK(.9/r)/(Pi*r);> chlen9 := (.901510*.9)*ch9;> ch10 := int(.2^2*BesselJ(0, t)*BesselJ(0, r*t)/(t+.2), t = 0 ..
20)-(.2*2)*EllipticK(1/r)/(Pi*r);> chlen10 := (1/2*.918239)*ch10;> funct := (chlen1+chlen2+chlen3+chlen4+chlen5+chlen6+chlen7+chlen8+chlen9+chlen10)*.1;> plot(funct, r = 1.1 .. 4);> restart; ch1 := int(1^2*BesselJ(0, .1*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.1/r)/(Pi*r); chlen1 :=(.6835683280*.1)*ch1;> ch2 := int(1^2*BesselJ(0, .2*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.2/r)/(Pi*r);> chlen2 := (.6853014820*.2)*ch2;> ch3 := int(1^2*BesselJ(0, .3*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.3/r)/(Pi*r);156(продолжение А3)> chlen3 := (.6882841701*.3)*ch3;> ch4 := int(1^2*BesselJ(0, .4*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.4/r)/(Pi*r);> chlen4 := (.6926736957*.4)*ch4;> ch5 := int(1^2*BesselJ(0, .5*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.5/r)/(Pi*r);> chlen5 := (.6987347986*.5)*ch5;> ch6 := int(1^2*BesselJ(0, .6*t)*BesselJ(0, r*t)/(t+1), t = 0 ..
20)-2*EllipticK(.6/r)/(Pi*r);> chlen6 := (.7069071689*.6)*ch6;> ch7 := int(1^2*BesselJ(0, .7*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.7/r)/(Pi*r);> chlen7 := (.7179540821*.7)*ch7;> ch8 := int(1^2*BesselJ(0, .8*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.8/r)/(Pi*r);> chlen8 := (.7333406724*.8)*ch8;> ch9 := int(1^2*BesselJ(0, .9*t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(.9/r)/(Pi*r);> chlen9 := (.7565263966*.9)*ch9;> ch10 := int(1^2*BesselJ(0, t)*BesselJ(0, r*t)/(t+1), t = 0 .. 20)-2*EllipticK(1/r)/(Pi*r);> chlen10 := (1/2*.8056100709)*ch10;> funct1 := (chlen1+chlen2+chlen3+chlen4+chlen5+chlen6+chlen7+chlen8+chlen9+chlen10)*.1;> plot(funct1, r = 1.1 ..
4);> restart; ch1 := int(1.8^2*BesselJ(0, .1*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.1/r)/(Pi*r);chlen1 := (.6117393612*.1)*ch1;> ch2 := int(1.8^2*BesselJ(0, .2*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.2/r)/(Pi*r);> chlen2 := (.6138606728*.2)*ch2;> ch3 := int(1.8^2*BesselJ(0, .3*t)*BesselJ(0, r*t)/(t+1.8), t = 0 ..
20)-(1.8*2)*EllipticK(.3/r)/(Pi*r);> chlen3 := (.6174944825*.3)*ch3;> ch4 := int(1.8^2*BesselJ(0, .4*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.4/r)/(Pi*r);> chlen4 := (.6227999047*.4)*ch4;> ch5 := int(1.8^2*BesselJ(0, .5*t)*BesselJ(0, r*t)/(t+1.8), t = 0 ..
20)-(1.8*2)*EllipticK(.5/r)/(Pi*r);> chlen5 := (.6300567332*.5)*ch5;157(продолжение А3)> ch6 := int(1.8^2*BesselJ(0, .6*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.6/r)/(Pi*r);> chlen6 := (.6397507039*.6)*ch6;> ch7 := int(1.8^2*BesselJ(0, .7*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.7/r)/(Pi*r);> chlen7 := (.6527773158*.7)*ch7;> ch8 := int(1.8^2*BesselJ(0, .8*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.8/r)/(Pi*r);> chlen8 := (.6709927329*.8)*ch8;> ch9 := int(1.8^2*BesselJ(0, .9*t)*BesselJ(0, r*t)/(t+1.8), t = 0 .. 20)-(1.8*2)*EllipticK(.9/r)/(Pi*r);> chlen9 := (.6992079006*.9)*ch9;> ch10 := int(1.8^2*BesselJ(0, t)*BesselJ(0, r*t)/(t+1.8), t = 0 ..
20)-(1.8*2)*EllipticK(1/r)/(Pi*r);> chlen10 := (1/2*.7658327908)*ch10;> funct2 := (chlen1+chlen2+chlen3+chlen4+chlen5+chlen6+chlen7+chlen8+chlen9+chlen10)*.1;> plot(funct2, r = 1.1 .. 4);.