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

DiesIrae

2011.10.09 21:05:51

조회 55 추천 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/07	-	-
공지	프로그래밍 갤러리 이용 안내 [88]	운영자	20.09.28	45293	65
2870653	나팔 ㅇㅅㅇ	헤르 미온느	04:24	8	0
2870652	태연 ㅇㅅㅇ	헤르 미온느	03:44	11	0
2870651	하루 한 번 헤르미온느 찬양	헤르 미온느	03:41	12	0
2870649	백신맞고 몸 망가진 2찍들은 국짐당 원망해야지 [5]	야옹아저씨	02:42	33	5
2870648	1일 1딸 멈추라고 하시면 멈추겠습니다.	도리스아(121.139)	02:32	12	0
2870647	나는 기후위기를 걱정한 적이 없었는데 [8]	아스카영원히사랑해	01:34	41	0
2870646	MCP 에이전트로 코딩하면서 나는 프로그래밍 책 읽고 있는데	ㅆㅇㅆ(124.216)	01:24	21	0
2870644	그러고보니 예전 프갤러들 싹 사라졌네 [1]	박민준	00:48	63	1
2870643	울산에서 일하는 31살 공장충인데 지금 입문 어떻게 생각해 [3]	ㅇㅇ(118.235)	00:25	57	0
2870642	바탕화면에 공룡 키워 보세요....	프갤러(118.41)	07.08	58	1
2870641	프갤하는거 들켰다	프갤러(220.65)	07.08	29	0
2870639	성추행해서 미안해 너무 아파서 그랬어...	넥도리아(223.38)	07.08	25	0
2870638	고1 프로그래밍에 손을 대보려합니다 [8]	프갤러(218.239)	07.08	87	0
2870636	가끔 나 따라 다니는 유동 저거 냥덩이 아닐까 한다 [1]	ㅆㅇㅆ(124.216)	07.08	30	0
2870635	조루디 ㅆㅇㅆ 잘어울려	ㅇㅇ(211.235)	07.08	32	2
2870634	블리치 정주행 예정 [1]	♥냥덩이♥	07.08	20	0
2870633	R 뭐 배우는데 OLD R베이스 문법? [2]	ㅆㅇㅆ(124.216)	07.08	41	0
2870632	결혼하면 프갤러들 축의금 줘 [3]	조루디(110.35)	07.08	44	0
2870631	R은 대학원 갈 거 아니면 뭐 쓸 일 거의 없을텐데	ㅆㅇㅆ(124.216)	07.08	25	0
2870630	제발 한숨쉬지마세요 노처녀드라 [1]	조루디(110.35)	07.08	41	0
2870629	문관데 제발 도와주세요 제발 제발제발 [15]	프갤러(118.47)	07.08	80	0
2870627	주식대박났다 [2]	조루디(110.35)	07.08	64	0
2870625	인공지능 슬기❤ [1]	♥냥덩이♥	07.08	56	0
2870624	면접보면서 당황했던 경험	ㅇㅇ(112.169)	07.08	26	0
2870622	17일 도쿄 여행인데	아스카영원히사랑해	07.08	42	0
2870620	PC-6001 페이지용 Contiki	발명도둑잡기(118.216)	07.08	18	0
2870618	국비들으면 쉽게 취업하는시대일때 국비 짤린거면 얼마나 병신새끼인거임??	뒷.통수한방(1.213)	07.08	61	1
2870617	좆같은 기술스택을 쓰고있다 <- 탈출해야함	프갤러(59.14)	07.08	40	0
2870615	어우 일나가기 실타..ㅇㅅㅇ	헤르 미온느	07.08	15	0
2870613	러스트 비판서 거의 다 썼다.	루비	07.08	24	0
2870610	냥덩이는 그냥 귀척하는 윤서인이라봐야 [2]	류류(115.136)	07.08	35	0
2870609	윤석열 무능한건 맞지않나 ㅇㅅㅇ	류류(115.136)	07.08	17	0
2870607	멍퀴님 제발 부탁드립니다 그만해주세요 [1]	♥냥덩이♥	07.08	40	0
2870605	실력상승은 인정에서 나온당⭐+ By 나님 [1]	♥냥덩이♥	07.08	37	0
2870603	별찍기를 무시하지마라 [1]	배구공(119.202)	07.08	27	0
2870602	이 사람이 대다나당⭐ [1]	♥냥덩이♥	07.08	41	0
2870600	❤✨☀⭐나님 시작합니당⭐☀✨❤	♥냥덩이♥	07.08	32	0
2870599	요즘 날이더워 그런강 피곤해서 일찍 자게됨 [5]	♥냥덩이♥	07.08	45	0
2870597	방화벽 지금 2주째 안 뚫어주는데 [12]	아스카영원히사랑해	07.08	102	0
2870595	퇴근하자 [8]	개멍청한유라	07.08	46	0
2870594	나는조현병이야 나는내향적이야 [5]	손발이시립디다	07.08	49	0
2870592	나이처먹고 피지컬은 떨어지는데 경험만 쌓여서 이거 저거 한번에 생각 [2]	프갤러(218.154)	07.08	60	0
2870590	프로펙트인가 지원해본 사람 있음? [1]	ㅇㅇ(118.36)	07.08	31	0
2870589	정좌불능증	재현	07.08	22	0
2870587	그냥 나가긴 싫고 실업급여라도.. [5]	ㅇㅇ(211.235)	07.08	39	0
2870585	와 근데 config 설정이 80여개 넘어가니까 수정하다 지치네 [4]	ㅆㅇㅆ(124.216)	07.08	64	0
2870584	deepseek 정말 못쓰겠다 쓰고싶어도 [1]	뒷통수한방(1.213)	07.08	42	0
2870583	와 근데 거북목 존나 많구나 [2]	루도그담당(211.184)	07.08	63	0
2870582	파이썬 똥쓰레기네	ㅇㅇ(106.101)	07.08	63	0
뉴스	‘나솔’ 22기 영숙 “결혼 2주 만에 이혼…‘이것’ 때문에 시부모 반대”	디시트렌드	07.08