%vary TQ

[lvect,gvect,pvect,peak(1),energy(1)]=qswitchsim(1,3,1000);
[lvect,gvect,pvect,peak(2),energy(2)]=qswitchsim(50,3,1000);
[lvect,gvect,pvect,peak(3),energy(3)]=qswitchsim(100,3,1000);
[lvect,gvect,pvect,peak(4),energy(4)]=qswitchsim(200,3,1000);
[lvect,gvect,pvect,peak(5),energy(5)]=qswitchsim(300,3,1000);
[lvect,gvect,pvect,peak(6),energy(6)]=qswitchsim(400,3,1000);
[lvect,gvect,pvect,peak(7),energy(7)]=qswitchsim(500,3,1000);
[lvect,gvect,pvect,peak(8),energy(8)]=qswitchsim(600,3,1000);
[lvect,gvect,pvect,peak(9),energy(9)]=qswitchsim(700,3,5000);
[lvect,gvect,pvect,peak(10),energy(10)]=qswitchsim(800,3,5000);
[lvect,gvect,pvect,peak(11),energy(11)]=qswitchsim(1000,3,5000);
[lvect,gvect,pvect,peak(12),energy(12)]=qswitchsim(2000,3,6000);
[lvect,gvect,pvect,peak(13),energy(13)]=qswitchsim(3000,3,6000);
TQs=[1 50 100 200 300 400 500 600 700 800 1000 2000 3000];
figure(2)
subplot(2,1,1)
plot(TQs, energy)
title('energy vs. T_Q')
ylabel('energy')
subplot(2,1,2)
plot(TQs, peak)
title('peak vs. T_Q')
xlabel('Normalized Time')
ylabel('peak')