clc;?clear?all;?close?all;?clf;
cputime=0;
tic;
ln=1;
i=sqrt（-1）;
Po=.00064;?%input?pwr?in?watts
alpha=0;?%?Fiber?loss?value?in?dB/km
alph=alpha/（4.343）;?%Ref?page#55?eqn?2.5.3?Fiber?optic?Comm?by?GP?Agrawal
gamma=0.003;?%fiber?non?linearity?in?/W/m
to=125e-12;?%initial?pulse?width?in?second
C=-2;?%Input?chirp?parameter?for?first?calculation
b2=-20e-27;?%2nd?order?disp.?（s2/m）
Ld=（to^2）/（abs（b2））;?%dispersion?length?in?meter
pi=3.1415926535;
Ao=sqrt（Po）;?%Amplitude
%----------------------------------------------------------
tau?=-?4096e-12:1e-12:?4095e-12;%??dt=t/to
?dt=1e-12;
?rel_error=1e-5;
?h=1000;%?step?size
for?ii=0.1:0.1:1.5?%the?various?fiber?lengths?can?be?varied?and?this?vector?can?be?changed
z=ii*Ld;
????u=Ao*exp（-（（1+i*（-C））/2）*（tau/to）.^2）;%page#47?G.P.AGrawal
????figure（1）
???plot（abs（u）‘r‘）;
???title（‘Input?Pulse‘）;?xlabel（‘Time‘）;?ylabel（‘Amplitude‘）;
grid?on;
hold?on;
l=max（size（u））;??
%%%%%%%%%%%%%%%%%%%%%%%
fwhm1=find（abs（u）>abs（max（u）/2））;
fwhm1=length（fwhm1）;
dw=1/l/dt*2*pi;
w