레포트 > 공학기술

수치해석 lu분해법 c언어 코드

등록일 : 2011-11-21
갱신일 : 2011-11-21


  lu분해법.hwp   [size : 32 Kbyte]
  1,000   5 Page
 
  67%
 
수치해석 lu분해법 c언어 코드
문제 10.6 다음과 같은 연립 방정식의 해를 부분피벗화를 사용한 LU 분해법으로 구하라. (결과를 원식에 대입해 맞는지 확인할 것)





#include
#include
using namespace std;
double ABS(double t)
{
if(t<0)
return t=-t;
else return t;
}


int main()
{
int i,ii,j,k,p,er=0,st=0,stt=0;
double a[3][3]={0,},b[3]={0,},s[3]={0,},x[3]={0,},o[3]={0,};
double big,dummy,tol=0.0001,factor,sum=0;
for(i=0;i<=2;i++)
{
for(j=0;j<=2;j++)
{
cout<<"벡터 A"< cin>>a[i][j];
}
}
cout< for(i=0;i<=2;i++)
{
cout<<"벡터 B"< cin>>b[i];
}
//decompose
for(i=0;i<=2;i++)
{
o[i]=i;
s[i]=ABS(a[i][0]);
for(j=1;j<=2;j++)
{
if(ABS(a[i][j])>s[i])
s[i]=ABS(a[i][j]);
}
}
for(k=0;k<=1;k++)
{
//pivot
p=k;
st=o[k];
big=ABS(a[st][k]/s[st]);
for(ii=(k+1);ii<=2;ii++)
{
st=o[ii];
dummy=ABS(a[st][k]/s[st]);
if(dummy>big)
{
big=dummy;
p=ii;
}
}
dumm…(생략)

∴Tip Menu

lu분해법   수치해석   c언어   코드