The Azimuth Project
Experiments in predator-prey with Sage (Rev #4, changes)

Showing changes from revision #3 to #4: Added | Removed | Changed



Showing an easy way of doing predator-prey modeling in Sage. Right now it is a generic version.Sage. Right now it is a generic version

Plot and Code


# Predator prey in Sage. default for ode_solver runga-kutta-felhberg (4,5)
T = ode_solver()
T.function = lambda t, y: [y[0]-y[0]*y[1], -y[1]+y[0]*y[1]]
sol_lines = Graphics()
for i in srange(0.1,1.1,.1):
    y = T.solution
    sol_lines = sol_lines + line([x[1] for x in y], rgbcolor = (i,0,1-i))
show(sol_lines+point((1,1),rgbcolor=(0,0,0)), figsize = [6,6], xmax = 6, ymax = 6)