x=[0.0 0.2 0.8 1.1 1.3]
f=inline('exp(3*x).*cos(x)','x');
y=f(x)
p=polyfit(x,y,4);
xp=0.9;
polxp=polyval(p,xp);
yp=f(xp);
error1=abs(polxp-yp);
xx=0:0.13:1.5;
xx=[xx 1.5];
yy=f(xx);
sxp=spline(xx,yy,xp);
error2=abs(sxp-yp);
%Imprimimos la solucion:
fprintf('SOLUCIO:\n');
fprintf('(a) VALOR DEL POLINOMIO INTERPOLADOR EN EL PUNTO XP = %f:\n',xp)
fprintf('\t\tP(%f) = %23.15e\n',xp,polxp)
fprintf('    ERROR EN LA INTERPOLACION POLINOMIAL :\n')
fprintf('\t\t|P(%f)-f(%f)| = %23.15e\n',xp,xp,error1)
fprintf('(b) ERROR EN LA INTERPOLACION POR SPLINES:\n');
fprintf('\t\t|S(%f)-f(%f)| = %23.15e\n',xp,xp,error2);