x=[-10:0.01:10];
n=size(x,2);
y=zeros(n,1);
for i=1:n
  y(i)=exp(-x(i)^2/2)*(1-x(i)^2);
end;
plot(x,y,'black');
axis([-10 10 -1 1.5]);
axis on;
