w = [3.5 4 6 25 103 117 200 252 300 437 1340 2000 2700 5000 22500 30000 33000 50000 70000 100000 415000 450000 500000 3000000];
r = [787 660 588 670 347 300 420 352 300 269 251 205 187 120 100 85 81 70 72 70 45 38 40 48];



x = w.^(-1/3)'
y = r'

[m,b] = ls_normal (x,y)