function [x,y]=MysteryPeak
x=[45:1/80:70]'; % time
a0= 1; 
a1= 53; 
a2=0.95; 
a3=1.3; 
EMG1=(1/2)*(a0/a3)*exp((a2^2/(2*a3^2))+(a1-x)/a3).*(erf((x-a1)/(sqrt(2)*a2)-a2/(sqrt(2)*a3))+(a3/abs(a3)));
b0= 1; 
b1= 55; 
b2=1.2; 
b3=1.3; 
EMG2=(1/2)*(b0/b3)*exp((b2^2/(2*b3^2))+(b1-x)/b3).*(erf((x-b1)/(sqrt(2)*b2)-b2/(sqrt(2)*b3))+(b3/abs(b3)));
c0= 1; 
c1= 57.5; 
c2=1.6; 
c3=1.3; 
EMG3=(1/2)*(c0/c3)*exp((c2^2/(2*c3^2))+(c1-x)/c3).*(erf((x-c1)/(sqrt(2)*c2)-c2/(sqrt(2)*c3))+(c3/abs(c3)));
Noise=.00;
y=EMG1+EMG2+EMG3+Noise.*randn(size(x));