SOR Method를 사용해 연립방정식 근 구하기

소개글

SOR Method를 사용해 연립방정식 근 구하기에 대한 보고서 자료입니다.

목차

없음

본문내용

if(i==j+1)A[i][j]=-0.5;
}
}
for(i=0;i b[0]=0.5;
for(i=0;i for(i=0;i while(norminf!=0){
norminf=0;
for(i=0;i for(j=0;j if(j!=i)t+=A[i][j]*x[j];
}
x[i]=(1-w)*xp[i] + w/A[i][i]*(b[i]-t);
t=0;
}
for(i=0;i t=absol(x[i]-xp[i]);
if(t>norminf)norminf=t;
}
for(i=0;i iteration++;
}
printf("interation = %d\n",iteration);
for(i=0;i printf("P%d = %f ",i+1,x[i]);
if(i%5==4)printf("\n");
}
getchar();
}
float absol(float x)
{
if(x>=0)return x;
else return (-1) * x;
}
n=10
n=50
n=100
/*
SOR method for 7.5.4
*/
#include
#define n 9
float absol(float x);
main(){
int i,j,iteration=0;
float A[n][n],b[n],t=0;
float x[n],xp[n], norminf=1,w=1;
for(i=0;i for(j=0;j A[i][j]=0;
if(i==j)A[i][j]=1.0;
if(i==j-1)A[i][j]=1.0/3.0-1.0;
if(i==j+1)A[i][j]=-1.0/3.0;
}
}
for(i=0;i b[0]=0.5;
for(i=0;i for(i=0;i while(norminf!=0){
norminf=0;
for(i=0;i for(j=0;j if(j!=i)t+=A[i][j]*x[j];
}
x[i]=(1-w)*xp[i] + w/A[i][i]*(b[i]-t);
t=0;
}
for(i=0;i t=absol(x[i]-xp[i]);
if(t>norminf)norminf=t;
}
for(i=0;i iteration++;
}
printf("A[9][9]\n");
for(i=0;i for(j=0;j printf("%.3f ",A[i][j]);
}
printf("\n");
}
printf("interation = %d\n",iteration);
for(i=0;i printf("P%d = %f ",i+1,x[i]);
if(i%5==4)printf("\n");
}
getchar();
}
float absol(float x)
{
if(x>=0)return x;
else return (-1) * x;
}
n=10
n=50
n=100

키워드

SQR,   method,   연립방정식,   수치해석,   Gauss,   가우스