procedure problem...

i wrote a procedure for determine a lambda symmetries and first integrals of ODEs.

Procedure is running but commands is getting error.

Can  anybody help me.

unexpected string and procedure error...

i wrote a procedure to  find some functions but at the end  of the procedure gives unexpected string error. i didn't find any way out.

Best regards.

Suggestions on how to check the progress of a proc...

Hello

I wonder if there is a way to see the progress of a specific procedure.  I have checked the forum for some answers and found that some of them propose a progress bar. However the progress bar doesn't seem to fit within a procedure.

I don't need anything fancy but I need to see the progress of the loop inside of the procedure (over 3,000,000 times).

Any suggestions?

Many thanks

PS. I am running Maple 2016 too but I would rather see a solution that works on Maple 14 too, if possible.

How to retrieve only the monomials (no coefficient...

Hello

I wonder if there is a function to retrieve only the monomials from a multivariable polynomial in x, y and z. Below is one such polynomail.

theta[2]*x3*x0-theta[8]*x2*x1*x0-theta[2]*x2*x1-theta[6]*theta[2]*x2*x0^2-theta[4]*theta[2]*x2*x0-theta[2]*theta[7]*theta[10]*x1*x0^3+theta[8]*theta[2]*theta[9]*x1*x0^2-theta[2]*theta[5]*theta[10]*x1*x0^2-theta[2]^2*theta[3]*theta[10]*x0^3+theta[6]*theta[2]^2*theta[9]*x0^3+theta[4]*theta[2]^2*theta[9]*x0^2

and the result

[x3*x0,x2*x1*x0,x2*x0^2,x2*x0,x1*x0^3,x1*x0^2,x0^3,x0^2]

Many thanks

Ed

Plot in a same graph ( a sequence of functions in ...

Hello,

I have a sequence of functions :  solution[i] , i = 1..n

I have a sequence of times:          Time[i], i = 1..n+1

I need help to plot in the same graph:

plot(  Sol[1] , t = Time[1]..Time[2] )

plot(  Sol[2] , t = Time[2]..Time[3] )

plot(  Sol[3] , t = Time[3]..Time[4] )

etc...

plot(  Sol[n] , t = Time[n]..Time[n+1] )

Thank you

using worksheet mode, is there a way to automatically color any comments after (#) with a different color than the default red?

I have been using Maple for years and did not need to significantly document or comment on my worksheets before because I only needed to share my worksheets with close collegues who have a lot of experience with Maple too. However, now I need to share my code with a general audience that might not be a Maple user. So I need to add enough comments after each line, and I have been manually changing the color of comments after the # to green to give the reader the indication that this is a comment and not part of the code, like this for example:

 > restart;
 > f:=(x,y)->sin(sqrt(x^2+y^2))/sqrt(x^2+y^2); # Define f as a function of the variables x and y
 (1)
 > df:=(x,y)->eval(diff(f(a,y),a),a=x); # Define df as the partial derivative of f with respect to x
 (2)
 > df(1,3); # Evaluate df at x=1, y=3
 (3)
 >

It would be great if Maple can automatically color comments like typical editors do. Is this possible? if not, does anyone have an advice for a nice easy way to add proper documentation to worksheets?

Modeling of DC Motor. Problem with DiffEquation...

Hello,

I'm modeling the simple DC motor system in Maple.
The equations describing the system;

eq1:=J*diff(theta(t),t,t)+b*diff(theta(t),t)=K*i(t):
eq2:=L*diff(i(t),t)+R*i(t)=V(t)-K*diff(theta(t),t):
DCMotor:=[eq1,eq2];

First, I create the system using DiffEquation:

Sys:=DiffEquation(DCMotor,[V(t)],[theta(t)]);

And now I have problem. The input var is V(t) (input voltage) and the output var is theta(t) (position of the rotor).

But I wont to have in output var not position of the rotor but speed of the fotor - diff(theta(t),t)

How to set output var for diff(theta(t),t) (the speed of the motor)?

Best

Rariusz

Problem with procedure...

I tried to make a procedure that would find the determinant of any 3x3 matrix but I keep getting unterminated procedure what should I change??

Code:
DeterminantMat:=proc(matA::Array)
local  i::integer,  j::integer,  x::integer,  y::integer,  A::integer,  B::integer,  indice::integer,  S::integer:
A:=0:  B:=0:  S:=0: i:= 1:  j:=1:  x:=1:  y := 1:  indice:=1:

for x from 1 to 3 do
indice := 1:
for i  from 1 to  3 do
if x = i then
end do:
else if x <>i then
for j from 1 to 3 do
if y = j then  end do:
else if indice = 1 then  A := matA(i,j):  indice := 2:
else if indice = 2 then  B:=matA(i,j):  indice :=3:
else if indice = 3 then  B :=B * matA(i,j):  indice :=4:
else if indice = 4 then  A := A * matA(i,j):
if x = 2 then  S := S + (-1 * matA(x,y)) * (A-B):
else  S := S + matA(x,y) * (A-B):   end if:   end if:
end do:
end if:
end do:
end do:
end proc:

Trace in Physics package at Maple14...

I try to repeat lines (25)-(28) at

http://www.maplesoft.com/support/help/maple/view.aspx?path=Physics%2fTrace#commentform

I use Maple 14. However, instead of (28) I get the following result:

It means that Maple 14 does not perceive p_\mu, k_\nu and m as scalar quantities. I would like to ask how to define these variables correctly.

Eigenvalues and Characteristic Polynomial...

Hi!

I am comptuing the eigenvalues and the characteristic polynomial of a 8 by 8 symmetric matrix, say M. Thus, we define the matrix M, and compute its charast. plynm. by

and its eigenvalues with the command

Well, Maple returns the charast. polynm. an dthe eigenvalues. But, if we compute p(E[k]), for k=1,...,8, thats is, the values of the polynomial p(x) in the eingenvalues, Maple not turns cero!!! I'm really confused ... anyone know what could be happening?

Maple attached file with this example. Thank very much for your help!!