program vaplap
             implicit double precision(a-h,o-z)
              dimension x(4000),y(4000),z(4000),t(4000),p(4000)
		
CCCCC nombres de points de quadratures		
		         n=1512
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='zs.txt',status='unknown')
			   open(unit=40,file='ts.txt',status='unknown')
	           open(unit=50,file='Gamma.txt',status='unknown')
		       do j=1,n
		       read(10,*) x(j)
			   read(20,*) y(j)
			   read(30,*) z(j)
			   read(40,*) t(j)
			   read(50,*) p(j)
		        enddo
		       close(10)		 
               close(20)	  
               close(30)	
               close(40)	
			   close(50)
cccccc deux exemples sur des fonctions régulieres			 
              s=0.d0
			  sol=256*(dexp(1.d0/4)-dexp(-1.d0/4))**4
			  s1=0.d0
			  sol1=16.d0*dsin(1.d0)**4
            do i=1,n
             s=s+p(i)*dexp((x(i)+y(i)+z(i)+t(i))/4.d0)
			 s1=s1+p(i)*dcos(x(i))*dcos(y(i))*dcos(z(i))*dcos(t(i))
            enddo
Comparaison sol approchées sol exactes			
          print*,s	,sol,abs(sol-s)/s
          print*,s1	,sol1,abs(s1-sol1)/s1	
          
                        end