format long;
clear;
f0=193.1e12;
lambda0=c/f0;
gamma= 0.2; # 2W^-1/km;
P0=15.06921253847083;
alpha = 0;
D = -0.16e-12/1e-9/1e3;
?"beta2";
?beta2 = -lambda0^2*D/(2*pi*c);
width=0.0003430999009831653;
timeWindow = 4e-10;
Tfwhm=width*timeWindow;
?"T0g";
?T0g=Tfwhm/(2*sqrt(log(2)));
?"Ld";
?Ld=T0g^2/abs(beta2);
?"Lnl";
?Lnl=1/(gamma*P0);

?Leff = 0.2*Ld;

if (alpha ==0.0) {
L=Leff;}
else{
L=-log(1-Leff*alpha)/alpha;
}
?"L";
?L;

?Aeff=1;
?"n2";
?n2=gamma*c*Aeff/(2*pi*f0);
?"N";
?N=sqrt(gamma*P0*T0g^2/abs(beta2));
targetN = 10;
?"P0new";
?P0new=P0*targetN^2/N^2;
?"s";
?s=1/(2*pi*f0*T0g);
targets = 0.01;
?"T0gnew";
?T0gnew=1/(2*pi*f0*targets);
?"widthnew";
?widthnew=T0gnew*(2*sqrt(log(2)))/timeWindow;

?"-------------------------";