clear all
close all
[t1, y1]=expeuler_sajat(@zh21de,0,2,0.5,8);
[t2, y2]=expeuler_sajat(@zh21de,0,2,0.001,4000);
u=2*exp(-3*t2);
subplot(2,1,1)
plot(t1,y1,':',t2,u,'r')
xlabel('t')
ylabel('y')
title('Megoldás h = 0.5 lépésközzel')
legend('Explicit Euler-módszer','Pontos megoldás')
subplot(2,1,2)
plot(t2,y2,':',t2,u,'r')
xlabel('t')
ylabel('y')
title('Megoldás h = 0.001 lépésközzel')
legend('Explicit Euler-módszer','Pontos megoldás')