residY_original function

Compute the residuals from the observational equation