조 밑에 올린 복소수 짜기전에 DEV에선 이걸 짰었어요

DiesIrae

2011.10.09 21:05:51

조회 62 추천 0 댓글 5

#include <fstream>

#include <iostream>

#define L 3

#define N 3

int main()

{

using namespace std;

int i,j,k;

double **a, **b, **c, **d, **e, **f, **g, **A, **B, **h;

a=new double *[L+1];

b=new double *[L+1];

A=new double *[L+1];

B=new double *[L+1];

c=new double *[L+1];

d=new double *[L+1];

e=new double *[L+1];

f=new double *[L+1];

g=new double *[L+1];

h=new double *[L+1];

//행렬정의

for(i=0;i<=L;i++)

{

a[i]=new double [N+1];

b[i]=new double [N+1];

c[i]=new double [N+1];

d[i]=new double [N+1];

e[i]=new double [N+1];

f[i]=new double [N+1];

g[i]=new double [N+1];

A[i]=new double [N+1];

B[i]=new double [N+1];

h[i]=new double [N+1];

}

cout.setf(ios::fixed);

cout.precision(5);

a[1][1]=5; a[1][2]=-2; a[1][3]=3;

a[2][1]=4; a[2][2]=7; a[2][3]=-1;

a[3][1]=3; a[3][2]=4; a[3][3]=8;

A[1][1]=5; A[1][2]=-2; A[1][3]=3;

A[2][1]=4; A[2][2]=7; A[2][3]=-1;

A[3][1]=3; A[3][2]=4; A[3][3]=8;

cout << "Matrix A:" << endl;

for(i=1;i<=L;i++)

{

for(j=1;j<=L;j++)

cout << a[i][j] << "";

cout << endl;

}

//a 정의 및 출력

b[1][1]=-4; b[1][2]=8; b[1][3]=1;

b[2][1]=6; b[2][2]=-2; b[2][3]=7;

b[3][1]=4; b[3][2]=2; b[3][3]=5;

B[1][1]=-4; B[1][2]=8; B[1][3]=1;

B[2][1]=6; B[2][2]=-2; B[2][3]=7;

B[3][1]=4; B[3][2]=2; B[3][3]=5;

cout << endl << "Matrix B:" << endl;

for(i=1;i<=L;i++)

{

for(j=1;j<=L;j++)

cout << b[i][j] << "";

cout << endl;

}

//b 정의 및 출력

g[1][1]=1; g[1][2]=0; g[1][3]=0;

g[2][1]=0; g[2][2]=1; g[2][3]=0;

g[3][1]=0; g[3][2]=0; g[3][3]=1;

h[1][1]=1; h[1][2]=0; h[1][3]=0;

h[2][1]=0; h[2][2]=1; h[2][3]=0;

h[3][1]=0; h[3][2]=0; h[3][3]=1;

// I행렬을 g로 정의

cout << endl << "Matrix C(A multiplied by A):" << endl;

for(i=1;i<=L;i++)

{

for(j=1;j<=N;j++)

{

c[i][j]=0;

for(k=1;k<=L;k++)

c[i][j]+=a[i][k]*a[k][j];

cout << c[i][j] << "";

}

cout << endl;

}

//a행렬의 제곱 겅의

double Sum, m;

for(k=1;k<=N-1;k++)

for(i=k+1;i<=N;i++)

{

m=a[i][k]/a[k][k];

for(j=1;j<=N;j++)

{

a[i][j] -= m*a[k][j];

g[i][j] -= m*g[k][j];

}

for(i=N;i>=1;i--)

for(j=1;j<=N;j++)

{

Sum=0;

for(k=1;k<=N;k++)

Sum += a[i][k]*d[k][j];

d[i][j]=(g[i][j]-Sum)/a[i][i];

}

cout << endl << "Matrix d(inverse a):" << endl;

for(i=1;i<=L;i++)

{

for(j=1;j<=L;j++)

cout << d[i][j] << "";

cout << endl;

}

//a역행렬을 d로 정의

double Sum2, q;

for(k=1;k<=L-1;k++)

for(i=k+1;i<=L;i++)

{

q=b[i][k]/b[k][k];

for(j=1;j<=L;j++)

{

b[i][j] -= q*b[k][j];

h[i][j] -= q*h[k][j];

}

for(i=L;i>=1;i--)

for(j=1;j<=L;j++)

{

Sum2=0;

for(k=1;k<=L;k++)

Sum2 += b[i][k]*e[k][j];

e[i][j]=(h[i][j]-Sum2)/b[i][i];

}

cout << endl << "Matrix e(inverse b):" << endl;

for(i=1;i<=L;i++)

{

for(j=1;j<=L;j++)

cout << e[i][j] << "";

cout << endl;

}

//b역행렬을 e로정의

cout << endl << "Matrix F(total):" << endl;

for(i=1;i<=L;i++)

{

for(j=1;j<=N;j++)

{

f[i][j]=0;

for(k=1;k<=L;k++)

f[i][j]+=c[i][k]*e[k][j]+d[i][k]*B[k][j]-A[i][k]*B[k][j];

cout << f[i][j] << "";

}

cout << endl;

}

//토탈결과를 f로정의

for(i=1;i<=L;i++)

delete a[i];

delete b[i];

delete c[i];

delete d[i];

delete e[i];

delete f[i];

delete g[i];

delete A[i];

delete B[i];

delete h[i];

delete A, B, a, b, c, d, e, f, g, h;

system("pause");

return 0;

}

시팔 이렇게 긴게 돌아가는데 저 짧은게 왜 안돌아가 시발 고인 개같은 dev

고정닉 0

전체 댓글 0개

등록순 최신순 답글순

본문 보기

타인의 권리를 침해하거나 명예를 훼손하는 댓글은 운영원칙 및 관련 법률에 제재를 받을 수 있습니다.
Shift+Enter 키를 동시에 누르면 줄바꿈이 됩니다.

갤러리 리스트
번호	제목	글쓴이	작성일	조회	추천
설문	탈모 걱정 없어 보이는 머리숱 금수저 스타는?	운영자	25/07/14	-	-
공지	프로그래밍 갤러리 이용 안내 [88]	운영자	20.09.28	45619	65
2873043	여름철 내 몸냄새 꼭 확인해봐	ㅇㅇㅇㅇ(124.198)	10:34	3	0
2873041	뀨❤+	♥삣삐냥덩♥	10:20	9	0
2873040	원래 nl계 빨갱이들이 노무현 존나 싫어했는데 같은꼴 날 운명	♥삣삐냥덩♥	10:12	5	0
2873039	근본없는 2찢명 ㅋㅅㅋ 노무현 꼴 날듯	♥삣삐냥덩♥	10:10	8	0
2873038	역시 2찢명은 그냥 바지사장 이었네	♥삣삐냥덩♥	10:09	6	0
2873037	epicgames에서 문명6 무료 배포 중	발명도둑잡기	10:07	9	0
2873036	윤석열 대통령이 전쟁 막았던거네..	♥삣삐냥덩♥	10:06	12	0
2873034	나님은 사람을 좋아하지 않음	♥삣삐냥덩♥	09:42	15	0
2873032	원종이의 앱 ㅇㅅㅇ	따당(106.101)	09:36	25	0
2873031	게임개발자하고 싶은데	프갤러(106.101)	09:23	15	0
2873030	클로드 코드 미쳤다	노력	09:19	19	0
2873028	나무삼보	♥삣삐냥덩♥	08:59	13	0
2873027	❤✨☀⭐⚡☘♥+나님 시작합니당♥+☘⚡⭐☀✨❤	♥삣삐냥덩♥	08:40	14	0
2873025	고양이 ㅇㅅㅇ	헤르 미온느	08:28	18	0
2873024	태연 ㅇㅅㅇ	헤르 미온느	08:27	10	0
2873022	하루 한 번 헤르미온느 찬양	헤르 미온느	08:26	13	0
2873020	tiobe index 에서도 ADA 언어순위 떳더라 [3]	ㅇㅇ(183.101)	08:08	26	0
2873018	[논평] 방위비분담금 폐지하고 미군 주둔비 받아야	발명도둑잡기	08:07	13	0
2873015	루비 멋있어	ㅇㅇ(183.101)	08:01	21	0
2873013	나님 문제해결능력 ㄱㅆㅅㅌㅊ !!	♥삣삐냥덩♥	07:50	16	0
2873011	Ada 문법은 옛날의 PASCAL 지금의 DELPHI 가 느껴진다.	프갤러(59.16)	07:29	25	0
2873010	❤✨☀⭐⚡☘♥+나님 시작합니당♥+☘⚡⭐☀✨❤	♥삣삐냥덩♥	07:00	17	0
2873007	Ada 언어가 앞으로 유망할 예정 [1]	나르시	05:29	37	0
2873005	러빨러, c++빨러 왜 안 보이냐 [1]	나르시	05:10	35	0
2872997	난 프갤러들 대부분 서울 • 경기에 살줄 [2]	루도그담당(58.239)	01:58	65	0
2872990	KG에서 건수하나 올렸구나 [4]	개멍청한유라	01:38	57	0
2872989	전한길 표정 재밋지않냐? [7]	헬마스터	01:35	77	0
2872987	일하면서 틈틈히 개인용 프로그램만들었는데	ㅇㅇ(14.5)	01:05	35	0
2872985	모스탄이 사짜 냄새 나는 이유 2가지.JPG [1]	야옹아저씨	01:04	51	0
2872984	컴퓨터 프로그래밍 강의 4과목 양도 KG 아이티뱅크 [2]	프갤러(14.63)	00:39	48	1
2872983	아반떼CN7 페리 차량 5일 주행후기 [5]	야옹아저씨	00:28	95	7
2872981	유럽 복지국가 흑인 범죄율, 미국 흑인 범죄율 과학적 비교와 분석	발명도둑잡기	00:25	20	0
2872978	"너는 가장 많은 시간을 보내는 5명의 평균"의 문제 [1]	발명도둑잡기	00:19	43	0
2872976	비와서 그런가 인터넷이 간헐적으로 끊기네 [3]	ㅆㅇㅆ(124.216)	00:03	66	0
2872975	오늘 가족, 친척들 대인관계 상담을 잘 해줬더니 [3]	발명도둑잡기	07.19	52	0
2872974	ㅆㅇㅆ 클러스터랑 이중화차이좀	프갤러(106.101)	07.19	38	0
2872971	음기 충전	발명도둑잡기	07.19	23	0
2872969	천국 사이카 [1]	넥도리아(220.74)	07.19	35	0
2872968	우리신입뽑는대 너무귀찮타... [2]	프갤러(106.101)	07.19	79	0
2872966	돈까스시킴	프갤러(121.139)	07.19	21	0
2872964	코레일에서 이직할만한 기업 추천좀	프갤러(211.203)	07.19	23	0
2872963	이재명 정부가 부활시킨 특활비... '내로남불' 논쟁 끝내는 방법	발명도둑잡기	07.19	20	0
2872962	취미로 진짜 언어 배우는 사람 있음? [1]	ㅇㅇ(106.241)	07.19	45	0
2872961	"스윙스는 여자화장실에 들어갈 수 있다" 글 보니 드는 생각	발명도둑잡기	07.19	24	0
2872960	블로그글 써야하는데 최근 너무 바쁘다	ㅆㅇㅆ(124.216)	07.19	21	0
2872959	취업할라믄 무해야대지 [1]	무관	07.19	88	0
2872958	헤롱헤롱	발명도둑잡기	07.19	20	0
2872957	류류야 그만해라, 내 동창들 다 결혼했다. 나만 안했다. [2]	넥도리아(220.74)	07.19	32	0
2872956	엄마가 내 방 쓰레기래... 랜다. 낮에 [7]	넥도리아(220.74)	07.19	48	0
뉴스	소프라노 황수미, 다채로운 음악 선보이는 '사운드트랙' 콘서트 기획	디시트렌드	07.18