1
$\begingroup$

I have following code in my Mathematica:

P[k_?NumericQ, L_?NumericQ] := P[k, L] = NIntegrate[Exp[I*k*r]*1/r^(2/3), {r, 0, L}];

i = 1;

rhot = {};

a[k_] := RandomVariate[NormalDistribution[0, Sqrt[Re[P[k, 100]]]]] + 
   I*RandomVariate[NormalDistribution[0, Sqrt[Im[P[k, 100]]]]];

While[i < 3, 
      alist = Table[a[k], {k, 2*Pi/100, 2*Pi, 2*Pi/100}]; 
      rslist = Fourier[alist]; 
      rhotl = Total[rslist];
      Append[rhot, rhot1]; 
      i++]

Unfortunately, the while loop does not give any output. Any idea why this is so?

| improve this question | |
$\endgroup$
  • 4
    $\begingroup$ Try with Table or AppendTo. $\endgroup$ – Kuba♦ Aug 25 '16 at 9:01
  • $\begingroup$ You mean change the statement Append[rhot, rhot1]; to AppendTo[rhot, rhot1]; ? $\endgroup$ – titanium Aug 25 '16 at 9:05
  • 5
    $\begingroup$ rhot = Append[rhot, Total[rslist]]; Append[] doesn't actually dynamically change a list, you still need =. With AppendTo you don't need = $\endgroup$ – Feyre Aug 25 '16 at 9:19