
	subroutine integra(func,a,b,ris,npassi)
	ris=0
	dx=(b-a)/npassi
	xm=a+(dx/2.)
	do i=1,npassi-1
	ris=ris+func(xm)*dx
	xm=xm+dx
	end do
	return
	end
