# 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")