Курсовая работа: Моделирование рассеяния плоской упругой продольной волны на упругом однородном изотропном цилиндрическом слое

k=-1.0;

sum3=0.0;

s=0.0;

do{

k=k+1.0;

sum3=sum3+s;

s=pow(-1.0,k)/fact(k)/fact(n+k)*pow(x/2.0,2*k+n);

ss=0.0;

for(i=0.0;i<(k+1.0);i=i+1.0)

{

sss=0.0;

for(nn=1.0;nn<(n+i+1.0);nn=nn+1.0)

sss=sss+1.0/nn;

ss=ss+sss;

}

s=s*ss;

}while(module(s)>=eps);

sum3=-sum3/M_PI;

return(sum1+sum2+sum3);

}

}

}

complex Ne_der(double x,double n)

{ return((Ne(x,n-1.0)-Ne(x,n+1.0))/2.0); }

complex Ne_der_der(double x,double n)

{ return((Ne_der(x,n-1.0)-Ne_der(x,n+1.0))/2.0); }

complex H1(double x,double n)

{ return(J(x,n)+iii*Ne(x,n)); }

complex H1_der(double x,double n)