def approx(x,A,F,y0,hp):
	S = x + hp * A * A #10-10
	Y = x+hp*S*S #10-09
	return Y

def step(F,y0, h, p):
	for i in range(p):
		y0 = [approx(y0[0],A[0],F,y0,h/p), approx(y0[1],A[1],F,y0,h/p)]
		print y0

A = [0.,2.]	 #example Moore1965interval p. 97
y0 = [1.,1.]
F = 'EXPLICITLY_DEFINED'
p = 2
step(F,y0, 0.25, p)