Forums before death by AOL, social media and spammers... "We can't have nice things"
|    sci.math.symbolic    |    Symbolic algebra discussion    |    10,432 messages    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
|    Message 10,193 of 10,432    |
|    Richard Fateman to All    |
|    optimal integral results    |
|    19 Dec 21 12:01:24    |
      From: fateman@gmail.com              There's more than just "size" to determine if the result of an integration is       optimal. Indeed, size might be misleading.        Consider the two results for integration of x^n wrt x.              r1= x^(n+1)/(n+1)              and       r2= (x^(n+1)-1)/(n+1)              Clearly r2 has more characters, and by some measure is more complex.              r1-r2 is 1/(n+1), a constant wrt x. so we should realize both answers are       correct "within a constant".              Arguably r2 is better because the limit of r2 as n approaches -1 is ln(x).       The limit as r1 approaches -1 is indeterminate.              Looking at a very small sample of the Rubi tests, it seems that other measures       of optimality might be missed. For example, an answer that has a structure of       a series in sin(x), sin(2*x), sin(3*x) .... vs powers of sin(x), cos(x)...              Another aspect is whether the integration program has as its responsiblity the       reduction of the answer -- or is the user likely to apply a simplification       program (say, to canonicalize rational functions, expand polynomials or       factor them, apply        trigonometric simplifications, etc.)              I can imagine a situation in which one program returns an answer that includes       an integral of sin(x)/x [thereby failing the Rubi test] where another       returns Si(x), the sine integral, defined as the integral of sin(x)/x. Which       looks mostly like a        nomenclature issue. There is also the possibility that (I see this in       Maxima), the result is expressed in another form entirely, e.g. via       Incomplete Gamma functions. I suppose that it depends on what your favorite       numerical library contains, Si,        Gamma, Hypergeometric functions etc.              If it is not already considered, perhaps the simplification issue can be       separated out somewhat. For instance, the Rubi testing might consider a       result R from Maxima to be too complicated. What about trigsimp(R), or       ratsimp(R) or for that matter,        ratsimp(trigsimp(R))? There are a number of additional plausible routines.       I don't see an easy way of handling the possibilities though.        In fact, given g = diff ( integrate(f,x),x), if the programs are correct, we       know that f and g are equivalent. Yet computing simplify(f-g) to get 0, for       some simplification program, may be quite difficult.              This testing of Rubi is valuable --- I look forward to Rubi 5 for Maxima,        the system I most often use.       RJF              --- SoupGate-Win32 v1.05        * Origin: you cannot sedate... all the things you hate (1:229/2)    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
(c) 1994, bbs@darkrealms.ca