function [f g]=intpolyfunc(x);

f=1/6*x.^6 + x.^5 - 3/4*x.^4 + 2/3*x.^3 - 1/2*x.^2 - 20*x;
g=x.^5 + 5*x.^4 -3*x.^3 + 2*x.^2 - x - 20;