Finite elements in time domain for optimal control using higher-order shape functions