source term in time2.m