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

v[i][j]=a[i][j];

w[i]=f[i];

}

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

{

maxk=k;

max=fabs(a[k][k]);

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

if(fabs(a[i][k])>max)

{

maxk=i;

max=fabs(a[i][k]);

}

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

{

e=a[k][i];

a[k][i]=a[maxk][i];

a[maxk][i]=e;

}

e=f[k];

f[k]=f[maxk];

f[maxk]=e;

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

{

c=a[i][k]/a[k][k];

f[i]=f[i]-f[k]*c;

for(j=k;j<N;j++)

a[i][j]=a[i][j]-a[k][j]*c;

}