clear
x=-3:.1:3;
figure(1)
plot(x,f1(x))
grid
xlabel('x')
ylabel('y')
title('y=f1(x)=x^3-4*x')