スポンサーサイト

2018.06.30 Saturday
0

    一定期間更新がないため広告を表示しています

    category:- | by:スポンサードリンク | - | -

    実行と修正とこころ強さと

    2018.06.20 Wednesday 23:48
    0

      【劇団NONNY102回公演】

      『プラヌラ』

      脚本: 高石紗和子 演出: 三木駿

      日時:6/30() 14:00

         7/1() 13:00

      場所:studio CoCoLo(千葉駅徒歩5)

       

      しばらくのお付き合いを願います。本日のブログ担当はあるふぃこと、2年目の盡澤です。

       

      本番が近づいてきましたね、、、、今日は2回目の通しでした。舞台美術や小道具がだんだんと揃ってきて、演技するのがますます楽しくなってます!
      ただ、まだまだ思うようにいかないところもあり・・・。本番まで気は抜けませんね!

       

      さて、突然ですが、ここで僕の力作を見ていただこう。

       

      #include <stdio.h>

      #include <math.h>

      #define DIM 4

      void bsort(double array_A[][DIM],double copy[][DIM],double vec_b[], int k)

      {

       int i, bigest = k, j;

       double max = fabs(array_A[k][k]);

       double box;

       for(i = k ; i < DIM; i++){

                       if(max < fabs(array_A[i][k])){

                                        max = fabs(array_A[i][k]);

                                        bigest = i;

                       }

       }

       for(j = 0; j < DIM; j++){

                       box = array_A[bigest][j];

                       array_A[bigest][j] = array_A[k][j];

                       array_A[k][j] = box;

       }

       box = vec_b[bigest];

       vec_b[bigest] = vec_b[k];

       vec_b[k] = box;

       for (i = 0; i < DIM; i++){

                       for (j = 0; j < DIM; j++){

                                        copy[i][j] = array_A[i][j];

                       }

       }

      }

      void new_array(double array_A[][DIM], double copy[][DIM],double vec_b[], int i, int j, int k)

      {

       double lamd;

       if(i > k && j >= k){

                       lamd = -copy[i][k] / copy[k][k];

       } else {

                       lamd = 0;

       }

       array_A[i][j] = array_A[i][j] + lamd * array_A[k][j];

       if (j == DIM - 1) {

                      vec_b[i] = vec_b[i] + lamd * vec_b[k];

       } else {

                       vec_b[i] = vec_b[i];

       }

      }

      double sigma(double array_A[][DIM], double vec_x[], int k)

      {

       int i;

       double sum;

       for(i = k + 1; i < DIM; i++){

                       sum += array_A[k][i] * vec_x[i];

       }

       return sum;

      }

      int main(void)

      {

       double array_A[DIM][DIM] = {{2,4,1,-3},{-1,-2,2,4},{4,2,-3,5},{5,-4,-3,1}};

       double copy[DIM][DIM] = {{2,4,1,-3},{-1,-2,2,4},{4,2,-3,5},{5,-4,-3,1}};

       double vec_b[DIM] = {0,10,2,6};

       double vec_x[DIM];

       double sum[DIM];

       int i, j, k, l;

       for(k = 0; k < DIM - 1; k++){

       

                       if(array_A[k][k] == 0){

                                        l = k;

                                        break;

                        }

       

                       for(i = 0; i < DIM; i++){

                                        for(j = 0; j < DIM; j++){

       

                                                        new_array(array_A, copy, vec_b, i, j, k);

       

                                        }

                      }

                       for (i = 0; i < DIM; i++){

                                        for (j = 0; j < DIM; j++){

                                                        copy[i][j] = array_A[i][j];

                                        }              

                       }

       }

       if (array_A[k][k] == 0){

                        for(k = l; k < DIM ; k++){

                      

                                        bsort(array_A,copy, vec_b, k);

       

                                        for(i = 0; i < DIM; i++){

                                                        for(j = 0; j < DIM; j++){

                                                                         new_array(array_A, copy, vec_b, i, j, k);

                                                        }

                                       }

                        }

       }

       vec_x[DIM - 1] = vec_b[DIM - 1] / array_A[DIM - 1][DIM - 1];

       

       for (i =DIM - 2; i >= 0; i -= 1){

                       sum[i] = sigma(array_A, vec_x, i);

                       printf("¥n");

                      vec_x[i] = (vec_b[i] - sum[i]) / array_A[i][i];

       }

       printf("x = %f¥ny = %f¥nz = %f¥nw = %f¥n",vec_x[0], vec_x[1], vec_x[2], vec_x[3]);

       return 0;

      }

       

      これはC言語で書かれたプログラムで、実行すると、なななんと!!!

      2x+4y+z+w=0

      -x-2y+2z+4w=10

      4x+2y-3z+5w=2

      5x-4y-3z+w=6

      の連立一次方程式の解を得ることができるのだ!! !!(゜□゜;ノ)ノ エエエッ

       

      ・・・・・・・・・・・・・意味わかんねぇwwww

       

       

      結局なにが言いたいかというと・・・季節の変わり目で、体調崩しやすいから、何事も無理しないでね!!ってことである。

      category:- | by:劇団NONNYcomments(0)trackbacks(0)

      スポンサーサイト

      2018.06.30 Saturday 23:48
      0
        category:- | by:スポンサードリンク | - | -
        Comment








           
        Trackback
        この記事のトラックバックURL

        PR
        Calender
              1
        2345678
        9101112131415
        16171819202122
        23242526272829
        30      
        << September 2018 >>
         
        profilephoto
        公演情報
        第94回公演
        第94回公演
        公式サイト
        Selected entry
        Category
        Archives
        Recent comment
        • 「雨の街」(終)
          古庄
        • お邪魔します、はじめまして
          two
        • お邪魔します、はじめまして
          o-2ka
        • お邪魔します、はじめまして
          FCarter
        • 新歓公演に向けて
          tokky
        • 観劇感激
          two
        • 演出その2より
          セキネ
        • 演出その2より
          その1
        • 七夕公演を終えて 〜今何時だと思ってんだこみやまっw〜
          セキネ
        • 七夕公演を終えて 〜今何時だと思ってんだこみやまっw〜
          こめ
        Recent trackback
        Recommend
        Link
        Profile
        Search
        Others
        Mobile
        qrcode
        Powered
        無料ブログ作成サービス JUGEM