디시인사이드 갤러리

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

갤러리 본문 영역

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

DiesIrae갤로그로 이동합니다. 2011.10.09 21:05:51
조회 63 추천 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/14 - -
277118 빨간머리 좋아하는 사람 있어? (살짝 은꼬르) [3] blackd갤로그로 이동합니다. 11.09.30 155 0
277117 프로그램짜봤는데요 딱 한부분이 이해가 안가는데 도와주세요 ㅠ 푸훼훼휄(210.101) 11.09.30 56 0
277116 블디형은 봐줍니다 [15] 三didas갤로그로 이동합니다. 11.09.30 92 0
277115 병특 자바개발쪽가면 경럭인정해줍니까? [1] 一ㅡ갤로그로 이동합니다. 11.09.30 93 0
277114 집안 전체가 공대 출신이라 [8] 일광면(119.198) 11.09.30 154 0
277113 헐.. 나 뜬금없이 진급했어.. [12] McHello갤로그로 이동합니다. 11.09.30 175 0
277112 이번주에 한것 [10] +어게인갤로그로 이동합니다. 11.09.30 154 0
277111 황금연휴라니? 뭔날이냐? [2] 一ㅡ갤로그로 이동합니다. 11.09.30 71 0
277110 좋은 아버지 는 봅니다. [2] 좋은아버지갤로그로 이동합니다. 11.09.30 62 0
277109 형들은 머리식힐때 무슨게임해?ㅇ [3] 三didas갤로그로 이동합니다. 11.09.30 82 0
277108 centos nvidia 드라이버 삽질 성공 [2] 이모군(58.228) 11.09.30 68 0
277107 vs2010에서 for( int R1 ; ; )안되나여 [1] 이힝(121.177) 11.09.30 56 0
277106 우리도 돈 좀 만져보자 [6] 막장갤신학생갤로그로 이동합니다. 11.09.30 86 0
277105 짬트랏슈 형은 봅니다. [27] 좋은아버지갤로그로 이동합니다. 11.09.30 197 0
277104 언니들 리눅스에서 getch구현(화면출력없이 바로 입력) 평가 좀 [5] dot(125.128) 11.09.30 89 0
277103 입갤과동시에 나의 랩실생활 [2] 三didas갤로그로 이동합니다. 11.09.30 86 0
277102 나도 나으 랩실 생활... [1] Kanon갤로그로 이동합니다. 11.09.30 247 0
277101 C++ 답변좀 해주세요 ㅜㅜ [7] ㅇㅇ(183.106) 11.09.30 95 0
277100 횽들 내 후배가 뭔가 만들었는데 돈이 될까? [34] 생물학(125.248) 11.09.30 229 0
277098 연봉 2200이라고치면요 [21] 짬트라슈갤로그로 이동합니다. 11.09.30 346 0
277097 내친구의 대학교 연구실 생활 들려줄게. [3] 좋은아버지갤로그로 이동합니다. 11.09.30 194 0
277096 대학교 연구실 원래 이렇게 줫같음?? [27] 므느으르(210.107) 11.09.30 283 0
277095 그런데 요즘 이도경은 뭐하고 지내는지 알아 ? [2] 좋은아버지갤로그로 이동합니다. 11.09.30 99 0
277094 네이버캐스트 - 인류의 삶을 바꾼 위대한 기계 컴퓨터 궁금이(121.129) 11.09.30 46 0
277093 오늘 좆밥인 내 집중력에 놀랏다 [3] 로하로하알로하갤로그로 이동합니다. 11.09.30 97 0
277092 앱등은 앱등앱등해 [6] monoless갤로그로 이동합니다. 11.09.30 88 0
277091 그런데 밑에 기사 보니까 이번에 MS + 삼성 이자나? [2] 좋은아버지갤로그로 이동합니다. 11.09.30 72 0
277090 오늘은 책만 보다 끝나는구나. [2] 좋은아버지갤로그로 이동합니다. 11.09.30 73 0
277089 황금 연휴 시작이네 [3] monoless갤로그로 이동합니다. 11.09.30 75 0
277088 아스트라로테 첨 봤을때 부터 떠오른건데.. [3] blackd갤로그로 이동합니다. 11.09.30 114 0
277087 아 잉여로운 휴일이 되겠군ㅋ 나리링갤로그로 이동합니다. 11.09.30 26 0
277086 엉엉 밑에 for문 예제 소스 올려주다가 도배로 차단 당했어요 엉엉 [3] 수상한김토끼갤로그로 이동합니다. 11.09.30 86 0
277084 리눅스 좀 하는 형 만 .. [10] 뿔테갤로그로 이동합니다. 11.09.30 139 0
277083 와 오라클에 dual 이라는 테이블 있는데 장난 아니넹 [1] 좋은아버지갤로그로 이동합니다. 11.09.30 124 0
277082 win+D가 바탕화면보기인데 명령어같은건 없나요? 이힝(121.177) 11.09.30 54 0
277081 한 내년쯤 내 비지니스 사이트를 돌릴까 하는데 [1] 쵸쵸(112.219) 11.09.30 45 0
277076 Rendering 3D with CSS and JavaScript [3] 김토끼(203.244) 11.09.30 75 0
277074 형들 간단한 For문 예젠데 도와줘. ㅜ.ㅜ;; [14] Chic.(211.223) 11.09.30 153 0
277073 모노폴리 게임 AI 프로그램을 리서치 프로젝트로 잡을까 생각중 [1] 멸자(130.126) 11.09.30 75 0
277069 자꾸 신경쓰이니까 은근 짜증나네... [2] 천재해커(59.14) 11.09.30 104 0
277068 야호 2d게임 만들기 막겜개발 3탄 포스팅완료 McHello갤로그로 이동합니다. 11.09.30 94 0
277067 샘숭&마소 크로스! [2] ,.,.(222.121) 11.09.30 70 0
277065 임베디드XP에서 시리얼포트 잡고있는 기본프로세스 있나요? (조공 추가) 퉁퉁퉁갤로그로 이동합니다. 11.09.30 152 0
277063 오늘같은 날은 [1] monoless갤로그로 이동합니다. 11.09.30 34 0
277062 나도 오늘부터 밀당을 좀 해볼까 한다 [4] 분당살람갤로그로 이동합니다. 11.09.30 92 0
277061 형들 ID3DXMesh쓰는게 좋아? [2] 늅늅(115.23) 11.09.30 107 0
277060 근데 매일 7만원 씩 쓴다고 하면 [8] 좋은아버지갤로그로 이동합니다. 11.09.30 123 0
277059 미국과 한국의 차이가 뭘까???? [28] ㅠㅠ(169.229) 11.09.30 425 0
277058 n드라이브써봣냐? [3] ㅇㅇ(175.214) 11.09.30 120 0
277057 와~ 이기사 댓글수 쩐다... [4] blackd갤로그로 이동합니다. 11.09.30 147 0
뉴스 '더 시즌즈-박보검의 칸타빌레' 슈퍼주니어·선우정아·폴 블랑코·지소쿠리클럽, 한여름 밤에 떠나는 시간여행 디시트렌드 07.19
갤러리 내부 검색
제목+내용게시물 정렬 옵션

오른쪽 컨텐츠 영역

실시간 베스트

1/8

뉴스

디시미디어

디시이슈

1/2