"Nasser M. Abbasi" schrieb:   
   >   
   > I have translated Timofeev to sympy, and I can run the integrals now   
   > via sympy.   
   >   
   > But there were many many changes needed. Not only that, I just found   
   > that sympy changes all exact numbers to floating point. For example:   
   >   
   > >python   
   > >>> from sympy import *   
   > >>> x=symbols('x')   
   > >>> integrate(1/4/sin(1/3*x),x)   
   >   
   > 0.375*log(cos(0.333333333333333*x) - 1) -   
   > 0.375*log(cos(0.333333333333333*x) + 1)   
   >   
   > The above is #4 in Timofeev. Notice how the output is all   
   > decimal. It turns out, one has to tell sympy not to do this   
   > by changing the input to   
   >   
   > >>> integrate(S('1/4')/sin(S('1/3')*x),x)   
   >   
   >      3*log(cos(x/3) - 1)/8 - 3*log(cos(x/3) + 1)/8   
   >   
   > So, unless there is a way to tell sympy not to change exact   
   > numbers to floating point number, I would have to do change   
   > each number n in the test file to become S('n')?  63,000   
   > integrals?  Not practical.   
   >   
   > Will try to ask if there is an option to sympy, so that it   
   > will do this automatically. This is very strange for a computer   
   > algebra system to do this.   
   >   
      
   Or perhaps run:   
      
   integrate(sympify("integrate(1/4/sin(1/3*x),x)"),x)   
      
   The definition of symbols may be unnecessary then.   
      
   Martin.   
      
   --- SoupGate-Win32 v1.05   
    * Origin: you cannot sedate... all the things you hate (1:229/2)