Comparison of floating-point evaluation with subs and eval
> sag_time2 := evalf( subs( DOvals, tcrit ), 3 );
> evalf( subs( sag_time2, DOeqn2 ), 3 );
>
> sag_time3 := eval( tcrit, DOvals );
> eval( DOeqn2, sag_time3 );