function y=func2(x)
y=1+x-x^2/4;