The update of the reference solution was necessary as the problem description did not fit the implementation. The current time in an implicit Euler context has to be t+1.