%vary Lf
clear
figure(1)
[lvect,gvect,pvect,peak(1),energy(1)]=qswitchsim(100,3, .01,250);
[lvect,gvect,pvect,peak(2),energy(2)]=qswitchsim(100,3, .1,250);
[lvect,gvect,pvect,peak(3),energy(3)]=qswitchsim(100,3, .3,250);
[lvect,gvect,pvect,peak(4),energy(4)]=qswitchsim(100,3, .4,250);
[lvect,gvect,pvect,peak(5),energy(5)]=qswitchsim(100,3, .45,250);
[lvect,gvect,pvect,peak(6),energy(6)]=qswitchsim(100,3, .47,250);
Lf=[.01 .2 .3 .4 .45 .47];
figure(2)
subplot(2,1,1)
plot(Lf, energy)
title('energy vs. L_f')
ylabel('energy')
subplot(2,1,2)
plot(Lf, peak)
title('peak vs. L_f')
xlabel('Normalized Time')
ylabel('peak')