Hello, I am trying to plot this one and I did not get the graph, I don't know where the problem. Could you please help. thank you

restart;
with(plots):

alpha := 0.005;
beta := 0.01;
                          beta := 0.01

mu := 0.005;
                          mu := 0.005

m := 0.003;
                           m := 0.003

u := 0.3;
                            u := 0.3

w := 0.3;
                            w := 0.3

v := 1;
                             v := 1

gamma = 0.1;
                          gamma = 0.1

Gamma := 0.01;
                         Gamma := 0.01

Eq := sigma^4 + sigma^3*[2*(mu*v + mu) + 2*alpha - beta*(u + w) + 2*beta*v] + sigma^2*[4*mu*m + 2*alpha*(2*mu + m) + alpha^2 - 2*beta*v(mu - 2*m) - 2*beta*(u + w)*(mu + m) + alpha*beta*(4*beta + (-u - w)) - beta^2*v*(u + w - v) + k^2*(Gamma^2 + gamma^2)] + sigma*[2*alpha*mu*(2*m - beta*(u + w)) + 2*alpha^2*(beta*v + mu) - 4*beta*mu*m*(u + w + v) + 4*alpha*beta*v*(mu + m) - 2*beta^2*v*(u + w) + 2*beta^2*v^2*(m + alpha) - alpha*beta^2*v*(u + w) + 2*Gamma^2*beta*v*k^2 + 2*k^2*(Gamma^2*mu + gamma^2*m) + gamma^2*k^2*(2*alpha - beta*(u + w))] + [Gamma^2*k^2*(2*beta*mu*v + gamma^2) + Gamma^2*beta^2*k^2*v*(Gamma*v - gamma*(u + w)) - beta*gamma^2*k^2*(u + w)*(alpha + 2*m) + alpha^2*(beta^2*v^2 + gamma^2*k^2) + 2*alpha*beta^2*v*(v*m - mu*(u + w)) + 4*beta*mu*m*v*(alpha - beta*(u + w))];
p1 := implicitplot(Eq, k = 0 .. 10, sigma = -0.1 .. 0.1);
           4
Eq := sigma 

          3                                                       
   + sigma  [2 mu v + 2 alpha + 2 mu - beta (u + w) + 2 beta v] + 

       2 [                                   2
  sigma  [4 mu m + 2 alpha (2 mu + m) + alpha 

   - 2 beta v(mu - 2 m) - 2 beta (u + w) (mu + m)

                                       2              
   + alpha beta (4 beta - u - w) - beta  v (u + w - v)

      2 /     2        2\]         [               
   + k  \Gamma  + gamma /] + sigma [2 alpha mu (2 m

                            2              
   - beta (u + w)) + 2 alpha  (beta v + mu)

   - 4 beta mu m (u + w + v) + 4 alpha beta v (mu + m)

           2                   2  2            
   - 2 beta  v (u + w) + 2 beta  v  (m + alpha)

               2                    2         2
   - alpha beta  v (u + w) + 2 Gamma  beta v k 

        2 /     2           2  \
   + 2 k  \Gamma  mu + gamma  m/

          2  2                         ]   [     2  2 /           
   + gamma  k  (2 alpha - beta (u + w))] + [Gamma  k  \2 beta mu v

          2\        2     2  2                            
   + gamma / + Gamma  beta  k  v (Gamma v - gamma (u + w))

               2  2                      
   - beta gamma  k  (u + w) (alpha + 2 m)

          2 /    2  2        2  2\
   + alpha  \beta  v  + gamma  k /

                 2                     
   + 2 alpha beta  v (v m - mu (u + w))

                                         ]
   + 4 beta mu m v (alpha - beta (u + w))]

[k1,0], [k2,0], [k3,0], ...) , where k_{n}=2*n*Pi/L and the domain [0,L] 
L:=10.0;
p2:=pointplot([seq([2*n*Pi/L,0],n=0..10)],color=red,symbolsize=15):
display(p1,p2);
                           L := 10.0

       /            /     4
display\implicitplot\sigma 

          3                                                       
   + sigma  [2 mu v + 2 alpha + 2 mu - beta (u + w) + 2 beta v] + 

       2 [                                   2
  sigma  [4 mu m + 2 alpha (2 mu + m) + alpha 

   - 2 beta v(mu - 2 m) - 2 beta (u + w) (mu + m)

                                       2              
   + alpha beta (4 beta - u - w) - beta  v (u + w - v)

      2 /     2        2\]         [               
   + k  \Gamma  + gamma /] + sigma [2 alpha mu (2 m

                            2              
   - beta (u + w)) + 2 alpha  (beta v + mu)

   - 4 beta mu m (u + w + v) + 4 alpha beta v (mu + m)

           2                   2  2            
   - 2 beta  v (u + w) + 2 beta  v  (m + alpha)

               2                    2         2
   - alpha beta  v (u + w) + 2 Gamma  beta v k 

        2 /     2           2  \
   + 2 k  \Gamma  mu + gamma  m/

          2  2                         ]   [     2  2 /           
   + gamma  k  (2 alpha - beta (u + w))] + [Gamma  k  \2 beta mu v

          2\        2     2  2                            
   + gamma / + Gamma  beta  k  v (Gamma v - gamma (u + w))

               2  2                      
   - beta gamma  k  (u + w) (alpha + 2 m)

          2 /    2  2        2  2\
   + alpha  \beta  v  + gamma  k /

                 2                     
   + 2 alpha beta  v (v m - mu (u + w))

                                         ]               
   + 4 beta mu m v (alpha - beta (u + w))], k = 0 .. 10, 

                     \                                         
  sigma = -0.1 .. 0.1/, pointplot([[0., 0], [0.6283185308, 0], 

  [1.256637062, 0], [1.884955592, 0], [2.513274124, 0], 

  [3.141592654, 0], [3.769911184, 0], [4.398229716, 0], 

  [5.026548246, 0], [5.654866778, 0], [6.283185308, 0]], 

                               \
  color = red, symbolsize = 15)/

Hello.

I plot  F1:=1/(x + μ) + x - 1 and on the graph i get the line x = -μ. How can I get rid of it? I try the discont option but I get an error. Do I need, maybe, a specific with command?

The command I use is the next one:

diagd1_1 := implicitplot(F1 = 0, mu = -8 .. 8, x = -4 .. 4, gridrefine = 5, crossingrefine = 15, rational, color = blue, thickness = 2, numpoints = 100000, resolution = 3000, rangeasview)

The error I get:

Error, (in plot/options2d) unexpected option: discont = True

Thanks a lot !!!

I recently found this piece of code while using showstat:

if type(i,{'symbol', not ('symbol')}) then
...
end if

At first glance, this felt strange to me, because it looked like if we were checking whether "i" was of the "symbol" type or of a type other than "symbol", and therefore the test was still true.

Doing a simple test I obtain this

if type(s,{'symbol', not ('symbol')}) then "TRUE" else "FALSE" end if;
                             "TRUE"
if type(cos(c),{'symbol', not ('symbol')}) then "TRUE" else "FALSE" end if;
                            "FALSE"

Could you explain me what  not ('symbol') really means?

TIA

I have this equation: 
cpd = 2/(1 + exp(-0.004*200)) - 1

i want it to give me a new equation for cp (not solving the equation). 
what i've tried:

2/(1 + exp(-0.004*200)) - 1 = cpd = 0.379948962
q1 := 0.379948 = 2/(1 + exp(-0.004*cp)) - 1
isolate(q1, cp) = cp = 199.99 (this is the answer but i want the equation)

if this is not clear enough then i will give a simpe example: 

given equation: x-2 = 5
WHAT I DON'T WANT: x = 3

WHAT I WANT: x=5-2

For example, can a worksheet display a radius vector rotating through equally spaced points on a circle at a constantly reducing angular velocity?

i have a big code some times maple get stuck in a certain evaluation  
How can I issue a command to Maple to stop an command  if the time expands for a specific time?

Hello there, 

When I tried to use 'solve' command to solve an algebraic matrix equation, I got this error:

Would you please tell me the correct way to solve it?

eq6_3_sol := T__PhPh = solve(eq6_3_1, T__PhPh);
Error, invalid input: solve expects its 1st argument, eqs, to be of type {`and`, `not`, `or`, algebraic, relation(algebraic), ({list, set})({`and`, `not`, `or`, algebraic, relation(algebraic)})}, but received (Vector(3, {(1) = cos(omega__o*t), (2) = -cos(omega__o*t+(1/3)*Pi), (3) = -sin(omega__o*t+(1/6)*Pi)})) = T__PhPh . (Vector(3, {(1) = cos(omega__i*t), (2) = -cos(omega__i*t+(1/3)*Pi), (3) = -sin(omega__i*t+(1/6)*Pi)}))
 

The matrix equation I tried to solve is presented below:

(sorry for the duplication below, perhaps the Maple app which is supposed to display the content of the worksheet is not too happy with Microsoft Edge browser)

Download Q20210127.mw

hi, how to find the range in which root of Eq lies and why fsolve is not giving result>?

numeric_solution.mw

Hi, I hope this question is allowed

I wish to become well learned in maple but I find all the time, when I come across a new problem, the maple help function only takes me so far. In this community are some really helpful and bright characters which I would like to match in knowledge.

Is it possible to 'know it all' in the maple universe, or are there certain areas you're well versed in over others? 
 

what are some of the ways you have learned to use maple effectively?

kind regards, Ben

In using maple 17, for a for loop, the alias command is valid outside the for loop, but entering the for loop, maple verifies that alias is not recognized.The code is as follows:
with(PDEtools):
for k from 1 to 4 do
alias(u[k]=u[k](x,y,t));
alias(a[k]=a[k](x,y,t));
od:
print(has(a[1],x));
sol:=[]:
for i from 1 to 4 do
sol:=[sol[],a[i]=rhs(isolate(a[i]+diff(a[i],x)+i=0,diff(a[i],x)))];
print(has(a[i],x));
end do;
The output is as follows: thequestionofalias
true
[a[1] = a[1] + 1]
false
[a[1] = a[1] + 1, a[2] = a[2] + 2]
false
[a[1] = a[1] + 1, a[2] = a[2] + 2, a[3] = a[3] + 3]
false
[a[1] = a[1] + 1, a[2] = a[2] + 2, a[3] = a[3] + 3, a[4] = a[4] + 4]
false
Why is it outside the loop? A [1] ... A [4] has the variable x, but either a [ I ] or a [1] has no variable X.Send the question to the Forum and see which expert can answer it. See how you can make the alias command work inside a for loop! I'd appreciate it if you could answer that!

Is there a way to use the Button Component to stop a maple computation? I have two long functions running on nodes via grid in my app and I want to have them stopped via special button, is it possible to do in maple?

Hi everyone,

In some step of my code, I need to solve the following equation: M = Id, where, M is a square matrix of dimension n and each element of this matrix M is an expression with parameters. Id, is the identity matrix.

Could you please tell me how to solve such a matrix equation M = Id, whatever n ?

Example: i would like to solve M=Id, where M is given by:
M:=Matrix(3, 3, [[-1/3*a + b + c - 2*d, -1/4*a + 2/3*c, 1/15*a - 1/4*c + 2/3*d], [e + a - 2*c, 2/3*a, -1/4*a + 2/3*c], [f - 2*a, 0, 2/3*a]]);
Id:=IdentityMatrix(3);

Thanks in advance,
Best regards

Is there 

Using maple, approximate integral of 2/(x^2-4) from 0 to .30 by using

(a) the Trapezoidal rule: n=1, h=(b-a)/n, inegral of f(x) from a to b=(h/2)(f(x0)+f(x1))

(b) Simpson’s rule:

(c) the Midpoint rule: 

Hi,

I try to build a point, but a problem assigning the f function blocks me

Ideas? Thanks

ConiqueEllipseComplet.mw

As you can see in the code, I am trying to print the solution of the ODE using the Homotopy perturbation method for N = 10 ( f[0] till f[10]), but Maple is only printing f[0],f[1], and the rest of the terms are not printed out, what could cause this issue to take place?

``` 
 
restart;
N := 10:
  F(Y) :=  add(p^i*f[i](Y), i = 0..N);
  HPMEq := (1 - p)*diff(F(Y), Y $ 4) + p*(diff(F(Y), Y $ 4) + R*(diff(F(Y), Y $ 3)*F(Y)- diff(F(Y), Y $ 2)*diff(F(Y), Y $ 1))-G*diff(F(Y), Y $ 2));
    sol:=[]:
  for i from 0 to N do
      sol:= [ sol[],
              dsolve
              ( [ eval
                  ( coeff(HPMEq, p, i) = 0,
                    sol
                  ),
                  f[i](0) = -a,
                  f[i](1) = -b,
                  D(f[i])(0) = B * D(D(f[i]))(0),
                  D(f[i])(1) = -B * D(D(f[i]))(1)
                ]
              )
            ];
end do:
sol;

```
Please kindly try to run this code in your Maple version and tell me if all terms are printed.