for k = 1:101 x(k) = (k-1)*0.01; u(k) = sin(6*pi*x(k))*exp(-x(k)); v(k) = cos(6*pi*x(k))*exp(-x(k)); z(k) = 0; end plot(x,u,'k-',x,v,'r--',x,z,'k:')