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

A[2][2]=0.0;

A[2][3]=0.0;

A[2][4]=-iii*n/r1*J(k3*r1,n);

A[2][5]=-iii*n/r1*Ne(k3*r1,n);

A[2][6]=xi1*J_der(xi3*r1,n);

A[2][7]=xi3*Ne_der(xi3*r1,n);

F[2]=-pow(iii,n+1.0)*n/r1*J(k1*r1,n);

A[3][0]=0.0;

A[3][1]=0.0;

A[3][2]=iii*n/r2*J(k2*r2,n);

A[3][3]=-xi2*J_der(xi2*r2,n);

A[3][4]=-iii*n/r2*J(k3*r2,n);

A[3][5]=-iii*n/r2*Ne(k3*r2,n);

A[3][6]=xi3*J_der(xi3*r2,n);

A[3][7]=xi3*Ne_der(xi3*r2,n);

F[3]=0.0;

A[4][0]=2.0*M1*k1*k1*H1_der_der(k1*r1,n)-L1*k1*k1*H1(k1*r1,n);

A[4][1]=2.0*M1*iii*n/r1*(xi1*H1_der(xi1*r1,n)-H1(xi1*r1,n)/r1);

A[4][2]=0.0;

A[4][3]=0.0;

A[4][4]=-2.0*M3*k3*k3*J_der_der(k3*r1,n)-L3*k3*k3*J(k3*r1,n);

A[4][5]=-2.0*M3*k3*k3*Ne_der_der(k3*r1,n)-L3*k3*k3*Ne(k3*r1,n);

A[4][6]=-2.0*M3*iii*n/r1*(xi3*J_der(xi3*r1,n)-J(xi3*r1,n)/r1);

A[4][7]=-2.0*M3*iii*n/r1*(xi3*Ne_der(xi3*r1,n)-Ne(xi3*r1,n)/r1);

F[4]=-pow(iii,n)*(2.0*M1*k1*k1*J_der_der(k1*r1,n)-L1*k1*k1*J(k1*r1,n));

A[5][0]=2.0*M1*iii*n/r1*(k1*H1_der(k1*r1,n)-H1(k1*r1,n)/r1);

A[5][1]=M1*(-xi1*xi1*H1_der_der(xi1*r1,n)-n*n/r1/r1*H1(xi1*r1,n)+xi1/r1*H1_der(xi1*r1,n));

A[5][2]=0.0;

A[5][3]=0.0;