|
문제 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])
…(생략)
X1 = 128.214 , X2 = 456.643, X3 = 58.2143
|
by EES
C++을 이용해 구한 해는 근사적으로 맞는 것으로 생각되어진다. 확인을 위해 EES로 계산해 보았더니 거의 비슷한 결과값을 얻을 수 있었다. 결과가 맞는지 알아보기 위해 C++과 EES 로 구한 값들을 원식에 대입해보았다.
EES는 로 나왔고
C++은 로 나왔다.
대체적으로 EES가 좀 더 정확한 값을 나타내었는데 자세한 원인은 모르겠지만 EES는 공학용 계산에 적합하도록 많은 정확도를 가지고 계산을 하기 때문이라 추정할 수 있겠다. 그래서 좀 더 정확한 값을 얻기 위해 C++에서 변수들을 DOUBLE에서 LONG DOUBLE로 선언했지만 결과에는 변화가 없었다. 반올림에 의한 오차 때문에 약간의 부정확한 값을 얻었다고 생각할 수있다.
인 쇄
바로가기
lu분해법.hwp
[Size :
32 Kbyte
]
장바구니
