I am a transplant from to Mathematica to Maple and still getting used to things. Could somebody please help me finding the maximum of f with repect to the constraints and with the variables below. I keep getting  syntax error, no matter how many times I check the documentation page. Thank you very much

f:=1/((4 a - 3 M) (-4 b + 3 M)) (-9 M^3 - 16 a b x +
    12 a M x + 12 M^2 x + 16 b x^2 - 12 M x^2 + 12 b M y +
    12 M^2 y - 16 b x y - 16 M x y - 12 M y^2 + 16 x y^2)...

hi dear

what is this error and how i can fix it?

Error, (in gcd/NoExpandDegreeCheck) inputs must be polynomials over all variables

i want to solve this pde

Is there a way to do a repitition statement with two indices?  Something along the lines of:

for i,j from 1 by 1 to N do a[i][j]:=Re(b[i][3][1][j]); end do; ?

I know that the syntax used here doesn't work, but is there a way to do this?

proc5.mw

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...

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...