%PLOT_RAMAN3
%
a=zeros(50,1); y1=zeros(50,1); y2=zeros(50,1); y3=zeros(50,1);
y4=zeros(50,1); y5=zeros(50,1); y6=zeros(50,1); y7=zeros(50,1);
t1=1/100; t2=1/10; t3=1/2; t4=1; t5=2; t6=10; t7=100;
for i=1:50
  a(i)=.02*i;
  y1(i)=feval(symengine,'igamma',a(i),t1)/gamma(a(i));
  y2(i)=feval(symengine,'igamma',a(i),t2)/gamma(a(i));
  y3(i)=feval(symengine,'igamma',a(i),t3)/gamma(a(i));
  y4(i)=feval(symengine,'igamma',a(i),t4)/gamma(a(i));
  y5(i)=feval(symengine,'igamma',a(i),t5)/gamma(a(i));
  y6(i)=feval(symengine,'igamma',a(i),t6)/gamma(a(i)); 
  y7(i)=feval(symengine,'igamma',a(i),t7)/gamma(a(i));
end
plot(a,y6,'r','LineWidth',2)
set(gca,'FontSize',15)
axis([0 1.3 0 5e-5])
xticks([0 .2 .4 .6 .8 1])
yticks([0 1e-5 2e-5 3e-5 4e-5 5e-5])
text(1.05,4.5e-5,'t=10','FontSize',15)
