Question: How I can solve this problem to plot

 

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)/

Please Wait...