function y = fd(x)
y = (3*x-4)*x+1;