www.pudn.com > NumericalComputingwithMatlabCode.zip > powersin.m, change:2003-07-27,size:238b

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