Good day,<div><br></div><div>                I 'm a graduate student and I'm using molpro to write cube files for atoms, and then I've wrote a fortran code to read it, because i want to recover the radial solution.  But when I plot the cube diagonal it shows something like a gaussian function with the maximum near the nucleus position. I was expecting to see something like the Hydrogen radial solution with zero value in the nucleus position. Is there any transformation that I should carry out to recover the hydrogen radial solution or am I missing something? The fortran code is attached, as the input file for molpro.</div>
<div><br></div><div><br></div><div><br></div><div>PS. I've compiled the file using gfortran and to use it is just "./executable input" where input is the cube file.<br clear="all"><div><br></div>-- <br><div style="text-align:center">
<b><i style="font-family:garamond,serif"><span style="color:rgb(0,0,153)">Rene Felipe Keidel Spada</span><br><span style="color:rgb(0,0,153)">Físico por Formação</span><br style="color:rgb(0,0,153)"><span style="color:rgb(0,0,153)">Humano por Natureza</span><br>
<u style="color:rgb(0,0,153)">Instituto Tecnológico de Aeronáutica</u></i></b><br></div><br>
</div>