Implement time-dependant sources in solve()