総和が求められません。

Question

コード
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<time.h>

#define T 1000
#define N 10

int main(int argc, char **argv){
    int i, j, k, l, m;
    double x[N], xb[N];
    double dt;
    double D;
    int init_seed;
    double w1, w2;
    double u1, u2;
double nijo[N];
double wa = 0.0;
double var;

    init_seed = time(NULL);

    srand(init_seed);

    D = 1.0;
    dt = 1.0/(double)T;


    printf("%f ", 0.0);
    for(i=0;i<N;i++){
        xb[i] = 0.0;

        printf("%f ", xb[i]);
    }

    printf("\n");



    for(j=1;j<=T;j++){
        printf("%f ", (double)j*dt);

        for(i=0;i<N;i++){
            u1 =  ( double ) rand ()/( RAND_MAX +1.0);
            u2 =  ( double ) rand ()/( RAND_MAX +1.0);
            w1 = sqrt(-2.0*log(u1))*cos(2.0*M_PI*u2);
            w2 = sqrt(-2.0*log(u1))*sin(2.0*M_PI*u2);

            x[i] = xb[i] + sqrt(2.0*D)*w1*pow(dt, 0.5);



      }

        for(i=0;i<N;i++){
            xb[i] = x[i];
        }

        for(k=0;k<N;k++){
          nijo[k] = xb[k]*xb[k];
        }
        for(l=0;l<N;l++){
          wa = wa + nijo[l];
        }
        for(m=0;m<N;m++){
          var = wa/N ;
          printf("%lf", var);
        }


        printf("\n");
    }

	return 0;
}

おそらく二乗の計算まではできたのですが、そこから先ができません。どなたか教えていただけませんか

Comments

[BeatStar]

そもそも何がしたいの？
少なくとも単に「総和」だけじゃ伝わりません。

今の状態は『足し算をしたい』と言いながら『微積分をやろうとしている』感じです。

まず、そもそもなにがしたいのか、『どうなれば正解か』ぐらいは書きましょう。

回答者は質問者ではないので、提示されないとわかりませんよ。

[tackle.8]

一次元ブラウン運動の拡散現象について数値的に求めていました。
０＜時間ｔ＜１の間に０からの距離ｘに存在する粒子（N個）を表したのがxbです。
要するに各時間におけるxbの二乗の和を総数Nで割って分散を出したいのですが、二乗和がうまくできないといった状況です。
説明不足ですみませんでした。

[BeatStar]

説明不足なのはまだいいですが、質問は編集できるので編集してください。
回答者は質問本文を見て判断します。

[tackle.8]

わかりました。
それでは、この問題についてどうお考えですか

[BeatStar]

あと、『うまくいかない』というのはやめましょう。そういうのは『主観』でしかありません。

もしかすると、質問者さんの(暗算時の)計算ミスかもしれないし、環境によるものかもしれないし、コンパイルエラーのことかもしれないし、…と色んな可能性がありますので。

[Zuishin]

あなたが話しているのは同じ授業を受けている友だちではなく、全く接点のない相手です。違う言葉を使っていたのでは話にならないので、文脈や用語のすり合わせがまず必要になります。

質問はこうです。「二乗はできたのですが総和の求め方がわかりません。どうすればいいですか？」

「総和」というのは、出た数を全て足すことですが、そこまでは同じ意味で使っていますか？
もしそうならば、なぜそれができないのか説明してください。
第三者が見る上で、小学校一年生で習うはずの足し算ができないという質問の意味はわかりません。

Answer 1

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<time.h>
#define T 1000
#define N 10

int main(int argc, char **argv)
{
  int i, j, m;
  double xb[N] = {0};
  double dt = 1. / T, D = 1;
  double w1, w2, u1, u2;
  double wa = 0.0;
  double var;
  srand(time(NULL));
  for(m = 0; m <= N; m++){
    printf("%lf", 0);
  }
  for(j = 1; j <= T; j++){
    /* wa = 0; */
    printf("%f ", (double)j*dt);
    for(i = 0; i < N; i++){
      u1 = (double) rand () / (RAND_MAX +1.0);
      u2 = (double) rand () / (RAND_MAX +1.0);
      w1 = sqrt(-2.0*log(u1))*cos(2.0*M_PI*u2);
      w2 = sqrt(-2.0*log(u1))*sin(2.0*M_PI*u2);
      xb[i] += sqrt(2.0*D)*w1*pow(dt, 0.5);
      wa += xb[i]*xb[i];
    }
    for(m = 0; m < N; m++){
      var = wa/N;
      printf("%lf", var);
      /* printf("%lf", xb[m]); */
    }
    printf("\n");
  }
  return 0;
}

チキンマルナゲットな貴公にはヒントをやろう、貴公の書いたコードを我が書くとこうなる。何が足りていないのか考え給え。コメントもよく読むのだぞ。

Comments

[tackle.8]

分散varはfor文にいれなくてよかったです
他は分かりません。

[majiponi]

各ループごとにwaを0にしないで前の値を引き継ぐと何が起こるか、考えて見ましょう。

Answer 2

for(l=0;l<N;l++){

      wa = wa + nijo[l];
    }
    for(m=0;m<N;m++){
      var = wa/N ;
      printf("%lf", var);
    }

最初の for で確かに総和出してますが、
その次の for ではなにをしてるんでしょうか。

Comments

[tackle.8]

総和を総数Nでわっています。

[y_waiwai]

forループで繰り返しNで割ってるようですが、そこらへんはどんなもんでしょうか。

そもそも実行結果はどうなってるのか提示しようよ