(Subscript[ε, ref])=(((∑_(i=1)^N)((Subscript[z, i])-(OverBar[z]))((Subscript[y, i])-(OverBar[y])))/((∑_(i=1)^N)((((Subscript[z, i])-(OverBar[z])))^2))),