f'=1 init f=0 u(t)=int[.5]{1#f} aux ut=(4/3)*t*sqrt(t) done