Основная информация Дата рождения: 28 Марта 1989 Пол: Мужской Семейное положение: женат Контакты VK: soleniyalex | |
Жизненная позиция Главное в жизни: саморазвитие Отн. к курению: резко негативное Отн. к алкоголю: резко негативное Вдохновляют: Любовь и наука | Личные интересы Интересы: спорт Любимые цитаты: ++++++++++++++++++++++++++++++++
float slau(float **a2, int n, float *b2, float *x)
{int i,j,k,r;
float max,c,s,m;
float **a, *b;
b=new float[n];
a=new float *[n];
for(i=0;i<n;i++)
{a[i]=new float[n];
b[i]=b2[i];}
for(i=0;i<n;i++)
for(j=0;j<n;j++)
a[i][j]=a2[i][j];
for(k=0;k<n;k++)
{ r=k;
max=fabs(a[k][k]);
for(i=k+1;i<n;i++)
if(fabs(a[i][k])>max)
{max=fabs(a[i][k]);
r=i;}
for(j=0;j<n;j++)
{c=a[k][j];
a[k][j]=a[r][j];
a[r][j]=c;}
c=b[k];
b[k]=b[r];
b[r]=c;
for(i=k+1;i<n;i++)
{for(m=a[i][k]/a[k][k],j=k;j<n;j++)
a[i][j]-=m*a[k][j];
b[i]-=m*b[k];}}
if(a[n-1][n-1]==0)
{if(b[n]==0)
return(-2);
else
return(-1);}
else
for(i=n-1;i>=0;i--)
{for(s=0,j=i+1;j<n;j++)
s+=a[i][j]*x[j];
x[i]=(b[i]-s)/a[i][i];}
return(0);} |