// Mais algumas bricadeiras com o scilab
// Aproximacao do logaritmo:
// ln(1+x)= x -x^2/2 + x^3/3 -x^4/4 + x/5 ...
deff('y=ln(x,n)', 's=x-1; L=0; for i=1:n, L=L+(-1)**(i-1)*s**(i)/i; end; y=L;');
// usando a formula ln[(1+x)/(1-x)]=2(x+ x^3/3+x^5/5]
deff('y=ap_ln(x,n)', 's=(x-1)/(x+1); L=0; for k=1:2:n, L=L+s**k/k; end; y=2*L;');