1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
using Plots
using QuadGK
using LaTeXStrings
Plots.default(fontfamily="Computer Modern",
titlefontsize=20,
xtickfontsize=12,
ytickfontsize=12,
legendfontsize=12)
V(x) = 1/2 * (1 - cos(x))
η = .1
plot(title=L"Steady state $\rho_{\eta}$ with forcing")
plot(title="")
for η ∈ [0, .5, 1, 2]
ρ(q) = exp(-V(q)) * quadgk(y -> exp(V(q + y) - η*y), 0, 2π)[1]
Z = quadgk(ρ, - π, π)[1]
plot!(q -> ρ(q) / Z, -π, π, label=L"\eta = " * string(η), lw=2)
end
plot!(xlims = (-π, π), grid=false, framestyle=:box)
savefig("figures/invariant_perturbed_ol.pdf")
|