%% Use below (this is now in radial cooordinates) for comparison with Python.
nu = 10^-15;
chi = 7/3;
b = nu^(1/3)*sqrt(chi/(chi-2));
e = sqrt(3/8*(chi-2));
T = Ternary_model(2, 'FRAP', [-1, b, 0.5, e, 1, 0.5], linspace(0.0, 400, 5000));
T.t = 0:0.1:20;
T.solve_tern_frap()