JACOBI ITERATION program code © Khec

#include
#include
#include
#define E 0.0001
void main()
{clrscr();
int n,i,j,key;
float a[10][10],b[10],xo[10],x[10],sum;
printf("enter n");
scanf("%d",&n);
printf("enter the matrix A");
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
printf("a[%d][%d]",i,j);
scanf("%f",&a[i][j]);
}
printf("b[%d]",i);
scanf("%f",&b[i]);
xo[i]=b[i]/a[i][i];
}
l:
key=0;
for(i=1;i<=n;i++)
{
sum=b[i];
for(j=1;j<=n;j++)
{if(j!=i)
{sum=sum-a[i][j]*xo[j];
}
}
x[i]=sum/a[i][i];

if(key==0)
{
  if(fabs((x[i]-xo[i])/x[i])>E)
{key=1;
}
 }
}
if(key==1)
{
 for(i=1;i<=n;i++)
 {
xo[i]=x[i];
}
goto l;
}
for(i=0;i<=n;i++)
{
printf("x[%d]:%f",i,x[i]);
}
getch();
}