#EJEMPLO 1


from sage.calculus.desolvers import desolve_system_rk4
sage: A,B,t=var('A B t')
sage: P=desolve_system_rk4([120-(A/25),(A/25)-2*B/(2*t+100)],[A,B],ics=[0,0,0],ivar=t,end_points=200) #ics[tiempo,indiceSalinidad T1,indiceSalinidad T2]
sage: Q=[ [i,j] for i,j,k in P]
sage: LP=list_plot(Q)
sage: Q1=[ [i,k] for i,j,k in P]
sage: LP1=list_plot(Q1,color='orange')
show(LP+LP1)

from sage.calculus.desolvers import desolve_system_rk4
sage: A,B,t=var('A B t')
sage: P=desolve_system_rk4([0.2-(A/10),(A/10)-3*B/(3-2.8*t)],[A,B],ics=[0,0,0],ivar=t,end_points=20) #ics[tiempo,indiceSalinidad T1,indiceSalinidad T2]
sage: Q=[ [i,j] for i,j,k in P]
sage: LP=list_plot(Q)
sage: Q1=[ [i,k] for i,j,k in P]
sage: LP1=list_plot(Q1,color='green')
show(LP+LP1)

from sage.calculus.desolvers import desolve_system_rk4
sage: A,B,t=var('A B t')
sage: P=desolve_system_rk4([50-(A),(A)-2*B],[A,B],ics=[0,0,0],ivar=t,end_points=7,step=0.1) #ics[tiempo,indiceSalinidad TA,indiceSalinidad TB]
sage: Q=[ [i,j] for i,j,k in P]
sage: LP=list_plot(Q)
sage: Q1=[ [i,k] for i,j,k in P]
sage: LP1=list_plot(Q1,color='orange')
show(LP+LP1)

from sage.calculus.desolvers import desolve_system_rk4
sage: A,C,B,t=var('A C B t')
sage: P=desolve_system_rk4([10-(3*A/(2*t+100)),(3*A/(2*t+100))-3*C/100 , 3*C/100 - (5*B/(100-2*t)) ],[A,C,B],ics=[0,0,0,0],ivar=t,end_points=45) #ics[tiempo,indiceSalinidad TA,indiceSalinidad TC,indiceSalinidad TB]
sage: Q=[ [i,j] for i,j,k,W in P]
sage: LP=list_plot(Q)
sage: Q1=[ [i,k] for i,j,k,W in P]
sage: LP1=list_plot(Q1,color='red')
sage: Q2=[ [i,W] for i,j,k,W in P]
sage: LP2=list_plot(Q2,color='green')
show(LP+LP1+LP2)