A 5th-order method for 1D-device solution