Defined SetupRHS

this function will set up RHS of the linear system. But only defined for
first time step.