디시인사이드 갤러리

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

갤러리 본문 영역

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

DiesIrae갤로그로 이동합니다. 2011.10.09 21:05:51
조회 59 추천 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 디지털 액세서리 기간한정 세일! 운영자 25/07/11 - -
공지 프로그래밍 갤러리 이용 안내 [88] 운영자 20.09.28 45474 65
2871779 헬창들만 아는 고민 하… ㅇㅇㅇㅇ(222.105) 06:52 4 0
2871778 세상에 억까가 심해더 시련을 견디고, 시련 넥도리아(220.74) 06:25 10 0
2871777 너네 개발공부 얼마나해 [4] 프갤러(106.101) 06:18 38 0
2871776 ❤✨☀⭐⚡☘♥+나님 시작합니당♥+☘⚡⭐☀✨❤ ♥지나가던길냥덩♥갤로그로 이동합니다. 06:08 11 0
2871775 저녁 안 먹구 간식 먹구 잣더니 머근개 없어서 그런지 끙야가 잘 안나오네 ♥지나가던길냥덩♥갤로그로 이동합니다. 05:54 11 0
2871774 무제 ㅇㅅㅇ 헤르 미온느갤로그로 이동합니다. 05:53 12 0
2871773 비스크돌 2기 ost 는 별루넹 ♥지나가던길냥덩♥갤로그로 이동합니다. 05:51 9 0
2871772 태연 ㅇㅅㅇ 헤르 미온느갤로그로 이동합니다. 05:51 7 0
2871771 하루 한 번 헤르미온느 찬양 헤르 미온느갤로그로 이동합니다. 05:50 16 0
2871770 C++ 인기는 왜 죽지 않는가? [3] 프갤러(59.16) 05:38 22 0
2871769 나님 모닝 미라클⭐+ ♥지나가던길냥덩♥갤로그로 이동합니다. 05:18 11 0
2871768 ❤✨☀⭐⚡☘♥+나님 시작합니당♥+☘⚡⭐☀✨❤ ♥지나가던길냥덩♥갤로그로 이동합니다. 05:17 9 0
2871767 프갤 저지능자 보수들은 뭔 찾아볼 생각을 못하냐.. ㅆㅇㅆ(124.216) 04:23 23 0
2871766 어느 디씨인의 진심어린 고백 ㅇㅇㅇㅇ(49.254) 03:30 19 0
2871765 자연어로 죄다 코딩해버리니까 현타오네 [1] 프갤러(175.193) 02:34 60 0
2871763 한국 주식 8~90%는 상위 10% 부자가 갖고있다 발명도둑잡기갤로그로 이동합니다. 01:26 18 0
2871761 부자감세 패싱하고 초부자 감세? 배당소득 분리과세의 함정 발명도둑잡기갤로그로 이동합니다. 01:18 13 0
2871759 하 취업 [3] 프갤러(124.50) 01:09 52 0
2871755 노조가 없으니까 당연히 힘이 안 실리지. 현대사회서 노조가 없으면 [2] ㅆㅇㅆ(124.216) 00:51 39 2
2871753 해외개발자들 세금써서 들여온다는거 ㄹㅇ이가.. [5] ㅇㅇ(223.38) 00:43 62 5
2871752 디시애들 멍청함 보다보면 혀를 내두르지 않냐? 노벨 경제학상을 전면부정함 [3] ㅆㅇㅆ(124.216) 00:42 53 0
2871750 잘자요 여러분 [5] 개멍청한유라갤로그로 이동합니다. 00:27 43 0
2871748 국내 자동차 업계는 망하게되있음 [5] ㅇㅇ(211.235) 00:20 83 0
2871747 백준 ㅈㄴ쉬운 문제 푸는데 3시간 걸리면 그냥 접는게 맞음? [1] 프갤러(1.238) 00:20 40 0
2871744 전에 ㅈㄹㄷ가 에고라는 적 책 추천해줬던 것 같은데 hrin(220.120) 00:07 20 0
2871742 삼국지 성씨로 보면 재밌는게 [2] 프갤러(211.210) 07.14 36 0
2871740 다들 임베를 안하려고 하니까 [4] CANON갤로그로 이동합니다. 07.14 77 0
2871736 프론트개발잔데 백엔드로 이직하면 연봉 [4] 프갤러(58.227) 07.14 86 0
2871734 중소 it기업 신입은 평균 학벌이 어캐됨? 프갤러(183.109) 07.14 41 0
2871733 유니티 노무 재밌다 [3] 루도그담당(58.239) 07.14 50 0
2871732 동양에서 제일 오래된 노래책 시경에 "고양이는 귀한 신"이라고 적힘 [1] 발명도둑잡기갤로그로 이동합니다. 07.14 20 0
2871730 차분하게 생각해보면 인생에 무서울것도 없는 것 같아. [1] ㅇㅇ(223.38) 07.14 23 0
2871728 시사힙합 내란의 힘 발명도둑잡기갤로그로 이동합니다. 07.14 11 0
2871726 개인프로젝트 배포 클라우드로 할때 가성비 조합 추천좀 [2] 프갤러(59.26) 07.14 37 0
2871724 좃같네 퉤 [3] 조루디(1.222) 07.14 45 0
2871723 파인만이나 아인슈타인은 배우는게 재밌었을까?? [3] ㅇㅇ(223.38) 07.14 28 0
2871721 조선시대에 있던 의외의 병과 발명도둑잡기갤로그로 이동합니다. 07.14 14 0
2871720 요즘은 시바 ㅋㅋ ai한태 블로그 써달라고하네 [1] 프갤러(61.79) 07.14 43 0
2871718 친중좌파 ♥지나가던길냥덩♥갤로그로 이동합니다. 07.14 20 0
2871717 나님 떠올랐당 이미지가 풍경이 ♥지나가던길냥덩♥갤로그로 이동합니다. 07.14 21 0
2871715 서울 가서 살면 빚이랑 방세때문에 300스타트 아니면 못모아 [7] ㅆㅇㅆ(124.216) 07.14 55 1
2871714 아 이분 성대 교수되셨네 [1] 아스카영원히사랑해갤로그로 이동합니다. 07.14 39 0
2871713 섹스도 통장잔고도 없다. 나에겐 빚뿐 [9] ㅆㅇㅆ(124.216) 07.14 78 0
2871711 ㅆㅇㅆ 몇살이냐? ㅋ [1] 아스카영원히사랑해갤로그로 이동합니다. 07.14 48 0
2871710 꺼억 재현갤로그로 이동합니다. 07.14 23 0
2871708 나 아스카 목요일 도쿄 여행을 기점으로 [7] 아스카영원히사랑해갤로그로 이동합니다. 07.14 57 0
2871706 저 아랫놈 말대로 하려면 컴포넌트 자체를 제네릭 화시켜야함 [1] ㅆㅇㅆ(124.216) 07.14 34 0
2871705 유니티중인데 문제가 생김 루도그담당(58.239) 07.14 29 0
2871703 ai도 결국 지가 아는만큼 쓰는거야 [5] 프갤러(110.8) 07.14 48 0
뉴스 '우리는 잉꼬부부가 아닙니다' 국제 부부로 변신한 수라 부부! 잉꼬 하우스 최초 국제 고민 부부 등장! 디시트렌드 07.14
갤러리 내부 검색
제목+내용게시물 정렬 옵션

오른쪽 컨텐츠 영역

실시간 베스트

1/8

뉴스

디시미디어

디시이슈

1/2