Infinity Time

x(t) = (x(t)^2 + 1)/2
with * (0) = 3

from *** import *

end_time = 42.
h = 0.0001
num_steps = int(end_time / h)
times = h * numpy.array(range(num_steps) + 1)

x = numpy.zeros(num_steps + 1)
x[0] = 3.

def forward_euler():
	for step in range(num_steps):

	return x

x = forward_euler()

@show_plot
def plot_me():
	matplotlib.pyplot.plot(times, x)
	matplotlib.pyplot.show()

plot_me()

Sources of Errors
-model, parameters, initial data, numerics