blob: 1773e6cca666d63d11a38e343bb9af1271eb812e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
% % Use Euler's method to integrate simple one variable ODE
% % Use program as template for more interesting models
% % Note: produces slightly different output compared with
% % example shown in class. Different values of b, dt, tlast
a = 20 ;
b = 0.5 ;
c = 5 ;
dt = 0.2 ;
tlast = 20 ; % s
iterations = round(tlast/dt) ;
xall = zeros(iterations,1) ;
x = c ;
for i = 1:iterations
xall(i) = x ;
dxdt = a - b*x ;
x = x + dxdt*dt ;
end % of this time step
time = dt*(0:iterations-1)' ;
figure
plot(time,xall)
|