algebraic procedure taylorp(f,c,n,ax,bx,ay,by);
% Returns and plots the n-th Taylor polynomial centered at c
% of a function f in the interval [ax, bx].
begin scalar p;
p:=taylor(f,x,c,n);
p:=taylortostandard p;
plotreset;
plot(f,p,x=(ax .. bx),y=(ay .. by));
return p;
end;
end;