computes Mscale using Newton Raphson

computes Mscale using Newton Raphson

MscaleNR(u, c1, k1)

Arguments

u	vector of distances (positive numbers)
c1	constant for rho_1 function
k1	constant for rho_1 function