Excelente programación
Jacoby
# incluir ltiostream gt
Usar espacio de nombres estándar
# incluir ltmath.h gt
# Definición N 3
int main()
{
int i, j;
Doble a[N][ N]= {
{10,-1,-2},
{-1,10,-2},
{-1,- 1, 5 }};
doble b[N]={7.2, 8.3, 4.2};
doble x[N]={0}, TM;
doble xnew[N]= { 0 };
Doble máximo
int m = 0;
Y ( m)
{
cout lt ltm;
max = 0;
for(I = 0; i ltn; i )
{ p>
TM = 0;
for(j = 0; j ltn; j )
si (i! = j)
TM = x[j]* a[I][j];
xnew[I]=(b[I]-TM)/a[I][I];
cout lt lt" x[" lt; ltI lt lt"]= " lt; ltxnew[I];
}
cout lt ltendl
for(j = 0; j ltn; j )
{
if(fabs(x[j]-xnew[j]) gt; máximo)
max = fabs(x[j]-xnew[j]);
x[j]= xnew[j];
}
if(max lt; 0.000001)
Romper;
}
Regresar 1;
}
古西
# include ltiostream gt
Usar espacio de nombres std
# include ltmath.h gt
#Definición N 3
int main( )
{
int i, j;
Doble a[N][N]={
{10,-1 ,-2},
{-1, 10,-2},
{-1,-1,5}};
doble b[ N]={7.2, 8.3, 4.2};
doble x[N]={0}, TM;
int m = 0;
p >
Doble máximo
Y ( m)
{
cout lt ltm;
max = 0;
for(I = 0;iltn;i)
{
TM = 0;
for(j = 0;j ltn;j )
Si (yo! =j)
TM = x[j]* a[I][j];
if(fabs((b[I]-TM)/a[I][ I]-x[I]) gt;max)
max = fabs((b[I]-TM)/a[I][I]-x[I]);
x[yo]=(b[yo]
-TM)/a[I][I];
cout lt lt" x[" lt; lti lt lt"]= " ltx[I];
}
cout lt ltendl
if(max lt; 0.000001)
break;
}
return 1 ;
}