function [e] = StepError(s0,par)
delta_s = Step(s0,0,par)-s0;
e = sum(delta_s.^2);