home *** CD-ROM | disk | FTP | other *** search
/ Chip 2002 June / Chip_2002-06_cd1.bin / zkuste / derive / download / Setup.exe / %MAINDIR% / Users / Int_misc.txt < prev    next >
Encoding:
Text File  |  2002-05-15  |  3.3 KB  |  50 lines

  1. Continuous antiderivatives and change of variable in integrals
  2. ==============================================================
  3.  
  4. Author:  Valeriu Anisiu
  5. email:   anisiu@math.ubbcluj.ro
  6. address: Faculty of Mathematics and Computer Science
  7.          "Babes-Bolyai" University
  8.          Kogalniceanu 1 Street
  9.          3400 Cluj-Napoca, Romania
  10.  
  11.  
  12. The file INT_MISC.MTH contains three user functions which can be used to help or correct the integration process in DERIVE.
  13.  
  14. MAKECONT(y,x,v) := adds a locally constant function to y such that it becomes continuous at the points contained in the vector v.
  15.  
  16. CONT_INT(y, x, v) := INT(y,x), but it is "guaranteed" to be continuous at the points contained in the vector v. This function is useful in those (rare) occasions when DERIVE does not find a continuous antiderivative.
  17.  
  18. CONT_INT(y, x, v, a, b) := INT(y, x, a, b); the result is "guaranteed"  to be correct even if the antiderivative found by DERIVE is not continuous at the points contained in the vector v.
  19.  
  20. For example, the antiderivative found by DERIVE for ASIN(2*x/(x^2+1)) is discontinuous at the points x=-1 and x=1.
  21. To find a continuous antiderivative, use CONT_INT(ASIN(2*x/(x^2+1)), x, [-1,1]).
  22. The correct integral for this function on the interval [0,2] is CONT_INT(ASIN(2*x/(x^2+1)), x, [-1,1], 0, 2) or just CONT_INT(ASIN(2*x/(x^2+1)), x, [1], 0, 2) (because -1 does not belong to the interval [0,2]).
  23.  
  24.  
  25. INTSUBST(y, x, rx, t, rt) := the antiderivative of y(x) computed using the change of variable given by the relation rx.
  26. rx is a relation (equality) between x and t which can be solved (by SOLVE) with respect to t.
  27.    If rx in not a relation but an expression depending on x, rx=t is considered.
  28.    t denotes the new (intermediate) variable, and rt is a relation equivalent to the relation rx (in a certain interval) which can be solved with respect to t. It was necesary to deal with two relations because the SOLVE function is not powerful enough.
  29.    If the relation rt is absent, rt:=rx.
  30.    As a bonus, the INTSUBST can apply on demand the well known Euler substitutions for y = R(x,SQRT(a*x^2 + b+x +c)), where R is a rational function: in this case rx must have one of the special forms
  31.    a*x^2 + b+x +c = "euler1",  a*x^2 + b+x +c = "euler2"  or  a*x^2 + b+x +c = "euler3"
  32. [As a recommandation, use "euler1" when a>0, "euler2" when c>0 and "euler3" when b^2-4ac>0].
  33.    Furthermore, if the second argument x is replaced with the vector [x,a,b], the definite integral is computed.
  34.    The user is responsible for the equivalence of the given relations in the interval where he/she wants to integrate.
  35.  
  36. As one can see, the INTSUBST is very versatile, and can be even extended for other changes of variable (in order to eliminate the computation by hand of the inverse substitution). Please inspect its definition in the file INT_MISC.MTH.
  37.  
  38. For example, if 
  39. f:= 1/(SQRT(x^2+x+4)+2*x+3)
  40.  
  41. then an antiderivative for f can be computed with:
  42.  
  43. INTSUBST(f, x, x^2+x+4="euler1")        or
  44. INTSUBST(f, x, x^2+x+4="euler2")        or
  45. INTSUBST(f, x, SQRT(x^2+x+4)=x+t, t, x^2+x+4=(x+t)^2).
  46. Note that the integral cannot be compuded with INT or INT_SUBST.
  47.  
  48. INTSUBST(EXP(ASIN(x)),[x,0,1],ASIN(x)) computes the integral of EXP(ASIN(x)) between x=0 and x=1. 
  49.  
  50. For other examples, please load and run the demo file INT_MISC.DMO.