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

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

sum2=-sum2/M_PI;

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);

}

if(n==1.0)

{

n=1.0;

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

sum2=0.0;

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

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