x=0:pi/40:2*pi;
a=2;b=5;
y=b*exp(-a*x).*sin(b*x).*(0.012*x.^4 -0.15*x.^3 +0.075*x.^2 +2.5.*x);
z = y.^2;
w(:,1)=x;
w(:,2)=y;
w(:,3)=z;
format short g;
w
plot(x,y,'r--','LineWidth',1.5,...
'Marker','p','MarkerSize',14,...
'MarkerEdgeColor','w',...
'MarkerFaceColor','r');hold on;
plot(x,z,'b-s','MarkerFaceColor','g');



