%% 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()