home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #18 / NN_1992_18.iso / spool / sci / math / symbolic / 2220 < prev    next >
Encoding:
Internet Message Format  |  1992-08-17  |  1.9 KB
  1. Path: sparky!uunet!olivea!decwrl!mips!swrinde!cs.utexas.edu!torn!watserv2.uwaterloo.ca!watdragon.uwaterloo.ca!watdragon!drclark
  2. From: drclark@daisy.uwaterloo.ca (David R. Clark)
  3. Newsgroups: sci.math.symbolic
  4. Subject: Re: Maple question.
  5. Message-ID: <DRCLARK.92Aug17085123@daisy.uwaterloo.ca>
  6. Date: 17 Aug 92 13:51:23 GMT
  7. References: <pete.714023792@marsh>
  8. Sender: news@watdragon.uwaterloo.ca (USENET News System)
  9. Organization: University of Waterloo, Waterloo, Ontario, Canada
  10. Lines: 33
  11. In-Reply-To: pete@marsh.cs.curtin.edu.au's message of Mon, 17 Aug 1992 03:56:32 GMT
  12.  
  13.  
  14.    Newsgroups: sci.math.symbolic
  15.    From: pete@marsh.cs.curtin.edu.au (Pete Staples)
  16.    Subject: Maple question.
  17.  
  18.    ...
  19.  
  20.    I have written a procedure that uses Laplaces' method to solve integrals.
  21.    Maple gives me an integral representation of the answer if I 
  22.    enter the limits (s=0..infinity) in the int procedure. If I evaluate
  23.    the indefinite integral and use subs for each limit of integration
  24.    it still doesn't evaluate the answer, eg leaves it in terms of 
  25.    erf(0), erf(x*infinity). Is there anyway that x > 0 can be assumed 
  26.    in the expression, so expressions like erf(x*infinity) can be interpreted?
  27.  
  28. There are a couple of problems here.  The first is that "subs" does not 
  29. attempt to evaluate an expression after performing the substitution.  The
  30. "eval" function can be used to force this evaluation.  This will force
  31. "erf(0)" to 0.  To compute the other endpoint you will have to use the
  32. "limit" function as it will handle erf(x*infinity) as x goes to infinity.
  33. This may well be what you want anyway and is the method used by the
  34. integration routines.
  35.  
  36. In general, assuming x to be greater than 0 can be done with an assignment
  37. to the signum remember table "signum(x) := 1;" but it will not help in this
  38. case as erf(x*infinity) will still not evaluate.
  39.  
  40.    Hope this makes sense, thanks in advance..
  41.  
  42.    Peter.
  43.  
  44.  
  45. David Clark
  46.