module Quanc8FMM
!!=============================================================================
!! Procedura Quanc8 
!! "Computer Methods for Mathematical  Computations",
!! by George E. Forsythe, Michael A. Malcolm, and Cleve B. Moler
!! Prentice-Hall, 1977.
!!-----------------------------------------------------------------------------
!! Prepis pro Fortran 90:
!! Last updated: 28 February 2005 by Ralph Carmichael ralph@pdas.com 
!! Public Domain Aeronautical Software 
!! P.O. Box 1438 Santa Cruz CA 95061 USA 
!!-----------------------------------------------------------------------------
!! Drobne upravy pro g95
!! Jan Cely, UFKL,Prirodovedecka fakulta MU, 270307,130208
!!=============================================================================
!!
!! SUBROUTINE Quanc8 (Fun,a,b,abserr,relerr,hodnota,errest,nofun,flag)
!!
!!  Estimate the integral of FUN(X) from A to B to a user provided 
!!  tolerance. This is an automatic adaptive routine based on the 8-panel 
!!  Newton-Cotes rule. 
!!	VSTUP:
!!	Fun		= integrand deklarovany takto : FUNCTION Fun(x) RESULT(f) 
!!	a,b		= integracni meze
!!	abserr	= maximum pro absolutni chybu (>=0)
!!	relerr	= hranice pro relativni chybu (>=0)
!!	VYSTUP:
!!	hodnota	= vypocteny odhad integralu
!!	errest	= odhad velikosti skutecne absoltni chyby
!!	flag	= informace o spolehlivosti vysledku
!!			  0.0 = vysledek pravdepodobne vyhovuje zadanym pozadavkum na chybu
!!			  XXX.YYY = cast zakladniho intervalu, ktera zustala nezpracovana,
!!			  kdyz se procedura priblizila k limitni hodnote pro "nofun"
!!			  Zadano v procedure jako INTEGER,PARAMETER :: nomax = 5000
!!	nofun	= pocet vypoctu Fun pri behu procedury
!!=============================================================================
!