powersin.m

有趣的可视的数值方法 出自网站http://www.mathworks.com/moler

function s = powersin(x)
%POWERSIN  Power series for sin(x).
%  y = POWERSIN(x) tries to compute sin(x) from its power series.
s = 0;
t = x;
n = 1;
while s+t ~= s;
   s = s + t;
   t = -x.^2/((n+1)*(n+2)).*t;
   n = n + 2;
end