arcsin.m

来自「信号与系统分析以matlab实现」· M 代码 · 共 21 行

function y=arcsin(x)
n=1;
p=x;
for n=1:100
    j=1;
    l=1;
    m=2*n;
    for i=1:n
        j=j*i;
    end
    for i=1:m
        l=l*i;
    end
    p=p+(l*(x^(m+1)))/((2^m)*(j^2)*(m+1));
    p1=(l*(x^(m+1)))/((2^m)*(j^2)*(m+1));
    if ( p1<=0.002)
       break;
  end
end
y=p