Sorry if this has been already posted.

When print() is invoked from a proc into a module, non-English characters are not properly displayed with Maple 18.

It works ok if it is invoked from within the workbook.

Example:
print("Están en perspectiva")

Put this sentence in a proc into a module and the character "á" wont be displayed

Output: "Est�n en perspectiva"

Any hint about how to treat this issue?

Thank you very much.

César Lozada 

So im trying to write a maple script that computes the Jordan form of a given (3x3)- matrix
A. If {a,b,c} is a basis with respect to which A is in Jordan form, then I'm trying to make it
plot the three lines spanned by a, b and c, in the standard coordinate system. I was hinted to use plot3d here.

sidenote: I know how to compute the jordan matrix of A, such by find the eigen vectors and generalised eigen vectors and putting them in as columns in a 3x3 matrix say S,   where S is invertible    then  (S^-1)*(A)*(S) = (J).

Thanks in advance. <3

Hi:

i follow the code newton raphson for solve system nonlinear ODE in maple,every body have code for it?

When creating questions with Question Designer, how can I include notes that are visible to instructors only? I don't want students to see certain info. Is it in the "information fields"?   

I have a worksheet. I edit some variables' value. I then execute the whole worksheet to see how the graphs change. But the graphs generated by display command do not appear. All other calculations update fine. 

To overcome this. I save worksheet with my new variable values. Close the worksheet. Reopen it. Execute it. And now the graphs appear fine.

 I have a complex set of mechanical linkages and components, and I am trying to create a feedback control system for trajectory-tracking with the end-effector. They system is driven by a set of prismatic joints (representing ideal series-elastic actuators). At this stage I am simply trying to establish the control parameters and the required forces to achieve desired positions/velocities in the workspace, and have not yet inserted any realistic actuator characteristics.

An open-loop solution can be found using translational position control on the prismatic joints, this is no problem. However to close the loop properly, I need to convert the system to use a differential input, such as a velocity driver. At this point everything becomes very finicky and highly dependent on initial conditions. Generally Maplesim gives up, claiming that no solution can be found. Due to the complexity of the system, there is no way for me to guess a priori the appropriate initial conditions of all the joints and linkages. I was wondering if anybody had any tips for either pre-determining the initial parameters, or alternatively for relaxing the solver in some way to allow for a wider exploration of possible solutions. Or any other ideas!

So far: I have tried setting up an open loop position-driven system with the same trajectories and using Maple to read all the initial conditions, then transferring them as guesses to my closed-loop system. But this has not worked. As a temporary workaround I am using position drivers on the SEAs, which are in turn driven by a velocity controller (with an integrator term). However this is suboptimal. Any help would be much appreciated!

I  often find mself in the following quandry: I am investigating a problem, I do a few examples, run a do loop for j from1 to 5 to see if that works and then enter

for j from 1 to 200 do  a[j]:=blah, blah,...  end do:

and 15 minutes later nothing has happened. Is the loop finishing up with j = 198 or so or is it stuck at j = 23 and it will take forever to finish? Or perhaps there is an error that turns up for j > 23?

When I think about it I can insert a print(j) command to keep track of this but there are other situations where that doesn't work. For example,

Order := 35; mtaylor(....

When nothing happens is that because the calculation is almost done? or stuck? etc.

Is it possible to take a quick peek at the state of Maple and the decide whtether to continue the computation or abort? If the calculation of mtaylor(... is particularly long I might settle for the 15 coeficients computed so far rather than abort. Can I get at those?

Any advice?

Using the Fourier convolution theorem to solve f(t) =sin (t)

f(t)=R dJ(t)/dt+J(t)/C

Hi,

When I run two times my code, the results change.

Have you any idea, why???

Hi! 

When trying to find the fundamental solution of the Heat equation using Maple (software), I get the following Error message which seems to have no documentation available (?) :

Using :

          PDE := -(diff(f(x, t), t))+(diff(f(x, t), x, x))*Di = 0

assume(epsilon > 0);

pdsys := [PDE, f(x, 0) = Dirac(x-epsilon)];

pdsolve(pdsys, build)

"Error, (in casesplit/K) this version of casesplit is not yet handling the function: Dirac"

Anybody has an idea what that is? (Using Maple 17) . How can I solve this problem ? 

http://en.wikipedia.org/wiki/Heat_equation

Thanking you on Advance, 

Erez . 

Hi,

Please I need help in this subject. I would like to compare the numerical solution obtained by finite difference and pdsolve/numeric.

The equation considred is  diffusion Equation using Forward-time centered-space (FTCS) stencil
The code work well with Dirichlet boundary condition, but I want to let  x=-1  Dirichlet boundary condition but on x=1, we put a Neumann condition likeeval( diff(u(t,x),x),x=1)=1. Thank you very much to put the necessary in the attached code the changment.      
Many thinks.

Change_boundary_condition_in_procedure.mw    

assume f and g are unknown

and assume solve(f, x) = solve(g, x)

f -> a

g -> a

b -> f

b ->g

if assume f = (x+1)*(x+2), g = (x+2)*(x+3)

and a = (x+1)*(x+2)*(x+3)

would like to find map from (x+1)*(x+2) to (x+1)*(x+2)*(x+3)

is it the solution subs(x=(x+1)*(x+2),(x+1)*(x+2)*(x+3)) by composition?

subs(x=(x+1)*(x+2),(x+1)*(x+2)*(x+3))

subs(x=1, (x+1)*(x+2));
subs(x=2, (x+1)*(x+2));
subs(x=1, (x+1)*(x+2)*(x+3));
subs(x=2, (x+1)*(x+2)*(x+3));

6 -> 24
12 -> 60

subs(x=1, ((x+1)*(x+2)+1)*((x+1)*(x+2)+2)*((x+1)*(x+2)+3)); # not 24
subs(x=2, ((x+1)*(x+2)+1)*((x+1)*(x+2)+2)*((x+1)*(x+2)+3)); # not 60

it seems composition is wrong

more difficult and general case should be

f(x,t)  -> a(x,t)

g(x,t)  -> a(x,t)

b(x,t) -> f(x,t)

b(x,t) -> g(x,t)

solve(f(x,t), x) = solve(g(x,t), x) = in terms of t

Hi:

What methods to solve system nonlinear ordinary differential equations in maple?I follow a method to solve system nonlinear ode second order that very fast answer me?

Hi,

Please I need you to add in the output of my code the order of error defined in the procedure.

Thanks for helping me.

Here, the code.

QuestionNumber2.mw

Hi, I would like to solve this nonlinear problem :

with :

I use the NLPSolve command and i get this error : Error, (in Optimization:-NLPSolve) integer variables are not accepted.

Can you help me ?