r-x^4
u-(x^3)*y
v-x*y^3
w-y^4

Step 1
with(Groebner):
K := {r-x^4,u-(x^3)*y,v-x*y^3,w-y^4};
G := Basis(K, 'tord', degrevlex(r,u,v,w));

Step 2
G intersect k[r,u,v,w]

G2 := Basis({w*G[1],w*G[2],w*G[3],w*G[4],w*G[5],w*G[6],w*G[7],w*G[8],w*G[9],w*G[10],w*G[11],w*G[12],w*G[13],w*G[14], (1-w)*K[1], (1-w)*K[2], (1-w)*K[3], (1-w)*K[4]}, 'tord', degrevlex(r,u,v,w));

G2 has too much though include the correct solution

correct solution...

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

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 there i have function as:

f(B,F)=

where: F,B>0 and 

,

such that: 

h=B^2;

I need to plot the surface (contour plot) f(B,F) at 0 and -ve points i.e f(B,F)<0.

Thanks

Hi all I have the following ODE

ODE_1:=diff(w(r),r)+R*(diff(w(r), r))^3=K/r;

ODE_T_1:=collect(algsubs(w(r)=u(r)+(r-1)*(U-0)/(delta-1), ODE_1),diff(u(r), r)) ;

 and when I try to do the following integration,

eq1:=int(phi[i](r)*ODE_T_1,r=1..delta) assuming delta > 1;

it gives the following error

"

Error, (in assuming) when calling 'int'. Received: 'wrong number (or type) of arguments:

i am doing a Q Laplace and finally summation it

however, i find an AccurateQSummation which is not like Summation

Should i use it, if not, use SumTools' summation how to simplify it?

restart;
with(QDifferenceEquations):
qexp := QPochhammer(-(-p*t), q, infinity);
f:= qexp*t;
b := subs(t=q^n,f);
QSimplify((1-q)*AccurateQSummation(b*q^n,n=-infinity..infinity));

with(SumTools):
QSimplify((1-q)*Summation(b*q^n,n=-infinity..infinity));

Is it possible for someone to provide the roller coaster maplesim model?

Hi

i know to make an asymptotic expansion one uses asympt(1+x,x); how do you make it go to -infity rather than plus infinity?

Hi,

I'm trying to modify the size of the window where an animation is shown. Tipically, this window has axis going from xmin to xmax and ymin to ymax but this causes some scaling problems when xmax-xmin>>ymax-ymin. I would like to focus in a certain region of my animation and change this viewpoint depending on the position of the elements of the animation. Can anybody help me?

This is my code:

> tf := 5;
> nfr := 201;

I have a huge PDE system included reaction-diffusion system with 4 independent and 12 independent variables. I read on the maple help for pdsolve/numeric/method. I found it is only for single partial differential equations in two independent variables. Is there anyone knows the ability of Maple for such a system?

1.

qRacah := qhyperterm([q^(-x), alpha*q^(-x)/(gamma*delta)], [alpha*q], q, gamma*delta*q^(x+1)*t, k)*qhyperterm([beta*delta*q^(x+1), gamma*q^(x+1)], [beta*q], q, q^(-x)*t, k);

i use inverse laplace to transform ordinary generating function into exponential generating function

can not further evaluated

2.

i do not know whether x is main variable or q is main variable or both are.

ff := eval(subs(x=0, exponential_genfun));

Hi,

I am aplying a simple energie balance with some simplification, and i am implementing it on maplesim.

The component has the following equations:

Am i missing something? Do i need any aditional equation to maplesim know to where should the fluid be moving when Qjacto is positive or negative? Why do i have...

Hi,

I am trying to perform a simple model of the entrainment of air in a jet by aplying energy emand momentum balances to the injector and flute. But i am having doubts on how to prescribe pressure on a point. This is the equations of the component Injector:

In the component porti...

I want to solve 4 nonlinear equation using maple. I want to find many roots of this 4 nonlinear equation.

I have no idea about initial condition so I use 4 for initial condition. my program is as follow:

f1:=nonlinear equation1(include variable q1,q2,p1,p2);

f2:=nonlinear equation2(include variable q1,q2,p1,p2);

f3:=nonlinear equation3(include variable q1,q2,p1,p2);

f4:=nonlinear equation4(include variable q1,q2,p1,p2);