program vaplap
             implicit double precision(a-h,o-z)
              dimension x(4000),y(4000),z(4000),p(4000)
		
CCCCC nombres de points de quadratures		
		         n=495
ccccc  ouverture des fichiers de points et remplissage
               open(unit=10,file='xs.txt',status='unknown')
	           open(unit=20,file='ys.txt',status='unknown')
	           open(unit=30,file='z.txt',status='unknown')
	           open(unit=40,file='Gamma.txt',status='unknown')
		       do j=1,n
		       read(10,*) x(j)
			   read(20,*) y(j)
			   read(30,*) z(j)
			   read(40,*) p(j)
		        enddo
		     close(10)		 
             close(20)	  
             close(30)	
             close(40)	
cccccc deux exemples sur des fonctions régulieres			 
              s=0.d0
			  sol=27*(dexp(1.d0/3)-dexp(-1.d0/3))**3
			  s1=0.d0
			  sol1=8.d0*dsin(1.d0)**3
            do i=1,n
             s=s+p(i)*dexp((x(i)+y(i)+z(i))/3.d0)
			 s1=s1+p(i)*dcos(x(i))*dcos(y(i))*dcos(z(i))
            enddo
Comparaison sol approchées sol exactes			
          print*,s	,sol,abs(s-sol)/sol
          print*,s1	,sol1,abs(s1-sol1)/sol1	
          
                        end