另外,我试了用如下语句,其中a为方程组的解。但发现有时a里的数变了,不是刚解完方程组后的值。如正常a[1]=y1(t)=2t;但有时变成25t^2=-25t等莫名的数了。请高手指点迷津,谢谢!
> a:=dsolve(……):
> for i from 1 to 6 do
> Leqnsol:=lhs(a):Reqnsol:=rhs(a):
> if evalb(Leqnsol=y1(t)) then x1:=Reqnsol
> elif evalb(Leqnsol=y2(t)) then y1:=Reqnsol
> elif evalb(Leqnsol=y3(t)) then phi1:=Reqnsol
> elif evalb(Leqnsol=y4(t)) then x2:=Reqnsol
> elif evalb(Leqnsol=y5(t)) then y2:=Reqnsol
> else dphi2:=Reqnsol
> fi
> od;