<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def reg_score_function(X, y, mean, scale, shape, skewness):
 """ GAS Cauchy Regression Update term using gradient only - native Python function Parameters X : float datapoint for the right hand side variable y : float datapoint for the time series mean : float location parameter for the Cauchy distribution scale : float scale parameter for the Cauchy distribution shape : float tail thickness parameter for the Cauchy distribution skewness : float skewness parameter for the Cauchy distribution Returns - Score of the Cauchy family """

        return 2.0*((y-mean)*X)/(np.power(scale,2)+np.power((y-mean),2))