soliton_fission_supercontinuum_generation.m

lambda = [1.2	1.25	1.3	1.35	1.4	1.45	1.5	1.55	1.6	1.65	1.7	1.75	1.8	1.85	1.9	1.95	2];
lambda = lambda .* 1e-6;
neff = [3.35903905970416 	3.33765501184926 	3.31681134944427 	3.29635276597645 	3.27615314956330 	3.25610903639629 	3.23613470840525 	3.21615847619334 	3.19611982800696 	3.17596722007002 	3.15565634698692 	3.13514877575541 	3.11441085783588 	3.09341285564337 	3.07212823607959 	3.05053309538915 	3.02860568668973];
figure(1);
plot(lambda,neff,'r-'),hold on,grid on;
axis([1.2e-6,2e-6,2.8,3.5]);
dndlambda = [diff(neff)./diff(lambda),0];
ddnddlambda = [diff(dndlambda)./diff(lambda),0];
beta2 =1e24 .* lambda .^3 ./ (2*pi*C^2) .* ddnddlambda ;
figure(2);
plot(lambda,beta2,'b-'),hold on,grid on;
axis([1.2e-6,1.9e-6,-3,3]);