Program MONTECARLO_HIT_AND_MISS_10var
 		Implicit Double precision (A-H,o-z)
 		DOUBLE PRECISION g ,a,b,c,p,r,s ,u,v,pi
 		integer n , na

        a=0d0
        b=1d0
        c=1d0
        na=0
        n=10
        pi=acos(-1d0)
             open(60,File="ejerc10aciertfallos10var.xls")
               call cpu_time(start)

                do j=n,10000000000

        	do 1 i=1,n
		u1=rand()
                v=rand()
                u2=rand()
		u3=rand()
		u4=rand()
		u5=rand()
		u6=rand()
		u7=rand()
		u8=rand()
		u9=rand()
		u10=rand()
		
	if (g(a+(b-a)*u1,a+(b-a)*u2,a+(b-a)*u3,a+(b-a)*u4,a+(b-a)*u5,
     &  a+(b-a)*u6,a+(b-a)*u7,a+(b-a)*u8,a+(b-a)*u9,a+(b-a)*u10).gt.c*v)
     &  na=na+1

1 	        continue
		p= real (na)/n
		r=(b-a)*c*p
		s=sqrt(p*(1.-p)/n)*c*(b-a)
                write(60,20)n,r,s
                if (n.gt.1d7) go to 30
                n=n*10
                na=0
                enddo
30              call cpu_time(finish)
                write(60,*)"Temps = ",finish-start ,"s"
                close(60)
20              Format(I20,4X,F20.14,4X,E20.14)
                end

		Function g(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10)
                Implicit Double Precision (A-H,O-Z)
                g=exp(-x1**2-x2**2-x3**2-x4**2-x5**2-x6**2-x7**2-x8**2-
     &          x9**2-x10**2)*cos(x1*x2+x2*x3+x3*x4+x4*x5+x5*x6+x6*x7+x7
     &          *x8+x8*x9+x9*x10+x10*x1)
                end