Guest User

Untitled

a guest
Nov 13th, 2018
90
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.44 KB | None | 0 0
  1. Clear[r, re, p, pmax, delta, imagesize, delta]
  2. ClearSystemCache[]
  3. re[0, r_] := Sqrt[8/Pi]*((1 - r)/r)^(1/4)*1;
  4. re[1, r_] := Sqrt[8/Pi]*((1 - r)/r)^(1/4)*-1*2*(1 - 2*r);
  5. re[p_, r_] := re[p, r] = Sqrt[8/Pi]*((1 - r)/r)^(1/4)*(-1)^p*(re[1, r]*re[p - 1, r] - re[p - 2, r]);
  6. imagesize = 32;
  7. pmax = 10;
  8. delta = 2/imagesize;
  9. Table[r = Sqrt[x^2 + y^2]; re[pmax, r], {x, -1 + delta/2, 1 - delta/2, delta}, {y, 1 - delta/2, -1 + delta/2, -delta}];
Add Comment
Please, Sign In to add comment