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

k=-1.0;

sum=0.0;

s=0.0;

do{

k=k+1.0;

sum=sum+s;

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

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

return(sum);

}

}

}

complex J_der(double x,double n)

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

complex J_der_der(double x,double n)

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

complex Ne(double x,double n)

{

complex sum1,sum2,sum3,s,ss,sss;

double k,nn,i;

if(n<0.0) return(pow(-1.0,-n)*Ne(x,-n));

else

{

if(n>1.0) return(2.0*(n-1.0)/x*Ne(x,n-1.0)-Ne(x,n-2.0));

if(n==0.0)

{

n=0.0;

sum1=2.0/M_PI*(C+log(x/2.0))*J(x,n);

sum2=0.0;