include("../NuclearRMF.jl")

Z = 82
N = 126
r_max = 20.0
mesh_points = 400

s = system(Z, N, mesh_points, r_max)

solve_system!(s; monitor_plot=true)