function [M]=plotcurve1(x0)
M=[]
for x=-10:0.001:10
M=[M testfunction1(x0)+(x-x0)*(16*x0.^3+4*x0+1)+(1/2)*(x-x0)*(48*(x0.^2)+4)*(x-x0)];
end


