Hi,
I have a system of differential equations and sage give me the
solution as follows,
The system:
x10,x20,x30,t,f1,f2,c11,c22,c12,c23,k1,k2,x1,x2,x3=var('x10,x20,x30,t,f1,f2,c11,c22,c12,c23,k1,k2,x1,x2,x3')
v1=f1*(x1-k1*x2)/(1+c11*x1+c12*x2)
v2=f2*(x2-k2*x3)/(1+c22*x2+c23*x3)
x1= function('x1',t)
x2= function('x2',t)
x3= function('x3',t)
d1=(diff(x1,t)==0)
d3=(diff(x3,t)==0)
d2=(diff(x2,t)-v1+v2==0)
sol=desolve_system([d1,d2,d3],[x1,x2,x3],[x10,x20,x30],t);show(sol)
The result
[x1(t)=x1(0);x2(t)=¡(c22f1k1+c12f2)tx22c12c22x22+(c11c23x1+c12c23x2+c23)x3+(c11c22x1+c12+c22)x2+c11x1+1+((c12f2k2¡c23f1k1)x2+(c11f2k2+c23f1)x1+f2k2)tx3c12c22x22+(c11c23x1+c12c23x2+c23)x3+(c11c22x1+c12+c22)x2+c11x1+1¡((c11f2¡c22f1)x1+f1k1+f2)tx2c12c22x22+(c11c23x1+c12c23x2+c23)x3+(c11c22x1+c12+c22)x2+c11x1+1+f1tx1c12c22x22+(c11c23x1+c12c23x2+c23)x3+(c11c22x1+c12+c22)x2+c11x1+1+x2(0);x3(t)=x3(0)]
I need to use the value like I need to take nominator of the component
of the solution but I could not.
Any suggestion
Thanks
Doaa
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
