11.01.2015 Views

Geant4 Simulations for the Radon Electric Dipole Moment Search at

Geant4 Simulations for the Radon Electric Dipole Moment Search at

Geant4 Simulations for the Radon Electric Dipole Moment Search at

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

1 % LegendreP.m<br />

2 function P = LegendreP(n,x)<br />

3 if n == 0<br />

4 P = 1;<br />

5 elseif n == 2<br />

6 P = (1/2)*(3*(x)ˆ2−1);<br />

7 elseif n == 4<br />

8 P = (1/8)*(35*(x)ˆ4−30*(x)ˆ2+3);<br />

9 elseif n == 6<br />

10 P = (1/16)*(231*(x)ˆ6−315*(x)ˆ4+105*(x)ˆ2−5);<br />

11 elseif n == 8<br />

12 P = (1/128)*(6435*(x)ˆ8−12012*(x)ˆ6+6930*(x)ˆ4−1260*(x)ˆ2+35);<br />

13 elseif n == 10<br />

14 P = (1/256)*(46189*(x)ˆ10−109395*(x)ˆ8+90090*...<br />

15 (x)ˆ6−30030*(x)ˆ4−3465*(x)ˆ2−63);<br />

16 else<br />

17 error('Legendre polynomial not found.');<br />

18 return;<br />

19 end<br />

1 function out = Plot3DSurface(d<strong>at</strong>a)<br />

2 resTheta = 5;<br />

3 resPhi = 5;<br />

4 out = zeros(181,1);<br />

5 outRes = zeros((180)/resTheta,1);<br />

6 resSurface = zeros((180/resTheta),(360/resPhi));<br />

7 sum = 0;<br />

8 <strong>for</strong> i=0:1:180<br />

93

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!