rss_cpp function

Function for computing sigma_hat_sq