def METHOD_NAME(self,t):
    uh = self.uh0
    def f(x):
        return self.pde.source(x, t)
    b = doperator.source_vector(f, self.femspace, self.integrator, self.area)
    return self.M@uh + self.tau*b