Advanced Search

Question: Why is my procedure refusing to accept that this variable equals zero!

Question:Why is my procedure refusing to accept that this variable equals zero!

brian abraham 599 Maple
plot
July 07 2007
0
In the attached file I have a variable RT_FR which appears in the system equation a1. Download 2129_forMaple(RT_FR).mws
View file details If I do the following (mb,tr,Ic are parameters) >RT_FR:=0: > a1 := -(PSX*sin(s1)*tr-RT_FR+mb*tr*vHn*v1)/(mb*tr^2+Ic); a1 := -11.64824973 vHn v1 then the numerical procedure - RT_FR is declared global - which is called by odeplot behaves exactly as it should. However if I do the following > a1 := -(PSX*sin(s1)*tr-RT_FR+mb*tr*vHn*v1)/(mb*tr^2+Ic): >RT_FR:=0: > a1; -11.64824973 vHn v1 Even though a1 appears to have no RT_FR term and this is confirmed using lprint(a1) within the numerical solver, the solver seems to include the RT_FR term. Since the solver contains the command > RT_FR:=RTFR(v1): the wrong output is produced. But surely if RT_FR is set to zero that should not happen! What subtlety of understanding am I missing?
Please Wait...