g1 := x^2-w*y;
g2 := x*y-w*z;
g3 := y^2-x*z;
S12 := SPolynomial(g1, g2, plex(x, y, z, w));
eq1:= S12 = u1*g1 + u2*g2 + u3*g3;
with(LinearAlgebra):
T:=GenerateMatrix([eq1],[u1,u2, u3],augmented);
S:=LinearSolve(T,free=t);

(u1,u2,u3) should be (0,0,-w);

the goal is find a u1,u2,u3 to make eq1 equal S12, where u1,u2,u3 should be in terms of x y z w

g3 := y^2-x*z;
X3 := LeadingPowerOfProduct(g3, plex(x,y,z,w));

X3 should be y^2

I noticed in math mode the alt windows shortcut codes won't work but in text mode it will, understandably.

ie.  alt-0125 }   alt-0254 þ   etc...

Is there a way to get maple to output the list of alt codes?

I am doing an algorithm that include a " for " loop, inside that loop in the first statement I assigned a certain matrix to a variable G and then the algorthim changes some elements in G but when the algorithm repeat the second step of the loop it uses the last resulted matrix of the previous step NOT the one that is defined in the first statement of the loop. I want it to start each step of the loop with the same value of the matrix. what should I do?

3*rho1 - 2*rho2 + rho3 - rho4 = -1

4*rho1 +   rho2 - rho3        = 5

original without cost function:

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,(y1^1000)*(y2^1)*(y3^1)*(y4^100)- x1*x2*w + 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));

after have cost function 1000*rho1 + rho2...

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,x1*x2*w - 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));
(1,1,0,2)

subs({rho1=1,rho2=1,rho3=0,rho4=2},3*rho1 - 2*rho2 + rho3 - rho4);
subs({rho1=1,rho2=1,rho3=0,rho4=2},4*rho1 +   rho2 - rho3       );

can also reduce with G8,G7

R1 := Reduce((x2^(5+1))*(w^1), {G[7...

I was trying to solve for the equilibrium points for a system of differential equations.  Solving for the last equation, I encountered an error that I can't seem to find anywhere else.  "Error, (in Engine:-Dispatch) not implemented yet: 13". The code is below if anyone wants to take a go a figuring what's up.  Thanks!

A:=Asource/muA;
                            Asource
                            -------
                              muA  

(I'm french, scuse me about my bad english)

Hi everyone, 

I'm new with maple and I have a big question. I tried really hard for 2 days to find the answer but... 

My problem it's really simple, i think...

this is the end of the program where is my problem:

      dsn1 := dsolve(dfin, numeric); 

      dsn1(1);

      proc(x_rkf45)  ...  end;

Here is my sheet. non_linear_P2_last_q.mw

In this problem I have achieved every thing what I aimed, thnaks to the wonderful people here!!!

I am looking for a way to do it in a more structured way meaning the calculation of Jacobian in the attached sheet is somewhat manuel, though I wonder does maple has some other way?

Hi. I am creating a costum step and i did this

I am attempting to write a procedure for the numerical integration methods of the composite trapezoidal rule and Simpson's 1/3 rule, however the results that are being produced when I try to test the procedure are differing to the computed integral using the int command.

Here is my attempt so far:

> f:= x -> cos(x)*exp(-x/4):

> CompTrap:= proc(f,a,b,n)
local i, aa, h;
h:= (b-a)/n;
aa:= f(a)+f(b);

for i from 1 to n by 1 do

int(mul(sin(x/2^k)*2^k/x, k = 0 .. 123), x = 0 .. infinity) in closed form  with Maple?
I conjecture that the answer is Pi/2 and 123 can be replaced by any natural (i.e. ::posint) number.

Hello

This is actually somewhat the same question I asked some months ago. I found a way to go around it, but now unfortunately I'm back to the same problem, and I didn't have any answers when I asked.

Here it is a system that I solve numerically:

##############################################

eq1 := diff(W(r), r) = -(1-beta*(W(r)-W0))*(M(r)+4*Pi*r^3*p(r))/(beta*r*(r-2*M(r)));

eq2 := diff(M(r), r) = 4*Pi*r^2*rho(r);

W0 := solve(thetaR = theta0-W0, W0);

I just found out something which I don't understand.

Why am I getting this answer?

solve(a*x = 0, {x});
                                   {x = 0}

When I try some similar expressions I see it behaves differently for equality and inequality.

solve(a*x <= 0, {x});...

Hello everybody,

I'm completely new to Maple (installed it today) and have a question about the use of multi variabele expressions. I wanted to use a function B depending from the variables a, b, c and d without describing the expression explicit. 

For example: B:= (a, b, c, d);            [without the -> a....]

Then I wanted to take the taylor expansion of the first order for [a=e, b=f, c=g, d=h...