EllipsinRotaatio.jl

# Ellipsit origokeskinen ja (4,3)-keskinen
# Origon ympäri ja (4,3)-ympäri rotatoidut ellipsit
# Juhani Kaukoranta 24.10.2022
using Plots # original and rotated ellipse aroung origo
t = range(0, 2pi, length=100); # parametrinen kulma, 100 arvoa
a,b = 3,2 ; # ellipsin puoliakselit
x0,y0 = 4,3 # ellipsin keskipiste
x1,y1 = a * cos.(t) , b * sin.(t); # origokeskinen ellipsi
p1 = plot(x1,y1,xlabel="origokeskinen",title="a=3,b=2");
x2,y2 = a * cos.(t).+ x0 , b * sin.(t) .+ y0; # (x0,y0)-keskinen ellipsi
p2 = plot(x2,y2,xlabel="(4,3)-keskinen",title="non-rotated");
x3,y3 = (x1-y1)/sqrt(2) , (x1+y1)/sqrt(2); # origokeskinen rotated
p3 = plot(x3,y3,xlabel="origokeskinen",title="rotated");
x4,y4 = x2 .- x0, y2 .- y0; # keskipisteestä lasketut etäisyydet
x5,y5 = (x4-y4) / sqrt(2) .+ x0 , (x4+y4)/sqrt(2) .+ y0; # rotatoitu (x0,y0)-keskinen ellipsi
p4 = plot(x5,y5,xlabel="(4,3)-keskinen",title="rotated (4,3)");
plot(p1,p2,p3,p4,aspect_ratio=:equal,layout = (1, 4), legend = false)