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

complex A1_n[K+K+2],A2_n[K+K+2],A3_n[K+K+2],A4_n[K+K+2];

complex B1_n[K+K+2],B2_n[K+K+2],B3_n[K+K+2],B4_n[K+K+2];

complex A[MM][MM];

complex F[MM];

complex X[MM];

float a[N][N];

float f[N];

float x[N];

void Matrix_A_F(float n)

{

A[0][0]=k1*H1_der(k1*r1,n);

A[0][1]=iii*n/r1*H1(xi1*r1,n);

A[0][2]=0.0;

A[0][3]=0.0;

A[0][4]=-k3*J_der(k3*r1,n);

A[0][5]=-k3*Ne_der(k3*r1,n);

A[0][6]=-iii*n/r1*J(xi3*r1,n);

A[0][7]=-iii*n/r1*Ne(xi3*r1,n);

F[0]=-pow(iii,n)*k1*J_der(k1*r1,n);

A[1][0]=0.0;

A[1][1]=0.0;

A[1][2]=k2*J_der(k2*r2,n);

A[1][3]=iii*n/r2*J(xi2*r2,n);

A[1][4]=-k3*J_der(k3*r2,n);

A[1][5]=-k3*Ne_der(k3*r2,n);

A[1][6]=-iii*n/r2*J(xi3*r2,n);

A[1][7]=-iii*n/r2*Ne(xi3*r2,n);

F[1]=0.0;

A[2][0]=iii*n/r1*H1(k1*r1,n);