clbeta.f 1.3 KB

1234567891011121314151617181920212223242526272829303132333435363738
  1. *DECK CLBETA
  2. COMPLEX FUNCTION CLBETA (A, B)
  3. C***BEGIN PROLOGUE CLBETA
  4. C***PURPOSE Compute the natural logarithm of the complete Beta
  5. C function.
  6. C***LIBRARY SLATEC (FNLIB)
  7. C***CATEGORY C7B
  8. C***TYPE COMPLEX (ALBETA-S, DLBETA-D, CLBETA-C)
  9. C***KEYWORDS FNLIB, LOGARITHM OF THE COMPLETE BETA FUNCTION,
  10. C SPECIAL FUNCTIONS
  11. C***AUTHOR Fullerton, W., (LANL)
  12. C***DESCRIPTION
  13. C
  14. C CLBETA computes the natural log of the complex valued complete beta
  15. C function of complex parameters A and B. This is a preliminary version
  16. C which is not accurate.
  17. C
  18. C Input Parameters:
  19. C A complex and the real part of A positive
  20. C B complex and the real part of B positive
  21. C
  22. C***REFERENCES (NONE)
  23. C***ROUTINES CALLED CLNGAM, XERMSG
  24. C***REVISION HISTORY (YYMMDD)
  25. C 770701 DATE WRITTEN
  26. C 861211 REVISION DATE from Version 3.2
  27. C 891214 Prologue converted to Version 4.0 format. (BAB)
  28. C 900315 CALLs to XERROR changed to CALLs to XERMSG. (THJ)
  29. C***END PROLOGUE CLBETA
  30. COMPLEX A, B, CLNGAM
  31. C***FIRST EXECUTABLE STATEMENT CLBETA
  32. IF (REAL(A) .LE. 0.0 .OR. REAL(B) .LE. 0.0) CALL XERMSG ('SLATEC',
  33. + 'CLBETA', 'REAL PART OF BOTH ARGUMENTS MUST BE GT 0', 1, 2)
  34. C
  35. CLBETA = CLNGAM(A) + CLNGAM(B) - CLNGAM(A+B)
  36. C
  37. RETURN
  38. END