Program SIMPSON_1VAR
        IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        ! definimos el intervalo de integracion
        a=0d0
        b=1d0
        ! hacemos la aproximacion a la integral de la
        ! funcion g(x1) definida mas abajo
        
        r=(b-a)/6d0*(g(a)+4d0*g((a+b)/2d0)+g(b))
        
        ! abrimos un fichero excel que nos ense¤e el resultado
        open(20,File="ejerc10simp1var.xls")
        ! escribimos el resultado en el fichero
        write(20,*)r
        ! cerramos el fichero
        close(20)
        END

        ! definimos la funcion que queremos integrar
        FUNCTION  g(x1)
        Implicit Double Precision (A-H,O-Z)
        g=exp(-(x1*x1))*cos(x1*x1)
        END