Advanced Search

Question: system of equation

Question:system of equation

kagu 60
complex student
December 09 2009
0

this is my program to calcul pade approximant at n points complex or real. and i want to excute this program with L =40 or more but it can't.

 

L:=4:
P:=add(p[k]*x^k,k=0..L):
Q:=1+add(q[k]*x^k,k=1..L+1):
f:=(-1/x)*ln(1-x):

a:=seq(coeftayl(f*Q-P,x=I,k),k=0..L):
a1:=seq(coeftayl(f*Q-P,x=-I,k),k=0..L):


s:=fsolve({a, a1});
AP:=evalf(eval(P/Q,s)):
 please can you help to improve this program.
Please Wait...