% Taylor_sin.m % Example of Taylor polynomial : sin(x) % function y = Taylor_sin(n,x,a) function p = Taylor1(n,x,a)
p = sin(a)*ones(size(x)); for k=1:n Dkf = Dk_sin(a,k); p = p + Dkf*(x-a).^k/factorial(k); end subplot(1,2,1); plot(x,sin(x),x,p,'m',a,sin(a),'ro') tie = ['Taylor polynomial of degree ', num2str(n),' for f(x)=sin(x).']; title(tie); grid subplot(1,2,2); plot(x,sin(x)-p,a,sin(a),'ro') title('Error plot of sin(x) and its Taylor polynomial')
function y = Dk_sin(a,k) if rem(k,4)==0 y = sin(a); elseif rem(k,4)==1 y = cos(a); elseif rem(k,4)==2 y = -sin(a); else y = -cos(a); end
|
|