검색어 입력폼

PtLU Factorization

등록일 2005.01.01 한글파일한컴오피스 (hwp) | 2페이지 | 가격 1,000원

소개글

PtLU Factorization (Decomposition) 수치해석 기법입니다. 과제하실때 활용하세요.

컴파일 실행환경

Visual C++ 6.0
1.5G 노트북

본문내용

/*
PtLU factorization.
*/
#include<stdio.h>
main(){
int i,j,k;
float M[3][3]={{0,2,3},{1,1,-1},{0,-1,1}};
float PT[3][3]={{1,0,0},{0,1,0},{0,0,1}};
float L[3][3]={{0,0,0},{0,0,0},{0,0,0}},t;
for(i=0;i<2;i++){
if(M[i][i]==0){
for(j=i+1;j<3;j++){
if(M[i][i]==0 && M[j][i]!=0){
for(k=0;k<3;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
for(k=0;k<3;k++){
t = PT[k][i];
PT[k][i] = PT[k][j];
PT[k][j] = t;
}
}
}
}
for(j=i+1;j<3;j++){
t=M[j][i]/M[i][i];
for(k=i;k<3;k++){
M[j][k]=M[j][k]-M[i][k]*t;
}
L[j][i]=t;
}
}
for(i=0;i<3;i++)L[i][i]=1;
printf("Pt Matrix\n");
for(j=0;j<3;j++){
for(k=0;k<3;k++){
printf("%3.3f ",PT[j][k]);
}
printf("\n");
}
printf("\nL Matrix\n");
for(j=0;j<3;j++){
for(k=0;k<3;k++){
printf("%3.3f ",L[j][k]);
}
printf("\n");
}
printf("\nU Matrix\n");
for(j=0;j<3;j++){
다운로드 맨위로