Čo dokazuje tento graf na simulácií? Len toľko, že energia sa musí zachovať aj v nekonečne. Veď si predstavme simulované nekonečno chemickej reakcie, kde sa koncentrácia nukleárnej reakcie mení do okolia. Ak je teda toto pomyselné nekonečno, pýtajme sa inak, ako na obrázku. Á vidíme dotykové body hrubých čiar, ktoré toto potvrdili pre nepružný zráz aj pre delta en, či pomyselnú chemickú reakciu ^^^

using Pkg
Pkg.add("Plots")
using Plots
# konštanty
m1 = -1
v1 = 3.0
# os x
v2 = range(0.1, 10, length=10)  # vyhneme sa deleniu nulou
# funkcia m2(v2)
m2 = @. -m1 * (v1 + v2) / (v1 + 2v2)
# pôvodný graf (napr. Δₙ = 1)
n = 1:10
Δ = ones(length(n))
plot(
    n,
    Δ,
    label = "Δₙ=x2-x1",
    xlabel = "value x",
    ylabel = "value y",
    title = "Δₙ including inellastic collison, negative v2 for simmulating of infinity loop around collision velocity",
    size = (1920, 1080),
    guidefontsize = 12,   # osi (x, y)
    tickfontsize  = 12,   # čísla na osi
    legendfontsize = 14,  # legenda
    titlefontsize = 20,
    linewidth = 7
)
# pridanie tvojej rovnice
plot!(
    v2,
    m2,
    label = "m₂(v₂) negative m1",
    linewidth = 2
)

annotate!(
    8.0, -1.3,              # (x, y) v dátových jednotkách
    text("We know  the equation  of uncertainty  principle  is  
    delta x* delta p is greater  than or equal  to  h/ 2pi! 
    If Delta x* Delta p is equal  to  zero  then  
    the  particle  will be at rest!
    But  what  will  happen  
    if Delta x* Delta p  is equal  to  infinity?
    ANSWER: LOOK AT ENERGY, IT IS PRESERVED AT TOUCHING 
    POINTS OF THICK FIGURES! INFINITY DOES NOT MATTER!

    https://hrubos.tech/blogy/ (Mgr. Ján Hruboš)", 14)
)

m1 = 1
v1 = 3.0

# os x
v2 = range(0.1, 10, length=10)  # vyhneme sa deleniu nulou

# funkcia m2(v2)
m2 = @. -m1 * (v1 + v2) / (v1 + 2v2)

plot!(
    v2,
    m2,
    label = "m₂(v₂) positive m1",
    linewidth = 2
)


m1 = 1
v1 = 3.0

# os x
v2 = range(0.1, 10, length=10)  # vyhneme sa deleniu nulou

# funkcia m2(v2)
m2 = @. -m1 * (v1 - v2) / (v1 - 2v2)

plot!(
    v2,
    m2,
    label = "m₂(v₂) negative v2",
    linewidth = 2
)


m1 = 1
v1 = 1.0

# os x
v2 = range(0.1, 10, length=10)  # vyhneme sa deleniu nulou

# funkcia m2(v2)
m2 = @. -m1 * (v1 - v2) / (v1 - 2v2)
m = m1 .+ m2
v = v1 .+ v2
plot!(
    v,
    m,
    label = "m=m1+m2(v) NEGATIVE V2",
    linewidth = 7
)


m1 = 1
v1 = 3.0

# os x
v2 = range(0.1, 10, length=10)  # vyhneme sa deleniu nulou

# funkcia m2(v2)
m2 = @. -m1 * (v1 + v2) / (v1 + 2v2)
m = m1 .+ m2
v = v1 .+ v2
plot!(
    v,
    m,
    label = "m=m1+m2(v) positive v2",
    linewidth = 2
)

#porovnanie s koncentráciou v produkte chémie

# čas
t = range(0, 10, length=500)

# parametre reakcie
A = 1.8        # intenzita reakcie
k = 0.8        # rýchlostná konštanta
c∞ = 0.2       # ustálená koncentrácia

# koncentrácia produktu
c = @. A * t * exp(-k * t) + c∞

plot!(
    t,
    c,
    linewidth = 7,
    label = "chemical concentration reaction",
    guidefontsize = 18,
    tickfontsize = 14,
    legendfontsize = 14,
    titlefontsize = 20
)

# uloženie
savefig("delta7.png")
println("Figure saved as deltan.png")

---