3649.sagews

var('x, y')

(x, y)

A = matrix([[2, 1], [-1, 2]])
x0 = vector([0.1, 0.01])

A*x0

(0.210000000000000, -0.0800000000000000)

P = point(x0)
x_last = x0
for i in range(5):
    x_next = A*x_last
    P += point(x_next)
    x_last = x_next
P.show()

P = point(x0)
x_last = x0
for i in range(5):
    x_next = A*x_last
    P += line([x_last,x_next])
    x_last = x_next
P.show()