Hi,
Apparently I have a problem but I can't find it. Please advise what is the source of the error?
Please see the attached worksheet.
1.mw

Download error_in_plot.mw

What is the mistake?

in this worksheet

restart;
alias(u = u(x, z, t), f = f(x, z, t));
                              u, f
u := (f+sqrt(R))*exp(I*R*x);
                    /     (1/2)\           
                    \f + R     / exp(I R x)
pde1 := I*(diff(u, z))+diff(u, x, x)+diff(u, t, t)+u*abs(u)*abs(u)-(u*abs(u)*abs(u))*abs(u)*abs(u);
    / d   \              / d  / d   \\           
  I |--- f| exp(I R x) + |--- |--- f|| exp(I R x)
    \ dz  /              \ dx \ dx  //           

           / d   \                /     (1/2)\  2           
     + 2 I |--- f| R exp(I R x) - \f + R     / R  exp(I R x)
           \ dx  /                                          

       / d  / d   \\           
     + |--- |--- f|| exp(I R x)
       \ dt \ dt  //           

                                                            2
       /     (1/2)\                           2 |     (1/2)| 
     + \f + R     / exp(I R x) (exp(-Im(R x)))  |f + R     | 

                                                            4
       /     (1/2)\                           4 |     (1/2)| 
     - \f + R     / exp(I R x) (exp(-Im(R x)))  |f + R     | 

simplify(%);
         / d   \              / d  / d   \\           
       I |--- f| exp(I R x) + |--- |--- f|| exp(I R x)
         \ dz  /              \ dx \ dx  //           

                / d   \                 2             
          + 2 I |--- f| R exp(I R x) - R  exp(I R x) f
                \ dx  /                               

             (5/2)              / d  / d   \\           
          - R      exp(I R x) + |--- |--- f|| exp(I R x)
                                \ dt \ dt  //           

                                               2  
                                   |     (1/2)|   
          + exp(I R x - 2 Im(R x)) |f + R     |  f

                                               2       
                                   |     (1/2)|   (1/2)
          + exp(I R x - 2 Im(R x)) |f + R     |  R     

                                               4  
                                   |     (1/2)|   
          - exp(I R x - 4 Im(R x)) |f + R     |  f

                                               4       
                                   |     (1/2)|   (1/2)
          - exp(I R x - 4 Im(R x)) |f + R     |  R     
collect(%, exp(I*R*x));
  /  (5/2)       / d   \      2       / d   \   / d  / d   \\
  |-R      + 2 I |--- f| R - R  f + I |--- f| + |--- |--- f||
  \              \ dx  /              \ dz  /   \ dx \ dx  //

       / d  / d   \\\           
     + |--- |--- f||| exp(I R x)
       \ dt \ dt  ///           

                                          2  
                              |     (1/2)|   
     + exp(I R x - 2 Im(R x)) |f + R     |  f

                                          2       
                              |     (1/2)|   (1/2)
     + exp(I R x - 2 Im(R x)) |f + R     |  R     

                                          4  
                              |     (1/2)|   
     - exp(I R x - 4 Im(R x)) |f + R     |  f

                                          4       
                              |     (1/2)|   (1/2)
     - exp(I R x - 4 Im(R x)) |f + R     |  R     
 

I see that using units with the maximize() function causes the connection to the kernel to be lost and then Maple (v2018) must be restarted for things to work properly.  This is obviously not desired behavior - is there any known workaround for this issue? (other than forgoing the use of units?).  Attached is a simple worksheet that illustrates this problem.  It has one part without units that works properly and one part with units that causes the error:  Units_Lose_Kernel.mw

Please help with the bifurcation diagram for the system and parameter values below

Download Bifurcation_Equation.mw

I am using temperature units with thermophysical data and scientific constants but am getting inconsistent behavior when using these.  In some cases (e.g., when calling thermophysical data), it seems that it's best to use a Temperature object.  However I have tried to use scientific constants in calculations and it seems that temperature expressions work best (i.e., when using degrees K rather than deg F or deg C).  Simplifying/combining units don't seem to work when using a Temperature object (or when using expressions with deg F).  It may be something simple I am overlooking but I can't figure out the pattern of behavior yet.  I've attached a file that demonstrates the issue:  Temperature_Object_Use.mw

Thanks for any insight here.

Hi,
I want to find (w/k)^2 from the following Eq. by Maple. How do I do it?
(u0b, mu,deltab,sigma and A are fixed parameters)

Eq.mw

am attaching the worksheet of the problem please help me to solve not able to  compute coupled

error.mw

Download AGM.mw
 

eq1:=( d)/(dt)u+(d^(2))/(dy^(2))u + s*( d)/(dy)u + delta * theta = 0;

eq2:=( d)/(dt)theta + (d^(2))/(dy^(2))theta + s*Pr*( d)/(dy)theta +lambda* exp(theta/(1 +(epsilon*theta))) = 0; 

initial and boundary conditons   

t <=0; u = theta = 0, for 0 <= y  <= 1   

t> 0;  u =0, theta = 0   at  y = 0;  

t> 0;  u =1, theta = 0   at   y = 1  ;

where, s, epsilon, Pr, lambda, delta are arbitrary parameters

Hello 

How to remove this error 

3d_plots.mw

Hello,

Can we impliment Artificial Neural Network for nonlinear coupled ODE equation with boundary conditions.? In maple

I wont seen any post regarding ANN in mapleprime.

Maple code for solving system of ODE using forward-backward sweep method.

When we specify a set (a sequence of objects enclosed in curly braces), Maple removes duplicates, since the elements of the set must be unique, that is, they cannot be repeated. See below for 2 examples. With the first example  {a<=b  and  b>=a}, everything is in order, since they are one and the same. But Maple treats the same equality, written in two ways  {a=b, b=a} , as different objects. It seems to me that this is not very convenient:

restart;
{a<=b, b>=a}; # OK
{a=b, b=a}; # not OK
is((a=b)=(b=a)); # not OK