491
z*R[0]/ k[1]^(1/2))-BesselJ(1, z*R[0]/k[1]^(1/2))*BesselY(1, z*R[1]/
k[1]^(1/2)))* k[2]^(1/2)*BesselY(0, z*R[1]/ k[2]^(1/2))-
k[1]^(1/2)*(BesselJ(0, z*R[1]/k[1]^(1/2))* BesselY(1, z*R[0]/k[1]^(1/2))-
BesselJ(1, z*R[0]/k[1]^(1/2))*BesselY(0, z*R[1]/ k[1]^(1/2)))*
h[2]*BesselY(1, z*R[1]/k[2]^(1/2)))+(-h[1]*(BesselJ(1, z*R[1]/
k[1]^(1/2))*BesselY(1, z*R[0]/k[1]^(1/2))-BesselJ(1,
z*R[0]/k[1]^(1/2))*BesselY(1, z*R[1]/ k[1]^(1/2)))*k[2]^(1/2)*BesselJ(0,
z*R[1]/k[2]^(1/2))+ k[1]^(1/2)* (BesselJ(0, z*R[1]/k[1]^(1/2))*BesselY(1,
z*R[0]/k[1]^(1/2))-BesselJ(1, z*R[0]/ k[1]^(1/2))* BesselY(0,
z*R[1]/k[1]^(1/2)))*h[2]*BesselJ(1, z*R[1]/ k[2]^(1/2)))*BesselY(0,
z*R[2]/k[2]^(1/2)))/k[3]^(1/2)):
> R:=array(0..3):R[0]:=r0: R[1]:=R[0]+d[1]:
R[2]:=R[1]+sum(d[n2],n2=2..n1-1):R[3]:=d[n1]+R[2]:
> h[1]:=s[1]:h[2]:=(R[2]-R[1])/sum(d[n3]/s[n3],n3=2..n1-1):h[3]:=s[n1]:
> k[1]:=a[1]:k[2]:=(R[2]-R[1])*h[2]/ sum(d[n4]*s[n4]/a[n4],n4=2..n1-
1):k[3]:=a[n1]:
> z:=array(1..5):
> i:=1:x:=0.01:
> p:= array(1..15):p[1]:=0.01:
> while i<13 do if
evalf(F(x,l[1],l[2],l[3],h[1],h[2],h[3],k[1],k[2],k[3],R[0],R[1],R[2],R[3],al)
*F(x+0.1,l[1],l[2],l[3],h[1],h[2],h[3],k[1],k[2],k[3],R[0],R[1],R[2],R[3],al))
>0 then x:=x+0.1 else p[i+1]:=x: p[i+2]:=x+0.1: i:=i+2: x:=x+0.1: fi: od:
> for n to 5 do
z[n]:=fsolve(F(z,l[1],l[2],l[3],h[1],h[2],h[3],k[1],k[2],k[3],R[0],R[1],R[2],R
[3],al)=0,z=p[2*n-1]..p[2*n+2]):od:
> z[1]:=z[1]:z[2]:=z[2]:z[3]:=z[3]:z[4]:=z[4]:z[5]:=z[5]:
> y:=array(1..5,1..3):