A = csv("a30");G = csv("g30");
Rm(1) = 300;
Rx(1) = radiusx(A, G, 2, 10);

A = csv("a25");G = csv("g25");
Rm(2) = 250;
Rx(2) = radiusx(A, G, 4, 12);

A = csv("a20");G = csv("g20");
Rm(3) = 200;
Rx(3) = radiusx(A, G, 4, 12);

A = csv("a15");G = csv("g15");
Rm(4) = 150;
Rx(4) = radiusx(A, G, 5, 12);

A = csv("a10");G = csv("g10");
Rm(5) = 100;
Rx(5) = radiusx(A, G, 4, 12);

A = csv("a05");G = csv("g05");
Rm(6) = 50;
Rx(6) = radiusx(A, G, 3, 11);

A = csv("a00");G = csv("g00");
Rm(7) = 0;
Rx(7) = radiusx(A, G, 4, 10);

A = csv("a-05");G = csv("g-05");
Rm(8) = -50;
Rx(8) = radiusx(A, G, 3, 10);

A = csv("a-10");G = csv("g-10");
Rm(9) = -100;
Rx(9) = radiusx(A, G, 3, 8);

A = csv("a-15");G = csv("g-15");
Rm(10) = -150;
Rx(10) = radiusx(A, G, 3, 10);

Fx = polyfit(Rm, Rx, 1);
plot(Rm, Rx, 'b.')
hold
plot(Rm,Fx(1)*Rm+Fx(2),'r-')
hold
legend('location', 'northwest')
legend('',['x_0 = ' num2str(-Fx(2)/Fx(1))])
legend('boxoff')
