http://www.mapleprimes.com/posts/42145-Solve-A-Problem en-us 2026 Maplesoft, A Division of Waterloo Maple Inc. Maplesoft Document System Fri, 12 Jun 2026 04:25:51 GMT Fri, 12 Jun 2026 04:25:51 GMT The latest comments added to the Post, solve a problem http://www.mapleprimes.com/images/mapleprimeswhite.jpg http://www.mapleprimes.com/posts/42145-Solve-A-Problem http://www.mapleprimes.com/posts/42145-Solve-A-Problem?ref=Feed:MaplePrimes:solve a problem:Comments#comment78767 First of all I think you should change <em>omega</em> variable to sth. else, e.g. <em>om</em>, since Maple translates <em>omega</em> as a built-in function [except, of course, you need to use <em>omega</em>]. You may change the code writing: restart: eq[1]:=add((-1)^(i+1)*cos(om*t[i]),i=1..5): eq[2]:=add((-1)^(i+1)*sin(om*t[i]),i=1..5): eq[3]:=add((-1)^(i+1)*cos(r*om*t[i]),i=1..5): eq[4]:=add((-1)^(i+1)*sin(r*om*t[i]),i=1..5): eq[5]:=t[1]: eq[6]:=t[4]-2*t[3]+t[2]: eq[7]:=t[5]-2*t[3]: s1:=fsolve({seq(eq[i],i=1..7)},{seq(t[i],i=1..5),r,om}); Solutions:=subs(s1,[seq(t[i],i=1..5), r, om]): t:=[seq(op(i,Solutions),i=1..5)]; r:=op(6,Solutions); om:=op(7,Solutions); <blockquote> You can check results by writing: seq(eq[i],i=1..7); pantole The latest comments added to the Post, solve a problem 78767 Wed, 13 Dec 2006 11:59:29 Z pantole pantole http://www.mapleprimes.com/posts/42145-Solve-A-Problem?ref=Feed:MaplePrimes:solve a problem:Comments#comment86092 Maple 10.06 should not confuse omega with the built-in function Wrightomega, unless it has been given as an alias for the latter. AFAIK there were no omega functions in previous versions of Maple. Hope this helps, J. Tarr The latest comments added to the Post, solve a problem 86092 Wed, 13 Dec 2006 15:27:17 Z Mariner Mariner http://www.mapleprimes.com/posts/42145-Solve-A-Problem?ref=Feed:MaplePrimes:solve a problem:Comments#comment78766 Thank you for your numeric way. But, I try to get a solution symbolically. t[1],t[2],...t[5] are time variables. oemga, r are constants. Could you give me help? Thanks, The latest comments added to the Post, solve a problem 78766 Wed, 13 Dec 2006 17:25:11 Z sungyg sungyg http://www.mapleprimes.com/posts/42145-Solve-A-Problem?ref=Feed:MaplePrimes:solve a problem:Comments#comment78765 J. Tarr is right when he points that it should work with variable <em>omega</em> without problem. I checked it out, but another paradox emerged: although, fsolve gave results with both omega and alternative variable, when used with omega the calculations delayed a little. Again with time function: < restart: st:=time(): Digits:=10: eq[1]:=add((-1)^(i+1)*cos(omega*t[i]),i=1..5): eq[2]:=add((-1)^(i+1)*sin(omega*t[i]),i=1..5): eq[3]:=add((-1)^(i+1)*cos(r*omega*t[i]),i=1..5): eq[4]:=add((-1)^(i+1)*sin(r*omega*t[i]),i=1..5): eq[5]:=t[1]: eq[6]:=t[4]-2*t[3]+t[2]: eq[7]:=t[5]-2*t[3]: s1:=fsolve({seq(eq[i],i=1..7)},{seq(t[i],i=1..5),r,omega}): Solutions:=subs(s1,[seq(t[i],i=1..5), r, omega]): t:=[seq(op(i,Solutions),i=1..5)]; r:=op(6,Solutions); omega:=op(7,Solutions); seq(eq[i],i=1..7); time()-st; > If <em>omega</em> is substituted for <em>om</em> for instance, calculations do not delay (maybe this happens at random). Still, since Maple output generates the greek omega (something like a <em>w</em>) for <em>omega </em>variable, then I guess that one should be sceptical for using this variable. But, maybe he shouldn't. pantole The latest comments added to the Post, solve a problem 78765 Wed, 13 Dec 2006 17:52:32 Z pantole pantole http://www.mapleprimes.com/posts/42145-Solve-A-Problem?ref=Feed:MaplePrimes:solve a problem:Comments#comment86087 Thanks anyway, I am trying to solve it for time variables. Could you make some code for arbitrary constants of omega and r for time t[2] and t[3] solutions? Those should be in positive format due to time variables. I think that you have to use "solve" command. Thanks The latest comments added to the Post, solve a problem 86087 Fri, 15 Dec 2006 22:08:27 Z sungyg sungyg http://www.mapleprimes.com/posts/42145-Solve-A-Problem?ref=Feed:MaplePrimes:solve a problem:Comments#comment91064 If you know that the values of the variables fall within a certain range, you can insert this as option in fsolve - please see ?fsolve,details. Something like this: restart: Digits:=15: eq[1]:=add((-1)^(i+1)*cos(omega*t[i]),i=1..5)=0: eq[2]:=add((-1)^(i+1)*sin(omega*t[i]),i=1..5)=0: eq[3]:=add((-1)^(i+1)*cos(r*omega*t[i]),i=1..5)=0: eq[4]:=add((-1)^(i+1)*sin(r*omega*t[i]),i=1..5)=0: eq[5]:=t[1]=0: eq[6]:=t[4]-2*t[3]+t[2]=0: eq[7]:=t[5]-2*t[3]=0: s1:=fsolve({seq(eq[i],i=1..7)},{seq(t[i],i=1..5),r,omega},0..4): Solutions:= [seq(t[i]=eval(t[i],s1),i=1..5),r=eval(r,s1),omega=eval(omega,s1)]; check_Solutions := seq(eval(eq[i],s1),i=1..7); Hope this helps, J. Tarr The latest comments added to the Post, solve a problem 91064 Sat, 16 Dec 2006 02:53:30 Z Mariner Mariner