Dear all

I would like to find the optimal parameters of the attached optimal control problem 

The system of ODEs is defined and boundary conditions and The hamilton Jacobi Bellman equations is written.

Optimal_control.mw

Many thinks

I'm currently working my thesis and can someone help me to write a code to solve this IVP

u_t + 2×u²×u_x−(u_x)²−½×u_{xx}×u=0

with initial condition

u(x,0)=-tanh(x)

Dear friends,

Greetings.

How to get the second solution.

how to change the guess value in maple.

figure 1 plot in Matlab with two different initial guesses.

TWOSOLUTION.mw

I am writing a maths books using maple now. It is fantastic to use maple for writing books in maths.
 

Download cannotCopyWhy.mw

I enclose a part of my document where in I made a particular line with text and maths formats combined.Then I made changes in the line. Now copy paste does work only for the later half (both text and maths formats). The corrected first part is not being copied.

How do I do the corrections properly so that copy paste is not a problem at laer stages.

Thanks for the answer.

Ramakrishnan V

The equation sin(9*x-(1/3)*Pi) = sin(7*x-(1/3)*Pi) can be solved easy by hand with solutions k*Pi and -Pi/48 + K*Pi/8. With Maple, I tried 
solve({sin(9*x-(1/3)*Pi) = sin(7*x-(1/3)*Pi)}, x, explicit, allsolutions)

I don't get the above solutions. How can I get these solutions?

I use patmatch to look for certain expression inside a larger expression.

I find sometimes I need to repeat the same code to check for  "... + ..."   and also ".... * .....", since I do not know to tell Maple to look for + or * in the same code. *Luckily, I do not have to check for "-" or "/" operators, since "+" match with "-" and "*" match with "/").

An example will make things more clear.

Suppose I want to see if sin(x)*sqrt(x*y) has sqrt(x*y) anywhere in it. So I first try

restart;
expr:= sin(x)*sqrt(x*y);
if patmatch(expr,a::anything+(b::anything*x*y)^(c::anything),'la') then
    assign(la);
    if c =1/2 or c=-1/2 then
       print("found sqrt(x*y)");
    else
       print("did not find sqrt(x*y)");
    fi;
 else
   print("did not find sqrt(x*y)");
 fi;

And this fails, since I used "+" inside the patmatch. Then I try '*" instead

if patmatch(expr,a::anything*(b::anything*x*y)^(c::anything),'la') then

And now it does match.

What I'd like to write, is something like this (which ofcourse does not work)

if patmatch(expr,a::anything (* or +) (b::anything*x*y)^(c::anything),'la') then

I looked at conditional in patmatch, but it does not seem to apply for the above.

Any suggestions?

Maple 2019.1 on windows 10

Dear Users!

Hope everyone should fine here. I want to know how to multiple the entries of a matrix with a matrix like If A is matrix as:

A:= Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 2})
I want nonzero entries should multiply with m*m identity matrix and zero entries multiply with null matrix of order m*m. for m=2 the desired results i calculated manually as:

Matrix(4, 4, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 2, (3, 4) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 2});

I want to know the general way. Thanks in advance

most people who post here seem to use .mw written in 2D, which I do not like to use.

Is there a tool to convert such a file to 1D worksheet that one can use from the command line before opening the document itself in Maple?

The reason I ask, sometimes opening the original file and trying to do this from inside Maple by selecting the code using the mouse, then  Format->ConvertTo->1D   does not work, and gives an error.

Also, sometimes, when I try to first create an empty worksheet document, and then try to copy/paste the code from the other document over, it also does not work. This happens when there are syntax errors in the original document. The error that comes up is

        Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

As an example, please see the attached file in the following question

https://www.mapleprimes.com/questions/227506--I-Need-Help-Trying-To-Write-A-Code

It will good to have a tool that converts such documents to 1D worksheet or even plain Maple code (.mpl) but I did not see such option under SAVE AS either. Also, when I tried to export it as .mpl file, I get the same error as above in the file. So I gave up.

Hello Anybody can help me to write codes for PDE to solve by Galerkin finite element method or any other methods can be able to gain results? parameter omega is unknown and should be determined.

I attached a pdf file for more .

Thanks so much

fem2
 

Download fem2

buchanan2005.pdf

Download Dtm.mw

I found this strange result:
When the central moment of order 2 and the variance are assessed on a sample of a random variable, they don't return the same value.

Let S be the sample, N its size and M its empirical mean.  The difference comes from the fact that 

• Variance(S)               = add( (S[n]-M)^2, n=1..N) / (N-1)
• CentralMoment(S, 2) = add( (S[n]-M)^2, n=1..N) / N

By definition the variance of a random variable X is its 2nd order central moment.
This definition should also apply when these statistics are calculated on a sample of X.

It seems to me this is a mistake.
 

Download StrangerThings.mw

Hello. Let's say I have a set A1 in which there are nested sets {5,6},{8,9}. How can I get A2 from A1? Expand nested sets so that all elements are in the same set

The eval command doesn't work on random variables.

For instance:

with(Statistics):
X := RandomVariable(Normal(mu, sigma)):
Y := eval(X, mu=1);
Mean(Y); # returns mu

I once found a Maple command whose syntax is the same as eval's and acts as eval
Y := COMMAND(X, mu=1)
Mean(Y); # returns 1

Could someone remind me of its name?

Thanks in advance

How can i plot a probability function such as cos(x-y)*cos(y-z)*cos^3(x-2z)=0.6 where

x=0..5, y=0..x, z=0..y.

please guide me.

Dear Maple users,

I am progressing, but one last hitch, see below. A want the invariants of the PDE below. However, the final expression is too general to be useful. I would like to insert specific values for F1(R), F2(R), F3(R), and F4(R):

F3=F4=0; and F2=1, and F3(R)=aR, where a is a constant. R is one of my independent variables (the Reynolds number). 

I would like to do this at the step where the Infinitiesimals are generated by

For example, the first entry would then be, infinies:=[ _xi[y](y, R, l, u) = y, ... ]. Then Invariants should give me much simplified expressions which I need.  How can i do this?

Thanks

Nadeem