{ diff(x1, t) = (1*k5*k1 + (-1)*k5*k2*x1^(k4 + 1)/(k3^k4 + x1^k4))/k5 }