디시인사이드 갤러리

갤러리 이슈박스, 최근방문 갤러리

갤러리 본문 영역

조 밑에 올린 복소수 짜기전에 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

0

댓글 영역

전체 댓글 0
본문 보기

하단 갤러리 리스트 영역

왼쪽 컨텐츠 영역

갤러리 리스트 영역

갤러리 리스트
번호 제목 글쓴이 작성일 조회 추천
설문 끝까지 다 본 걸 후회하게 만든 용두사미 드라마는? 운영자 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
갤러리 내부 검색
제목+내용게시물 정렬 옵션

오른쪽 컨텐츠 영역

실시간 베스트

1/8

뉴스

디시미디어

디시이슈

1/2