clear;

V = 150;
ETA = linspace(1.60,1.70);

A = (2 * V ./ (sqrt(3) * ETA)) .^ (1 / 3);

RMT1 = A ./ 2;
RMT2 = A ./ 2 .* sqrt(1/3 + (ETA .^ 2) ./ 4);

plot(ETA, RMT1, '-r');
plot(ETA, RMT2, '--b');
xlabel("c/a");
ylabel("MT radius (Bohr)");