Implementirali nekaj osnovnih metod za reševanje diferencialnih enačb

euler.m

rk2.m

Potem smo rešili sistem Lotke-Volterra: 

preypredator.m

S klicem

y0=[110;10];

[y,t]=rk2('preypredator',y0,0,10,1000,0);

plot(y(1,:),y(2,:))

dobimo sliko 

Last modified: Thursday, 11 January 2024, 10:07 PM