Subprogramas y Funciones

• Calcular la función de distribución radial (FDR) para los orbitales 1s, 2s y 2p y entonces normalizar el resultado al valor máximo de cada FDR. La FDR es el valor de la función de ondas al cuadrado multiplicado por 4r² 

• Las funciones de ondas son:

  1s=exp(-r/2)

  2s=(2-r)exp(-r/2)/sqrt(32)

  2p=(6-6r+r**2)exp(-r/2)/sqrt(972)