home *** CD-ROM | disk | FTP | other *** search
/ Chip 2002 June / Chip_2002-06_cd1.bin / zkuste / derive / download / Setup.exe / %MAINDIR% / Math / Numeric.mth < prev    next >
Encoding:
Text File  |  2002-05-15  |  1.5 KB  |  37 lines

  1. "File NUMERIC.MTH, copyright (c) 1990-2000 by Texas Instruments Incorporated"
  2.  
  3. DIF_NUMERIC(y,x,x0,h,n):=IF(NUMBER(n),SUM((-1)^k*COMB(n,k)*LIM(y,x,x0-n*h/2+k*~
  4. h),k,0,n)/(-h)^n,DIF_NUMERIC(y,x,x0,h,1),DIF_NUMERIC(y,x,x0,h,1))
  5.  
  6. DIF_NUM(y,x,x0,h,n):=DIF_NUMERIC(y,x,x0,h,n)
  7.  
  8. DIF2_NUM(y,x,x0,h):=DIF_NUMERIC(y,x,x0,h,2)
  9.  
  10. SMOOTH_VECTOR(v):=VECTOR(IF(m_=1 OR m_=DIMENSION(v),v SUB m_,(v SUB (m_-1)+v S~
  11. UB m_+v SUB (m_+1))/3),m_,DIMENSION(v))
  12.  
  13. SMOOTH_COLUMN_AUX(a,j):=VECTOR(IF(m_=j,SMOOTH_VECTOR(a SUB j),a SUB m_),m_,DIM~
  14. ENSION(a))
  15.  
  16. SMOOTH_COLUMN(a,j):=SMOOTH_COLUMN_AUX(a`,j)`
  17.  
  18. DIF_DATA_AUX(v1,v2):=(v1 SUB 1^2*v2 SUB 2-v1 SUB 2*v2 SUB 1^2)/(v1 SUB 1*v2 SU~
  19. B 1*(v1 SUB 1-v2 SUB 1))
  20.  
  21. DIF_DATA(a):=VECTOR([a SUB m_ SUB 1,IF(m_=1,DIF_DATA_AUX(a SUB 2-a SUB 1,a SUB~
  22.  3-a SUB 1),IF(m_=DIMENSION(a),DIF_DATA_AUX(a SUB (m_-1)-a SUB m_,a SUB (m_-2)~
  23. -a SUB m_),DIF_DATA_AUX(a SUB (m_+1)-a SUB m_,a SUB (m_-1)-a SUB m_)))],m_,DIM~
  24. ENSION(a))
  25.  
  26. DIF2_DATA_AUX(v1,v2):=2*(v1 SUB 1*v2 SUB 2-v1 SUB 2*v2 SUB 1)/(v1 SUB 1*v2 SUB~
  27.  1*(v2 SUB 1-v1 SUB 1))
  28.  
  29. DIF2_DATA(a):=VECTOR([a SUB m SUB 1,IF(m=1,DIF2_DATA_AUX(a SUB 3-a SUB 2,a SUB~
  30.  1-a SUB 2),IF(m=DIMENSION(a),DIF2_DATA_AUX(a SUB m-a SUB (m-1),a SUB (m-2)-a ~
  31. SUB (m-1)),DIF2_DATA_AUX(a SUB (m+1)-a SUB m,a SUB (m-1)-a SUB m)))],m,DIMENSI~
  32. ON(a))
  33.  
  34. INT_DATA(a):=[ELEMENT(a`,1),ELEMENT(ITERATES([1+ELEMENT(v,1),ELEMENT(v,2)+(ELE~
  35. MENT(a,ELEMENT(v,1),2)+ELEMENT(a,1+ELEMENT(v,1),2))*(ELEMENT(a,1+ELEMENT(v,1),~
  36. 1)-ELEMENT(a,ELEMENT(v,1),1))/2],v,[1,0],DIMENSION(a)-1)`,2)]`
  37.