Joulupipari.jl

# Juhani Kaukoranta 20.12.2022
# n-muhkurainen joulupipari
# f(x) = sin(n*x)/n + m
# plottaa (x1,y1) = (f(x)cos(x),f(x)sin(x))
# laskee pinta-alan ja reunan max/min-etäisyyden origosta
using Plots
function joulupipari(n,m) # piirtää piparin
    x = range(0, 2*pi, length=1000) # kulma radiaaneina
    f(x) = sin(n*x)/n+m
    x1 = @. cos(x)*f(x)
    y1 = @. sin(x)*f(x)
    A = (4*n*(2*pi*m^2*n^2+2*m+pi)-8*m*n*cos(2*pi*n)-sin(4*pi*n))/(8*n^3) #pinta-ala
    maxR =1/n + m # kuvaajan maksimietäisyys origosta
    minR = m # kuvaajan reunan minimietäisyys origosta
    A=round(A,sigdigits=5) # rajoitetaan tarkkuus 5 numeroon
    maxR,minR =round(maxR,sigdigits=5),round(minR,sigdigits=5)
    otsikko = "Kuvaajan max- ja min-etäisyydet origosta \n$maxR ja $minR  ja pinta-ala = $A\n(n,m)=($n,$m)"
    plot(x1,y1,aspect_ratio = :equal,title=otsikko,legend=false)
end

# esim arvoilla n = 7,m = 1 ajetaan komennolla joulupipari(7,1)