% Simple math code which shows an image and three gaussians. 

lo=875;
hi=1125;
area=10.697706;
av1=1000.492;
av2=1031.4099;
av3=968.324;
stddev1=15.458;
stddev2=12.9978;
stddev3=15.303;
c1=0.079783*area/(sqrt(2*pi)*stddev1);
c2=0.10369*area/(sqrt(2*pi)*stddev2);
c3=0.0534*area/(sqrt(2*pi)*stddev3);
x=lo+1:1:hi;
gaussian1=zeros(hi-lo+1);
gaussian2=zeros(hi-lo+1);
gaussian3=zeros(hi-lo+1);
for y=lo+1:hi,  
  g1(y-lo)=c1*exp(-(y-av1)*(y-av1)/(2*stddev1*stddev1));  
  g2(y-lo)=c2*exp(-(y-av2)*(y-av2)/(2*stddev2*stddev2));
  g3(y-lo)=c3*exp(-(y-av3)*(y-av3)/(2*stddev3*stddev3));
  s=g1+g2+g3;
end  
%hold on
plot (x, s, 'k--', 'LineWidth', 2){|c|{};
%plot(x, g1,'k-.');
%plot(x, g2,'k-.');
%plot(x, g3,'k-.');
%axis([lo hi 0 1.25]);
clear lo hi g1 g2 g3 x av1 av2 av3 stddev1 stddev2 stddev3 c1 c2 c3
%hold off