function Intensverteilung

%---Konst. Parameter

lambda = 1e-6;
z0     = 10e-6; %(Rayleigh-Länge)
z      = linspace(-5,5,1000);

%---Berechnung

w0 = lambda*z0/pi;
wz = w0^2*(1+(z./z0).^2);
h  = (w0./wz).^2;

figure(1)
subplot(2,1,1)
plot (z,h)
grid on
subplot(2,1,2)
plot (z,wz)
grid on

%---Berechnung Intensverteilung

I = 1;
w0 = 1;
%wz = rand(5);
x = rand(5);
y = rand(5);

In = I * (w0./wz).^2 * (exp(-2*(x.^2+y.^2)./wz.^2));

%I=1;
%[x,y] = meshgrid(0:0.5:10, 0:0.5:10);
%In = I*(w0./wz).^2 * (exp(-2*(x.^2+y.^2)./wz.^2));
%surf(x,y,In)