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

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?

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});...

I want to get a closed form of the nth order differentiation: diff(exp(-(x+sigma)^2/(2*sigma^2)), `$`(x, n)) where sigma is a positive constate. But maple gives an unevaluated answer. Howere, if calculating diff(exp(-x^2/sigma^2), `$`(x, n)),maple gives a closed form:x^(-n)*2^n*MeijerG([[0, 1/2], []], [[0], [1/2+(1/2)*n, (1/2)*n]], x^2/sigma^2).

Hey,

I have some expression of the form

sqrt((a+b)*(a-b))/(sqrt(a^2-b^2))

which maple unfortuntely does not simplify...

I have tried assumptions as a>0 b

how can one make this work?

Hi! Today I tried Maple for the first time and unfortunately encountered a silly but really annoying problem. When I write a simple calculation expression like 2+2; it works as expected. However, when I use an expression like sin, ln, Pi, etc., the return (enter) key stops working, the program just doesn't accept them. I'm beginning to think this has something to do with my (quite old) Windows setup since I couldn't find a similar error anywhere in the internet.

I'm working on;

The command

f := -.1*(x+3)*(x-1)*(x-3.5):
Student[Calculus1][FunctionChart](f,x=-4..4,concavity=[],title="",caption="",gridlines);

produces different results in Maple 13 and Maple 16.02.

In particular, the concavity=[] option seems to be ignored in 16.

Also, it may be my imagination, but execution seems slower in 16, but I didn't try timing it.

My apologies, if this has already been noted.

Hi

Could someone pls help write me code such that when text is entered, the first chemical element symbol in each word is identified and highlighed.

i/p: "the cat molly"

o/p: "tHe Cat Molly"

( like the opening title sequence on "Breaking Bad")

I got out the sledghammer and made a start

element.mw

Hello,

I am trying to get a shaded area in my plot but I could not.

First we solve ODE without randomness:

ode := diff(U(t), t) = -(A+B*U(t))*U(t);

Then we add randomness to ODE and solve:

ode2 := diff(U(t), t) = -(A+r(t)+B*U(t))*U(t);

A with randomness for r in R=( - 0.0001/365, 0.0001/365) is:

A(t,r)= A+r

Where A is constant =  0.0001/365

We plot both solution. For the plot I...

relativistic_7.mwrelativistic_7.mwrelativistic_7.mw

I want to evaluate the function A(xi(t),eta(t),diff(xi(t),t),diff(eta(t),t)) atthe values of xi(t), eta(t), diff(xi(t),t),diff(eta(t),t)  given above.

Thank you