clear;
clc;
 
sigma=input('sigma: ');
x=0:sigma:10
 
 y=exp(-0.2*x).*cos(x);
 y_poch=exp(-0.2*x).*cos(x)-exp(-0.2*x).*sin(x);
 y_poch_2=diff(y);
 
 plot(x,y,'m');
 ylabel('os y');
 xlabel('os x');
 hold on
 plot(x,y_poch,'y');