function z= BiRationalObjectiveFunction(p, p0, mea, time,T)
 
z= sum(abs(birational4parameters(p,p0,time,T)-mea)); %+1e+12*abs((k*b*d/c)-(k1*b1*d1/c1));

end