 function r = reseul(par, dt, z)
a = par(1);
b = par(2);
c=par(3);
g=par(4);
dz=diff(z);
n = length(dz);
z1 = z(1:n);
   a0 = a + b.*z1;
s = c*z1.^g;
r = (dz - dt*a0)./s;
