Hi! I have a problem that is driving me crazy. It is related to shooting methods for ODEs with Maple. Essentially what one needs is to find the solution to two equations simultaneously, defined in terms of the numerical solution to a second-order differential equation. I have used fsolve and also DirectSearch but nothing works. Graphically one can see the solution to the problem exists.

I include the Maple file. Thank you. 

Download shooting_question.mw

How can we collect coefficients of different powers of involved variables form Eq, to get system of algebraic equations, like x,y,z,t but these are also appear in JacobiCN function? Maple file link is attached below, please help me to solve this issue. 

collection_of_coefficients.mw

the function is true but i want to be sure when i use pdetest must give me zero, but there must be a way for checking such function, please if your pc not strong don't click the command pdetest, i want use explore for such function but i am not sure it work or not, becuase the graph are a little bit strange  and long , i want  a way for easy plotting and visualization of such graph , can anyone help for solve this issue?

 sol.mw

Hi,

I want to solve a first order ODE and (i) plot phi vs. x and then (ii) export data of the plotted curve as an T-shape ASCII file. I try to do this by Maple, but there are some errors and I couldn’t get the mentioned curve and ASCII file. Please, help me to remove errors:
(dphi/dx)**2+2*V(phi)=0,

where V(phi)= (1-alpha)*M^2 - (1-alpha)*M*sqrt(M^2 - 2*phi) + mu*(mu + beta*nu)*(1 - exp(phi/(mu + beta*nu))) + 
            (nu/beta)*(mu + beta*nu)*(1 - exp(beta*phi/(mu + beta*nu))) + (alpha/gamma1)*(1 - exp(-gamma1*phi)),

where we assume mu := 0.01, nu := 1 - mu, beta := 0.05, alpha := 0.3, M := sqrt(0.704), gamma1 := 0.001.
As it is seen from the attached figure, V(phi) and dV(phi)/dx =0 for phi=0 and phi=phi_m (two extreme points of V(phi)) and d^V(phi)/dphi^2<0 at phi=0 and phi=phi_m (phi and x are real). 

From the attached figure, it seems that phi=phi_m at x=0 and phi->0 as x-> infinity. (Hint: I expect the plot of V(phi) vs. phi and phi vs. x be similar to the curve "A" in the attached file).

Thanks.
rk4.mw

i did pdetest without conjugate like the paper did i get zero but when i did pde test  with conjugate i didn't where is my problem 
i will do without conjugate but how change p[2]=conjugate(p[1])

Download p1.mw

I am in the process of trying to learn Maple. As an exercise I am trying to plot an amplitude modulated waveform. No problem plotting the waveform of the carrier and the modulating signal but for some reason (user error on my part) I can't plot the combined waveform y given by. 

y = [1 + M*cos(2*Pi*f__m*t)]*A*sin(2*Pi*f__c*t)

Where the following values have been set M = 0.5, A = 1, f_m = 2, f_c = 10.

Plotting y using the following

plot(y, t = 0 .. 1, title = "Graph of Carrier", labels = ["time (seconds)", "Amplitude (Volts)"], color = red);

Gives the following error. Having tried various combinations of variables, values, single parameter functions etc I need a pointer as to where I am going wrong. 

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct
What I am expecting is a graph similar to the below derived in Excel? My version of Maple is 2023.2

VC1.mw

How to solve this error?

i don't know how apply conversation language to matlab in righ hand side  don't show up to do conversation language for short is come up but for this not 

Download convert-to-matlab.mw

I have plotted a 3D figure by maple 2023. but all numbers and units on axes are seen a black boxes. How can fix this problem?

Dear Maple Community,

I come to you with a question about the reduced involutive form (rif) package. Namely, I decided to try the classic example from the "LONG GUIDE TO THE STANDARD FORM PACKAGE", which dates back to 1993. Here is the link to the complete documentation:

https://wayback.cecm.sfu.ca/~wittkopf/files/standard_manual.txt

So, the example is the following:

2.1 SIMPLE EXAMPLES

EXAMPLE A

Consider the system of nonlinear PDEs:       

y Zxxx - x Zxyy  =  Zyy - y Zy

                        2     2    2
2 y x Zxxx Zxyy + x Zxxx + x y Zxyy  =  0

                  2    2
y Zxyy - x W + 2 x  y Z  =  0

                 2    2
Zyy - y Zy  + 2 x  y W  =  x W

where the dependent variables W and Z are functions of the
independent variables x and y, and Zxxx denotes the third partial
derivative of Z with respect to x etc.

After making computations back in 1993 with Maple V, they obtain the following involutive form:

In our original notation the (considerably) simplified system is:
                                            2
  Zxxx = 0, Zxy = 0, Zyy = y Zy, W = 2 x y Z

So, I tried the same system of PDEs in the modern Maple and the modern rifsimp() command. You can find the result below:

I demo_question.mw

The problem is that nowadays [Maple 2022.1] , I get only the trivial solution: z = 0 and w = 0.

Could someone clarify, please, where the truth is and what am I doing wrong?

Thanks a lot in advance for any help and clarification!

Best regards,

Dr. Denys D.
 

I have one random variable as follows:  r := sqrt(RandomVariable(Uniform(0,1)))

I have another as follows:  x := cos(RandomVariable(Uniform(0, 2*Pi)));

The r variable returns a PDF of 2*t, where t<1; 0 otherwise.  The variable x returns a PDF of 1/(Pi*sqrt(1-t^2)), where |t|<1; 0 otherwise.

For both random variables, I have validated through simulated values.  I am convinced the above PDF functions are correct.  When I multiply r*x in Maple and then ask for a PDF of the resultant random variable, I fail.  Maple seems to not like this.

My question is then as follows...can anyone help me get the PDF of r*x?

Which version of Maple contains its most optimal and elegant coding, regardless of utility,  and was it the product of someone's direct creative input (not auto-coded with AI for example).

As an aside, there's a spelling mistake in the word "separate" below within the "Tags" instruction:

"Tags are words are used to describe and categorize your content. Combine multiple words with dashes(-), and seperate tags with spaces."

I am making changes in my code in order to make it work with Grid. This involves making many modules I had export on them as local.

I have lots of code of this form

use toX=A:-some_very_long_module_name:-toX in
local B := module() 
       export foo:=proc()
            toX(...);
       end proc;
end module;
end use;

This is done so I do not have to type A:-some_very_long_module_name:-toX() all the time as toX() is a very common call I make all over the place in many different modules.

This was working fine, except when I now changed some_very_long_module_name module to be local.

I really do not want to change all my code and change toX(...) to the explicit fully qualified name  some_very_long_module_name:-toX(....)

I could probably look into using alias instead, but I do not like aliases.

But before doing this, any one knows why this now makes Maple not happy? And if it possible to make use....end use while keep this long named module local?

Below is work sheet showing 4 examples. The first one is how I had things before, where the module was export.

The second and third examples showing the problems that show up when changing the module to local.

The last one showing it works if I remove use...end use and just type the full long name.

my goal is to keep this long named module local, but still use  use...end use. around other modules which needs to make calls to it. Just to safe typing, that is all. 

Any ideas to try are welcome. 

Download trouble_with_use.mw

The following page of the documentation cites an article from the newsletter which I would be very interested in reading: 

- https://www.maplesoft.com/support/help/view.aspx?path=assume

The full citation for the article I'm interested in is as follows:
 

Download blogpost-algebra.mw