디시인사이드 갤러리

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

갤러리 본문 영역

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

DiesIrae갤로그로 이동합니다. 2011.10.09 21:05:51
조회 60 추천 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 - -
AD 휴대폰 액세서리 SALE 운영자 25/07/15 - -
280241 꼬꼬형은 몇년차야? [13] ^0^갤로그로 이동합니다. 11.10.12 109 0
280240 다들 퇴근하셨나?_? [2] Dawnwalkre갤로그로 이동합니다. 11.10.12 53 0
280239 컴공에서 배우는 것... [8] 534F444D61737465갤로그로 이동합니다. 11.10.12 167 0
280238 횽들 수고하심 [4] 어떡해갤로그로 이동합니다. 11.10.12 49 0
280237 우리 적어도 중앙대 이하급 대졸 뽑는편인데... [6] 꼬꼬월드갤로그로 이동합니다. 11.10.12 172 0
280236 C언어 할 줄 아는사람 이것 좀 봐줘 [3] 시개갤로그로 이동합니다. 11.10.12 114 0
280235 대학생으로서 프로그래밍언어 배운다는건 수치아니냐? [4] ㅇㅇㅇㅇ(175.204) 11.10.12 127 0
280234 집에[ 전기끊겨서 피ㅣㅆ방와서 프갤중 [1] ^0^갤로그로 이동합니다. 11.10.12 36 0
280233 이야 시발 슬레이트7대박이다 [3] 꿀레(14.33) 11.10.12 234 0
280232 근데 컴공과에서는 뭐 가르침? [23] 꼬꼬월드갤로그로 이동합니다. 11.10.12 251 0
280231 난 그래도 벌레 덕분에 즐겁다. [1] Ashenvale(220.90) 11.10.12 71 0
280230 일광면 형 글땜에 상처받았다... [3] 얼빠진호랑이갤로그로 이동합니다. 11.10.12 107 0
280229 대학교 문과 졸업해서 학원도 안다니고 개발자로 취직한 사람도 있음.. [6] 꼬꼬월드갤로그로 이동합니다. 11.10.12 177 0
280228 C언어에서 받은 충격2 [9] 도루도라갤로그로 이동합니다. 11.10.12 125 0
280226 Spring 프레임워크 언제까지 써먹을 수 있을까 꼬꼬월드갤로그로 이동합니다. 11.10.12 41 0
280225 학원생으로써 학원 후기 [35] 꼬꼬월드갤로그로 이동합니다. 11.10.12 506 0
280223 정말 시골은 정글인가? ㅎㅎㅎ [11] 하아(110.13) 11.10.12 108 0
280222 씨언어 공부하는데 진짜 완전 개골때린다 [13] 도루도라갤로그로 이동합니다. 11.10.12 106 0
280221 아까 숫자 두개입력 해결했는데요. 궁금한게 생겼어요. [2] 하양(61.73) 11.10.12 54 0
280220 kes오늘 개막했네 [1] 꿀레(14.33) 11.10.12 47 0
280219 솔직히 문과는 문과 공부만하세요 [9] 일광면(119.198) 11.10.12 148 0
280217 2의 n승 판별법 [16] ㄹㄴㅁ(183.98) 11.10.12 199 0
280216 프갤러는 야구 안봄???? [2] 없쥐갤로그로 이동합니다. 11.10.12 46 0
280215 얼빠진호랑이 언니야는 인증 했나? [6] iljeomobolt갤로그로 이동합니다. 11.10.12 95 0
280214 락은 저항이니까 다운 받아들으라구??ㅡㅡ [2] abbey road갤로그로 이동합니다. 11.10.12 62 0
280213 얘들아 고맙다 너희들 덕분에 해결했엌ㅋㅋㅋ [5] 없쥐갤로그로 이동합니다. 11.10.12 62 0
280212 웹프로그램 짜려면 보안도 같이 알아야 되는거야? [5] 아자(211.237) 11.10.12 91 0
280211 국카스텐 검정치마 못을 세계로 진출시키자 [3] abbey road갤로그로 이동합니다. 11.10.12 74 0
280210 프로그래밍했는데 질문 좀 드릴께요~ [1] 시개갤로그로 이동합니다. 11.10.12 40 0
280209 오페라 터보쓰미 느려서 안되것다 534F444D61737465갤로그로 이동합니다. 11.10.12 40 0
280208 우앙 갤로그 조회수 포텐터졌네 ㅋ_ㅋ [3] Dawnwalkre갤로그로 이동합니다. 11.10.12 86 0
280206 없쥐횽 봅니다 [4] 일광면(119.198) 11.10.12 57 0
280205 오늘도 하이킥 안할라나??? [4] 거칠게갤로그로 이동합니다. 11.10.12 70 0
280204 아까 올렸는데 웃기만 하고 자바 알려주세욬ㅋㅋ [7] 없쥐갤로그로 이동합니다. 11.10.12 62 0
280203 우왕 시발 new gay[max]횽이다 [4] 일광면(119.198) 11.10.12 73 0
280202 야 이 바보들아 아직도 패턴을 모르겠냐? [1] 534F444D61737465갤로그로 이동합니다. 11.10.12 69 0
280201 누구 여기 학원다니는 사람 없어? [12] 얼빠진호랑이갤로그로 이동합니다. 11.10.12 199 0
280200 ㅋㅋㅋㅋㅋㅋㅋㅋ min-MAX 알고리즘 다짰다 [1] McHello갤로그로 이동합니다. 11.10.12 75 0
280199 ㅋㅋㅋㅋㅋ 아웃겨 벌레도 여자한텐 안되는구나 ㅋㅋㅋㅋㅋㅋ [3] blackd갤로그로 이동합니다. 11.10.12 103 0
280198 으아아아 이상한 떡밥 물어가지고.. 534F444D61737465갤로그로 이동합니다. 11.10.12 43 0
280196 곽노현팀장형 봐라 [8] 얼빠진호랑이갤로그로 이동합니다. 11.10.12 129 0
280195 아오 500메가 요금제인데 하루만에 70메가 썼네 534F444D61737465갤로그로 이동합니다. 11.10.12 38 0
280194 얼호는 게이가 아니다 [1] Adelposs갤로그로 이동합니다. 11.10.12 49 0
280192 그럼 이제 프갤 공식 동성애자가 나왔구나 [3] 534F444D61737465갤로그로 이동합니다. 11.10.12 50 0
280191 ㅋㅋㅋㅋㅋ 이런 이야기에서 아스가 빠지다니 [1] Adelposs갤로그로 이동합니다. 11.10.12 58 0
280189 아스형 그럼 이런 짤 좋아할라나? [2] 얼빠진호랑이갤로그로 이동합니다. 11.10.12 105 0
280188 칼퇴근시켜준다고 날 현혹시켜놓고~~~~~~~~` [2] 꼬꼬월드갤로그로 이동합니다. 11.10.12 79 0
280187 근데 레즈는 왜 생기는걸까 [6] 534F444D61737465갤로그로 이동합니다. 11.10.12 83 0
280185 나도 사실 의느님에게 내 엉덩일 한번 허락한적있었어 [1] Adelposs갤로그로 이동합니다. 11.10.12 52 0
280184 왠지 길 걷다보면 프갤러 있을거 같다 [6] 로하로하알로하갤로그로 이동합니다. 11.10.12 68 0
뉴스 강철부대 멤버 ‘커밍아웃’… 동성연인과 럽스타그램 시작 디시트렌드 10:00
갤러리 내부 검색
제목+내용게시물 정렬 옵션

오른쪽 컨텐츠 영역

실시간 베스트

1/8

뉴스

디시미디어

디시이슈

1/2