blob: e4d00d1c6e6123f70706a417a077daeefb45e7db (
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
27
28
|
# Trivial Winston plotting example
using Winston
a=20
b=2
c=5
dt = 0.05
tlast = 2
iterations = int( 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
time = dt * [0:iterations-1]'
plot(time,xall)
title("Horray!")
xlabel("time")
ylabel("y")
|
