function v=rho_prime(x)
global DATAPOINTS
global P_VALUE

z = DATAPOINTS;
p = P_VALUE;
xx = x'*ones(size(z));
zz = ones(size(x))'*z;
diff = xx-zz;
diff2 = abs(diff).^(p-1).*sign(diff);
v = sum(diff2')';
