A = (xb - xa) / T - F T / (2 m) B = xa x = F t^2 / (2 m) + A t + B v = d(x,t) L = m v^2 / 2 + F x S = defint(L, t, 0, T) "Verify equation (1)" check(x == F t^2 / (2 m) + (xb - xa) t / T - F T t / (2 m) + xa) check(d(x,t,t) == F / m) check(eval(x,t,0) == xa) check(eval(x,t,T) == xb) "ok" "Verify equation (2)" check(v == F t / m + (xb - xa) / T - F T / (2 m)) "ok" "Verify equation (3)" check(S == m (xb - xa)^2 / (2 T) + F T (xb + xa) / 2 - F^2 T^3 / (24 m)) "ok" Run