I have changed the method used for calculating the w.f. by solving the
RK4 method several times within a r step saving only the last part
improving the accuracy of the method and removing anomalous behaviour
seen for negative Z values.