t=0:.01:40;
f=-(t.*cos(t)-sin(t))/2;
plot(t,f)



% a=4;
% b=20;
% t=0:.01:40;
% f=zeros(1,length(t));
% f(a*100+1:end)=sin(t(a*100+1:end)-a);
% f(b*100+1:end)=f(b*100+1:end)-sin(t(b*100+1:end)-b);
% plot(t,f)
% axis([0 40 -2.4 2.4])

% 
% t=0:.01:10;
% f=zeros(length(t),1);
% f(101:end)=exp(-(t(101:end)-1)).*(cos(t(101:end)-1)-sin(t(101:end)-1));
% plot(t,f)
% axis([0 10 -.5 1.3])