[Numpy-discussion] === I REALLY NEED YOUR HELP OUR PYTHON USERS ====

Happyman bahtiyor_zohidov@mail...
Thu Dec 20 12:32:01 CST 2012


  Hi Python users, 

First of all, Marry coming Cristmas!!! ALL THE BEST TO YOU AND YOUR FAMILY

I need solution of integration under trapz() rule:
There are following functions:
  
def  F1 (const1, x):
      """several calculations depending on bessel functions(mathematical functions) jn(), yv() exists in Python"""
       return a,b

def   F2(const1 ,const2, D) :

       """Several calculation process"""
        x = D / const2
        [a , b] = F1 ( const1, x)      # Where x - the same as in F1() function
        S= a*b                               # This is (a*b) just an example for simply explanation
       return S

def F3(D, R):
       
      """Here I also calculated some process. So:"""
       return arg1**arg3   # Just for example

def  Integrate_all(const1, const2, min1, step1, max1):
        
       R=arange(min1, max1, step1)      # This is for function "F3"
       D = arange ( 0.1, 7.0, 0.0001)
       M = zeros ( size(R) )
       for i in range(0,size(R)):
                 M [ i ] = integrate. trapz ( ( F2 ( const1, const2, D ) * F3 ( D ,R)) , x=D)
      return M

      const1=complex number, const2= float, 

The aim of the calculation is to use Integrate_all function for integration function above!!!!!!!

When I use those functions directly like one by one separately from python shell it works very accurately, BUT when I do it as shown above :
 ERROR OCCURS:    C:\calculation.py:194: RuntimeWarning: invalid value encountered in divide!!!!!! (I think this is occuring in F1()) --> bessel functions >>> jn(n,x) and yv(n,x)

HERE IS THE PROBLEM!!!!! I HAVE BEEN TRYING FOR MORE THAN 1.5 MONTHS, UNFORTUNATELY I AM LOOSING MY INTEREST, and TIME.

PLEEEEEAAAASEEEEEE HELP ME MY FRIENDS!!!!!!!!! 
                 
       


        
        
       
        
     
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/numpy-discussion/attachments/20121220/1c3ee622/attachment.html 


More information about the NumPy-Discussion mailing list