|
| |
|
|
수치해석 lu분해법 c언어 코드
파일 
lu분해법.hwp
[Size :
32 Kbyte
]
lu분해법.hwp
[Size :
32 Kbyte
]
분량 5 Page
가격 1,000 원
×
링크주소가 복사 되었습니다.
원하는 곳에 붙혀넣기 하세요
원하는 곳에 붙혀넣기 하세요
카트
|
다운받기
|
 네이버 ID로 다운 받기  카카오 ID로 다운 받기  구글 ID로 다운 받기  페이스북 ID로 다운 받기 |
뒤로
|
본문/내용
문제 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"<>a[i][j];
}
}
cout<>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;
}
}
dummy=o[p];
o[p]=o[k];
o…
…
…
| 📝 Regist Info |
|
I D : kwg3** Date : 2011-11-21 FileNo : 11032539 |
Cart  | |
| |