Hi,
I am currently trying to compute the hypergeometric function in Fortran on Microsoft Visual Studio.
I am using the FUNCTION hypgeo, but there seems to be a problem in the code as calling odeint requires y (the first input to odeint) to be a REAL number, while y is defined as COMPLEX in hypgeo.
My interpretation to the problem is that the input to odeint should be an array of y(1,4) with the following values: y(1)=REAL(y(1)); y(2)=IMAG(y(1)); y(3)=REAL(y(2)); y(4)=IMAG(y(2)) in the code, here y(1) and y(2) were computed by hypser, see code in book. Am I correct in altering the code this way?

Also, the book did not really specify what the 'right-hand side derivative' is when calling derivs in odeint. Can someone offer more explanation on what I should supply to the subroutine derivs when I am trying to us odeint in hypgeo?

Thank you very much.