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

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.10.25 Thursday 23:48
    0
      category:- | by:スポンサードリンク | - | -
      Comment








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

      PR
      Calender
          123
      45678910
      11121314151617
      18192021222324
      252627282930 
      << November 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