Revision history [back]

Try this:

ans=desolve_system_rk4(des, [s,i,r], ics=[0, 1, 0.00000127,0], ivar=t, end_points=20)
ts=[[a[0],a[1]] for a in ans]
ti=[[a[0],a[2]] for a in ans]
tr=[[a[0],a[3]] for a in ans]
line(ts)+line(ti,color='red')+line(tr,color='green')