home *** CD-ROM | disk | FTP | other *** search
- "File APPROX.MTH, copyright (c) 1990-2000 by Texas Instruments Incorporated"
-
- PADE_CF(a,m,n):=VECTOR(VECTOR(a SUB (m-n+i+j),j,1,n+1),i,1,n)
-
- PADE_DV(a,m,n):=ELEMENT(ROW_REDUCE(PADE_CF(a,m,n))`,n+1)
-
- PADE_NV(a,u,m,n):=VECTOR(a SUB (k+1)-SUM(u SUB (n+1-i)*a SUB (k+1-i),i,1,k),k,~
- 1,MIN(m,n))
-
- PADE_AUX3(a,v,u,x,c,m,n):=(a SUB 1+SUM(v SUB k*(x-c)^k,k,1,m))/(1-SUM(u SUB (n~
- +1-k)*(x-c)^k,k,1,n))
-
- PADE_AUX2(a,u,x,c,m,n):=PADE_AUX3(a,PADE_NV(a,u,m,n),u,x,c,m,n)
-
- PADE_AUX1(a,x,c,m,n):=PADE_AUX2(a,PADE_DV(a,m,n),x,c,m,n)
-
- PADE(y,x,x0,n,d):=PADE_AUX1(VECTOR(LIM(DIF(y,x,m_),x,x0)/m_!,m_,0,n+d+1),x,x0,~
- n,d)
-
- RECURRENCE(u,v,v0,m):=REVERSE(ITERATE(APPEND([u],v),v,v0,m-DIM(v0)))
-
- CHEBYSHEV_T_LIST(n,x):=RECURRENCE('(2*x*v_ SUB 1-v_ SUB 2),v_,[x,1],n)
-
- CHEBYSHEV_T_SERIES_AUX(i,n):=PROG(IF(i=0,RETURN(VECTOR(COMB(r,r/2)/2^r,r,0,n,2~
- ))),VECTOR(COMB(r,(r-i)/2)/2^(r-1),r,i,n,2))
-
- CHEBYSHEV_T_SERIES(f,x,a,b,k,n,T_list,T_Coefs,Taylor_Coeff_Even,Taylor_Coeff_O~
- dd,s_,t_,size,list,size_T_e,size_T_o):=PROG(IF(k>n,RETURN("k must be <= n")),f~
- :=SUBST(f,x,a+(b-a)/2*(x+1)),T_list:=CHEBYSHEV_T_LIST(k+1,x),T_Coefs:=VECTOR(C~
- HEBYSHEV_T_SERIES_AUX(s,n),s,0,n),Taylor_Coeff_Even:=VECTOR(SUBST(DIF(f,x,r_),~
- x,0)/r_!,r_,0,n,2),Taylor_Coeff_Odd:=VECTOR(SUBST(DIF(f,x,r_),x,0)/r_!,r_,1,n,~
- 2),s_:=1,t_:=1,c:=0,list:=[],size_T_e:=DIM(Taylor_Coeff_Even),size_T_o:=DIM(Ta~
- ylor_Coeff_Odd),LOOP(IF(c>k,exit),list:=APPEND(list,[T_Coefs SUB t_ . Taylor_C~
- oeff_Even SUB [s_,...,size_T_e]]),c:+1,IF(c>k,exit),list:=APPEND(list,[T_Coefs~
- SUB (t_+1) . Taylor_Coeff_Odd SUB [s_,...,size_T_o]]),t_:+2,s_:+1,c:+1),EXPAN~
- D(SUBST(list . T_list,x,-1+2*(x-a)/(b-a)),Rational))