Hello,

I need a bimodal distribution. Since I could not find any among the ones provided by Maple, I created a simple one:

with(Statistics):
U := Distribution(PDF = (proc (t) options operator, arrow; piecewise(t < -5, 0, t < 5, -(1/2000)*t^4+(9/1000)*t^2+7/80, 0) end proc)):
X := RandomVariable(U):

#Plotting PDF and CDF works fine:
plot(PDF(X, t), t = -infinity .. infinity);

plot(CDF(X, t), t = -infinity .. infinity)

However, plotting the quantile function does not work:

plot(Quantile(X, z), z = 0 .. 1);

it has a decreasing part for z<1/2 and a discontinuity at z=1/2.
I can plot it correctly as
plot('Quantile'(X, z), z = 0 .. 1);

but I wonder why the first option does not work for such a simple distribution.

I have a trouble of Maple.
I can't understand the solution of Maple for simplification of Dirac and Heaviside functions.

I wrote the below code,

"restart:
simplify(Dirac(x)*f(x))".

Then, Maple return the answer of this code, "Dirac(x)*f(0)".
I could understand this solution, then I rewrote the next code,

"restart:
_EnvUseHeavisideAsUnitStep:=true:
Heaviside(0)"

Then, Maple returned "1", because I set the value of Heaviside's step function H(x) at x=0.
Finally, I wrote the below code, but there was a problem, I think.

"restart:
_EnvUseHeavisideAsUnitStep:=true:
simplify(Dirac(x)*Heaviside(x))"

The solution of Maple was "0".

According to my first code, I think this solution is Dirac(x)*Heaviside(0), that is, Dirac(x).
I can't understand this result.
Someone help me, please.

The command

plots:-implicitplot(evalc(argument((1+x+I*y)/(1-x-I*y))) <= (1/4)*Pi, x = -5 .. 5, y = -5 .. 5, crossingrefine = 1, gridrefine = 2, rational = true, filled, signchange = true, resolution = 1000);

produces an incorrect result

in view of

evalf(argument((1-4+4*I)/(1+4-4*I)));
                          2.889038378

There is a workaround 

plots:-inequal(evalc(argument((1+x+I*y)/(1-x-I*y))) <= (1/4)*Pi, x = -5 .. 5, y = -5 .. 5);

Dear all,

I need to transforme these equation from time domain to frequency domain with fourier transforms and solve it in frequency domain but i received the flowing error

any helps

thank you !

Code :

Download Fourier_TRAns_MAPLEprime.mwFourier_TRAns_MAPLEprime.mw

Bear with me, while I've used Maple V, I've never done any actual programming in Maple as I've almost entirely used Matlab.  Right now I've got an issue with trying to fill a square matrix using a sinc function [(sin(Pi*x)/(Pi*x)].  Here's my procedure:

Note that the arguments for sin and for the denominator are written as they are so as to define the domain of  the sinc function as [-m/2,m/2] along a row and [-n/2,n/2] along a column of the matrix I want to fill and so that I the arguments can have decimal values.  Also note the conditional statements there because when the sine argument and the denominator yields Pi*0, the function equals 1 despite the zero in the denominator.

So I'm unsure how to fill a matrix of some size m x n using the procedure (or a function if that makes the programming easier).  Here's how I would write the code for Matlab:

Z=zeros(2001,2001);
X=(-5:.005:5)';
Y=(-5:.005:5)';

count1=1;

for j=-5:.005:5
    count2=1;
    for k=-5:.005:5
        Z(count1,count2)=sinc(k)^2*sinc(j)^2;
        count2=count2+1;
    end
    count1=count1+1;
end

Given my meager experience writing Maple code, I'm a bit flumoxed.

Thanks, and don't feel like you have to be gentle. No one ever accused me of writing elegant code.

Thanks

The command 

restart; st := time(): FunctionAdvisor(EllipticE); time()-st;

produces the result on my comp in 805.484 s. Too much time.

Dears;

Hope everyone is fine. I am try to find the numerical solutions of system of nonlinear algabric equation via newton's raphson method in the attached file but failed. Please see the attachment and try to correct. You can solve it least square method if possible. I am waiting your positive response. 

Help_in_Newton.mw

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

Email: muhammadusman@pku.edu.cn

Both the commands 

maximize(x*sin(t)+y*sin(2*t), t = 0 .. 2*Pi)assuming x>=0,y>=0;
minimize(x*sin(t)+y*sin(2*t), t = 0 .. 2*Pi)assuming x>=0,y>=0;

output 0. Simply no words.

The following three commands 

plots:-implicitplot(3*cos(x) = tan(y)^3, x = -Pi .. Pi, y = -(1/2)*Pi-1 .. (1/2)*Pi+1, thickness = 3, crossingrefine = 1, rational = true, signchange = true, resolution = 1000, gridrefine = 2);
plots:-implicitplot(3*cos(x) = tan(y)^3, x = -Pi .. Pi, y = -(1/2)*Pi-1 .. (1/2)*Pi+1, thickness = 3, crossingrefine = 1, rational = true, signchange = false, resolution = 1000, gridrefine = 2);
plots:-implicitplot(3*cos(x) = tan(y)^3, x = -Pi .. Pi, y = -(1/2)*Pi-1 .. (1/2)*Pi+1, thickness = 3, crossingrefine = 1, rational = true, resolution = 1000, gridrefine = 2);

produce the same incorrect plot 

It is clear the sraight lines given by y=Pi/2 and y=-Pi/2 are superfluous. It should be noticed that the Mmma's ContourPlot command without any options produces a correct plot.

Hi, I'm trying to use Maple to construct some examples of symmetry solutions for certain nonlinear PDE's.  As a warm up, however, I'm working through the commands just for the heat equation in 3d: u[t]-u[x,x]-u[y,y]-u[z,z]=0 

I've gotten Maple to produce both determining equations for the symmetry infinitesimal generators via the DeterminingPDE command.  I've also gotten the command Infinitesimals to work too.

However, when I next use PDETools Invariants command, it correctly outputs invariants for most of the generator output of Infinitesimals EXCEPT it won't output anything for the simple rotation generators yd[x]-xd[y].  It will, however, output invariants if the rotation is between an independent and the dependent coordinate.

An example:
with(PDETools)
S:=[_xi[x]=y, _xi[y]=-x, _eta[u]=0]
Invariants(S,u(x,y))

*Above returns nothing, But if you instead have _xi[x]=x and _xi[y]=y then it returns the right invariants.

Thanks in advance!

Up to http://www.maplesoft.com/support/help/Maple/view.aspx?path=solve&term=solve

 Let us consider 

solve({x > -Pi, (tan(x)-tan(x)^2)^2-cos(x+4*tan(x)) = -1, x < Pi}, [x]);
                               []

We see the command omits the solution x=0 without any warning. It should be noticed that Mathematica solves it, outputting

{{x -> 0}, {x -> 0}}

and the warning

Solve::incs: Warning: Solve was unable to prove that the solution set found is complete.

One may draw a conclusion on her/his own.

Hey, this is not the I've had this encounter. I want to open this saved document but when I open it and Maple starts up it just hits me with "A problem was encountered while opening the workbook. Database is not opened". How can I get to open it properly and see my math notes?

How can this be prevented?

When I try to upload the file in this message it says "Cant open a null file"

Any help?

Jacob

In manual solution is no problem, but i am interest to compute it with any software such maple since i am not familiar, How to find solution in term of parametric equation a(r), b(r), c(r) for r=-5..5 and also visualize this three derivative condition? d(a)=a+b+c, d(b)=a-b+c, d(c)=a-b-c

Does anyone know how to take a Maple plot and convert it to FITS format?  FITS is commonly used in astronomical images and stands for Flexible Image Transport System.

Thank you very much !