diff options
| -rw-r--r-- | transverse/LaguerreGaussianE.m | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/transverse/LaguerreGaussianE.m b/transverse/LaguerreGaussianE.m index f95985d..adc004a 100644 --- a/transverse/LaguerreGaussianE.m +++ b/transverse/LaguerreGaussianE.m @@ -175,14 +175,14 @@ w=w_(q,lambda); if min(size(r))>=2
for u=1:size(params,1)
z(:,:,u) = a(u)...
- .* sqrt(2*factorial(p(u))/(1+(m(u)==0))/pi/(factorial( m(u)+p(u) )))/w(u)...
+ .* sqrt(2*factorial(p(u))/pi/(factorial( m(u)+p(u) )))/w(u)...
.* (sqrt(2)*r/w(u)).^m(u) .*exp(1i*signm(u)*m(u).*theta).* LaguerrePoly([p(u),m(u)],2*r.^2/w(u).^2)...
.* exp( -1i*2*pi/lambda(u)*r.^2/2/q(u));
end
else
for u=1:size(params,1)
z(:,u) = a(u)...
- .* sqrt(2*factorial(p(u))/(1+(m(u)==0))/pi/(factorial( m(u)+p(u) )))/w(u)...
+ .* sqrt(2*factorial(p(u))/pi/(factorial( m(u)+p(u) )))/w(u)...
.* (sqrt(2)*r/w(u)).^m(u) .* exp(i*signm(u)*m(u).*theta).* LaguerrePoly([p(u),m(u)],2*r.^2/w(u).^2)...
.* exp( -1i*2*pi/lambda(u)*r.^2/2/q(u));
end
|