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