function [IV]=funcpl(P,x,y)
% Intensitätsverteilung einer Punkt-Linienkombination
% Durchm. 5 , Länge 15, Breite 5

ix=30; % Wertebreich in x
iy=15;  % Wertebereich in y
dx=0.1; % Pixelgröße in x in mm
dy=0.1; % Pixelgröße in y in mm

n=10000; % Supergaussfaktor

w0=2.5; % Radius Kreis
wx=7.5; % 1/2 Breite x in mm
wy=w0; % 1/2 Breite y in mm


vxp=wx+w0; % Verschiebung Punkt in x
vyp=0; % Verschiebung Punkt in y
vxl=0; % Verschiebung Linie in x
vyl=0; % Verschiebung Linie in y

[ix,iy]=meshgrid(-1*ix:dx:ix,-1*iy:dy:iy);

% Intensitätsverteilung normiert auf 1
% Verhältnis Punkt/Linie=2/1
punkt=2*exp(-2*((((ix-vxp-x).^2+(iy-vyp-y).^2).^0.5)./w0).^(2*n)); 
linie=exp(-2*(((ix-vxl-x)./wx).^(2*n)+(((iy-vyl-y)./wy).^(2*n))));

IV1=punkt+linie; 
Pges1=sum(sum(IV1)); % gesamte Leistung in normierter Intensitätsverteilung

% Leistungsinhalt anpassen
IV2=(P/Pges1).*IV1;
Pkontr=sum(sum(IV2)); % Kontrollsumme
Pkontr=num2str(Pkontr)

% Intensität in W/cm^2 berechnen
IV=IV2/((0.1*dx)*(0.1*dy));