clc for n=1:30000; t(n)=n/1000000; if t(n)<1.5/1000; y(n)=0; elseif t(n)<15/1000; y(n)=3.16*sin((120*pi*t(n))-(39/180)*pi)+1.69*exp(-342*t(n))-0.707*exp(-58.5*t(n)); else y(n)=4.27*sin((120*pi*(t(n)-0.015))-(94.3/180)*pi)+exp(-100*(t(n)-0.015))*0.918*cos(100*(t(n)-0.015))-1.37*sin(100*(t(n)-0.015)); %y(n)=4.27*sin((120*pi*t(n))-(58.3/180)*pi)+1.92*exp(-100*(t(n)))*cos(100*(t(n)+(47.6/180)*pi)); end end plot(t,y)