Adaptive quadrature using Filon's method

Please choose an integrand:

	g(x) = exp(-x2)*sin(π x)
	g(x) = cos(x π)2
	g(x) = 1/((x-1/2)^2+1/1000)+1/((x-7/8)^2+1/100)
	A function of your own: Please type the evaluation program of your function here using FORTRAN rules. Your final statement must be fu= some expression you computed before or just here depending on x You may use the constants pi, e(=exp(1)), sqrt2(=1.414...), the integer variables i,j,k, the logicals bool1,bool2,bool3 and the double precision variables sum,h1,h2,h3,h4,y(100),z(100),a(100,100) which are all intialized with zero resp. .false. . The routine has the parameters x (double, input) and fu (double out). never change x!. first is a local integer and set 0 before calling the function the first time. You may use this in order to initialize some local data and set it 1 afterwards to avoid multiple such initialization. Your settings of the local variables are preserved during program execution. C**** it computes the truncated McLaurin series of sin(pi*x) multiplied C**** by exp(-pi*x^2) if ( first .ne. 0 ) goto 100 C**** now compute the 40 first Taylor coefficients of sin(pi*x) in y(.) y(1)=pi do 50 i=2,40 y(i)=-y(i-1)*pi**2/dble(2*(i-1)*(2*i-1)) 50 continue 100 continue first=1 fu=0.d0 do 200 i=40,1,-1 fu=fu+y(i)*x**(i+i-1) 200 continue fu=fu*dexp(-pi*x**2)

Please specify the interval of integration (a < b):

a =

b =

We also need an upper bound for the stepsize:
h_max=

h_max <= (b-a)/10 is controlled internally .

Choose your desired precision (bound on the absolute error):
eps=
Important: 1.0e-12 <= eps <= 0.01

Choose ω > = 1
ω =

Click on "evaluate", in order to submit your input.

A message ''h reduced to h_min'' may be issued for a too small eps or if f is not smooth