y=@(u,t)-(5138804207313433*u*exp((1727*t)/(100*(t + 2373/10))))/(14073748835532800*(t +... 273)*((1411*u)/62500 + (25031*exp((1727*t)/(100*(t + 2373/10))))/(10*(t + 2373/10)^2) + 83/1250));>> y(1,5)ans = -0.0125