How to implement the following operation in Maple, as in MATLAB:

A(~isprime(A)) = 0

A =      0     3     2    13 
           5     0    11     0     
           0     0    ...

I'm trying to construct some tickmarks. I have so far found a way to do a "hat", where I use:

conjugate(alpha)

I also need to put a subscript and a superscript on that. Any ideas?

Cheers,

Christian

finding involutivedivision, autoreduction, involutivebasis for polynomial system,

which library contains these functions

I've a 3d surface built by the point's set;

It looks like this:

The question is:

1. how to build a surface based on this data?

2. how to define whether the arbitrary point (x,y,z) is inside this surface or not?

If I have a MapleSim file that successfully runs and I want to manipulate the data and equations in Maple, how would I do that? I know how to link a MapleSim file to a Maple document, and how to pull out the system of reduced DAE's, ODE's, parameters, etc. My question is whether or not there's a simple way to have Maple solve those equations and be able to display the same type of data (graphs, export the numerical solution to the system to an Excel file, etc.) that...

Hi,

I need to get inverse kinematics model of a double inverted pendulum. I got a model but it gives error when I run. It seems to me that the inverse kinematics model is not quite right. If you are able to help me, I will be very pleased.DoublePendulumI.msim 

Thanks for your help....

I would appreciate a hint why Maple 16 Classic under Win XP 32 bit crashes due to the following sequence of commands

> dir1:="C:\\Maple16\\P16\\test1\\";
> path:=cat("dir ",dir1);
> ssystem(path);

with the lasst ssystem command doing the crashing while the same sequence of commands under Maple 13 Classic does not: Instead, system returns to Maple 13 worksheet the content of the directory "dir1" as expected.

Is there another way to issue that ssystem command?

Hey,

So I have an array of plots [2,1] and I am trying to seperate them into two seperate plots to extract data using the getplot function (which only works on single plots).

I have looked all over the help guide but never found anything relating to my problem. Is there a way to seperate these two plots?

Thank you for the help.

Hello,

I have a system of two ODEs.  The solutions of the system (U(t) and V(t)) should never take on negative values, but may approach zero.  When using dsolve (numeric) to solve the system and odeplot to plot the system, I can see that dsolve gives slightly negative values to the solutions when close to zero.  This may result in the solutions becomes increasingly negative over time.

How can I prevent dsolve from assigning negative values to...

Hi everyone 

I want to solvethe following  differential equation system:

> equa1 := c[pa]*L*(diff(T[air](t), t)) = h*(T[sol](t)-T[air](t))+L*L[vap]*(diff(m[v0]-709.2750000*M*exp(1.000000000*10^(-8)*(1.311709131*10^9*T[air](t)-4.892675059*10^11)/T[air](t))/(R*T[air](t)), t));
> equa2 := c[ps]*l*(diff(T[sol](t), t)) = h*(T[air](t)-T[sol](t))-exp(-opt*L*(m[v0]-709.2750000*M*exp(1.000000000*10^(-8)*(1.311709131*10^9*T[air](t)-4.892675059*10^11)/T[air...

Hi guys,

i solve the equations of motion for a triple pendulum(Eq17,Eq27,Eq28). After that, I applied numerical values and initial conditions to solve these equations numerically. But unfortunatially i get the following error.

> g := 9, 81; l[1] := 5; l[2] := 3; l[3] := 1; m[1] := 5; m[2] := 3; m[3] := 1;
> ini := `ϕ`[1](0) = (1/3)*Pi, (D(`ϕ`[1]))(0) = 0, `ϕ`[2](0) = (2/3)*Pi, (D(`ϕ`[2]))(0) = 0, `ϕ`[3...

hello,

after defining N_1 N_2 and N_3 with long numbers, and defining M by N_1*N_2*N_3 i tried to step backwards.

By entering something like: M/N_1 (what is exactly like N_2*N_3) it does not show me full numbers, but instead 1/something*M.

After refreshing M and afterwards refreshing M/N_1 it totally worked out and showed me the product of N_2*N_3

Now my question is: how can i prevent this odd display format without refreshing everytime? (i tried...

I want to solve the inequality log[3](5*x+1)/(log[3](7*x-1)^2) <=0. I did

A:=log[3](5*x+1)/(log[3](7*x-1)^2):

solve(A<=0,x);

I got solution is x = 0. But this is a wrong anserw. Please help me. 

Another way, if i solve  the inequality log[3](5*x+1)/(log[3](7*x-1)^2) <=1, i got Warning, solutions may have been lost. Please help me.

hello all,

I work on an equation that represents the profit on a given horizon (L), according to three parameters (P, W, theta), and i want to find the triplet that minimizes this function. I even try to make a loop on the three parameter to find the minimum, but Maple does not give me a numerical result. Below is the function that I want to minimize (Prof(P,W,theta)):

restart;

beta:=1.2;

L:=2;

Q := t -> 50000/(2+23*exp(-120000*(W+1/10)^(1/10)*t/P^2)); ...