y

tmax = 5.5
tmin: ((-0.1) * tmax)
ymax = 2
ymin = -0.2
f(x) = (x * (1 - x))
y0actual = 0.3
ny = 10
dy: (ymax - ymin) / ny
y0possibilities = Sequence((0, y0actual + i), i, (ceil((ymin - y0actual) / dy) * dy), ymax - y0actual, dy)
Y0 = Point(y0possibilities)
dt = 1
y1actual: y0actual + (f(y0actual) * dt)
y1possibilities = Sequence((dt, y1actual + i), i, (ceil((ymin - y1actual) / dy) * dy), ymax - y1actual, dy)
Y1 = Point(y1possibilities)
y2actual: y1actual + (f(y1actual) * dt)
y3actual: y2actual + (f(y2actual) * dt)
y4actual: y3actual + (f(y3actual) * dt)
y2possibilities = Sequence(((2 * dt), y2actual + i), i, (ceil((ymin - y2actual) / dy) * dy), ymax - y2actual, dy)
Y2 = Point(y2possibilities)
y3possibilities = Sequence(((3 * dt), y3actual + i), i, (ceil((ymin - y3actual) / dy) * dy), ymax - y3actual, dy)
Y3 = Point(y3possibilities)
y4possibilities = Sequence(((4 * dt), y4actual + i), i, (ceil((ymin - y4actual) / dy) * dy), ymax - y4actual, dy)
Y4 = Point(y4possibilities)
a = Segment(Y0, Y1)
b = Segment(Y1, Y2)
c = Segment(Y2, Y3)
d = Segment(Y3, Y4)
nsteps = 5
y5actual: y4actual + (f(y4actual) * dt)
y5possibilities = Sequence(((5 * dt), y5actual + i), i, (ceil((ymin - y5actual) / dy) * dy), ymax - y5actual, dy)
Y5 = Point(y5possibilities)
e = Segment(Y4, Y5)
ylabeldisplayed: "" + (LaTeX[ylabel]) + ""
y0: y(Y0)
y1: y(Y1)
y2: y(Y2)
y3: y(Y3)
y4: y(Y4)
y5: y(Y5)