implicit real*8 (a-h,o-z)
        common/cnsts/pi,sq2,sq2pi,sq2pil
        external func
        call dflt
        open(unit=8,file='regfal.out')
        a=-3.
        b=3.
        iterl=100
        acc=1.d-10
        call regfal(a,b,sol,func,acc,iterl,ier)
        write (8,1000) sol
1000    format(' Solution = ',f12.8)
        stop
        end
        function func(x)
        implicit real*8 (a-h,o-z)
        common/cnsts/pi,sq2,sq2pi,sq2pil
        func=0.5+0.5*derf(x/sq2)
        return
        end
*************************************************************************
OUTPUT FILE FOLLOWS
*************************************************************************
 Solution =   0.00000000
 
<p align="center"> <font color="#400040">Return to</p>
<p align="center">|</font><a href="e.htm"><font color="#400040">Sect. E</font></a><font
color="#400040">|</font><a href="1.htm"><font
color="#400040">Beginning</font></a><font color="#400040">|</font></p>

</body>
</html>