old code, for refference/main_n_player.jl

include("./ring_np_num.jl")


ϕ_res = 200
ϕ = LinRange(0,2π, ϕ_res+1)[1:end-1]

p_0 = (zeros(ϕ_res) .+ 0.1)
p_0 .-= cos.(1 * ϕ) * 0.03

#config = Config(10, 0.1, 1.5, 0.1, 0.3, 5, 1, π, 2, 0.001, 1000, 30)
config_small = Config_small(1, 200, 0, π, 2, 0.01, 1000, 30)
p_end = solve_time_evolution(p_0, ϕ, config_small)

plot_result([p_end; config_small])