function P1()
% Checks out Quintic
randn('seed',0)
d = randn(4,1);
y = randn(4,1);
clc
format long
a = Quintic(d,y)
err = abs(y - (((((a(6)*d + a(5)).*d + a(4)).*d + a(3)).*d+a(2)).*d+a(1)));
dp1 = (((5*a(6)*d(1) + 4*a(5))*d(1) + 3*a(4))*d(1) + 2*a(3))*d(1) + a(2);
dp4 = (((5*a(6)*d(4) + 4*a(5))*d(4) + 3*a(4))*d(4) + 2*a(3))*d(4) + a(2);
ddp1 = ((20*a(6)*d(1) + 12*a(5))*d(1) + 6*a(4))*d(1) + 2*a(3);
ddp4 = ((20*a(6)*d(4) + 12*a(5))*d(4) + 6*a(4))*d(4) + 2*a(3);
fprintf('| p(d(1)) - y(1) | = %15.10e\n',err(1))
fprintf('| p(d(2)) - y(2) | = %15.10e\n',err(2))
fprintf('| p(d(3)) - y(3) | = %15.10e\n',err(3))
fprintf('| p(d(4)) - y(4) | = %15.10e\n\n',err(4))
fprintf('p''(d(1)) = %18.15f\n',dp1)
fprintf('p''(d(4)) = %18.15f\n\n',dp4)
fprintf('p"(d(1)) = %18.15f\n',ddp1)
fprintf('p"(d(4)) = %18.15f\n',ddp4)