added error calculation for the Harada with the step change depending of diff Ks and initial positions