```
using Plots
a = [1,3,5,3]
b = [1,5,7,4]
p = []
δ = π/2
t = LinRange(-π, π, 300)
for i ∈ 1:4
x = sin.(a[i] .* t .+ δ)
y = sin.(b[i] .* t)
push!(p, plot(x, y, legend = false, aspect_ratio= :equal))
end
plot(p..., layout=4)
```

```
┌ Info: Precompiling Plots [91a5bcdd-55d7-5caf-9e0b-520d859cae80]
└ @ Base loading.jl:1664
```

```
┌ Info: Precompiling GR_jll [d2c73de3-f751-5644-a686-071e5b155ba9]
└ @ Base loading.jl:1664
```