compute_resid function

Compute residual