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

for(i=0;i<N;i++)

x[i]=0.0;

for(i=N-1;i>=0;i--)

{

c=0.0;

for(j=i+1;j<N;j++)

c=c+a[i][j]*x[j];

x[i]=(f[i]-c)/a[i][i];

}

}

void Complex_Gauss(void)

{

int i,j;

complex sum;

for(i=0;i<MM;i++)

{

a[2*i][0]=real(A[i][0]); a[2*i][1]=-imag(A[i][0]);

a[2*i][2]=real(A[i][1]); a[2*i][3]=-imag(A[i][1]);

a[2*i][4]=real(A[i][2]); a[2*i][5]=-imag(A[i][2]);

a[2*i][6]=real(A[i][3]); a[2*i][7]=-imag(A[i][3]);

a[2*i][8]=real(A[i][4]); a[2*i][9]=-imag(A[i][4]);

a[2*i][10]=real(A[i][5]); a[2*i][11]=-imag(A[i][5]);

a[2*i][12]=real(A[i][6]); a[2*i][13]=-imag(A[i][6]);

a[2*i][14]=real(A[i][7]); a[2*i][15]=-imag(A[i][7]);

f[2*i]=real(F[i]);

a[2*i+1][0]=imag(A[i][0]); a[2*i+1][1]=real(A[i][0]);

a[2*i+1][2]=imag(A[i][1]); a[2*i+1][3]=real(A[i][1]);

a[2*i+1][4]=imag(A[i][2]); a[2*i+1][5]=real(A[i][2]);

a[2*i+1][6]=imag(A[i][3]); a[2*i+1][7]=real(A[i][3]);