# Kezdo parancsok minden script elejen

# Workspace kiuritese, valtozok torlese
clear

# Megnyitott abrak becsukasa
close all

display("Lehules modellezese, kulso homerseklet 17 fok, lehulesi egyutthato 1.2")
f = @(y, t) - 1.2 * (y - 17)

display("Egzakt megoldas:")
g = @(y0, t) (y0 - 17) * exp(-1.2 * t) + 17

# idoracs
t = linspace(0, 10, 10);

plot(t, g(35, t))

# hold on - megtartja aktivba a grafikont es lehet ra mast rajzolni, akar mas idoskalan is
hold on

t = linspace(0, 10, 1000);
plot(t, g(35, t))
title("Lehules egzakt megoldasai, 17fok kornyezet, 1.2-es wegyutthato, 35fokos kezdeti homerseklet")
legend("Egzakt 10 pontban", "Egzakt 1000 pontban")