Students discount...

Please may I know if you can offer ma student's discount ob the seleted version.

Thank you.

Fred.

How can I get Maple to find the Fourier transform ...

I'm a Maple novice.  I have two questions.  (1) I'm trying to use Maple to confirm some Fourier transforms of selected probability density functions.  In general I've succeeded, but Maple fails to find the Fourier transform of the pdf of a logistic random variable with mean 0.  Please explain how I can get Maple to carry out this request.  I've attached a Maple file with my work to this question.  (2) At several points in my computations, I wish to substitute 2*Pi*xi for omega.  I have an expression containing two omegas.  If I use algsubs(), only one of the two omegas is replaced.  I have to use subs() to replace both omegas.  Why is this?

Use Maple to confirm selected Fourier transform of logistic random variable

 >
 (1)
 >
 >
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)

Wikipedia's article on "Logistic distribution" gives a characteristic function that implies that the Fourier transform of this pdf should equal Pi*a*omega/sinh(Pi*a*omega).

Unit rectangular function

 >
 (8)
 >
 (9)
 >
 (10)
 >
 (11)
 >

Intersection of 3D-plane with a 3D-drop...

How do I calculate the intersection curve between a plane and a drop?
The "drop" is defined in the following way:

R1 := 3.;R2 := 1.0;DR := 4;g := R2 + DR;

f1 := h -> sqrt(R1^2 - h^2);
f2 := h -> sqrt(g^2 - h^2);
f3 := h -> (1 - h/g)*f1(h*R1/g) + h*f2(h)/g;
f4 := h -> sqrt(1/2*g - 1/2*h);
f5 := h -> (1 - h/g)*f3(h) + h*f4(h)/g;
gg := h -> piecewise(h < 0, f1(h), 0 <= h, f5(h));#Radius depending on the z-position h
cir := (h, phi, R) -> <sin(phi)*R, cos(phi)*R, h>; # a circle at the hight h with radius R
#The plane is placed inside the drop.
n := (x, y, z) -> <x, y, z>/sqrt(x^2 + y^2 + z^2);

# the following lines show, how it looks like:
with(plots);
with(plottools);
dro1 := plot3d(cir(h, phi, gg(h)), h = -R1 .. g, phi = 0 .. 2*Pi, scaling = constrained, orientation = [-60, 72, 0]);
plotDropWithPlane := (x, y, z) -> display(dro1, arrow(Vector([0, 0, 0]), 2*R1*n(x, y, z), 0.2, 0.4, 0.1, cylindrical_arrow, fringe = blue, color = "Green"), implicitplot3d(x*x1 + y*y1 + z*z1 = 0, x1 = -R1 .. R1, y1 = -R1 .. R1, z1 = -R1 .. g, color = blue));
plotDropWithPlane(3, 1, 2);

#I'm searching a function like
fintersect:=theta-><"?,?,?>"

Any idea how to solve?

Best regards,

Andreas

Derivative with numerical method...

 >
 >
 >
 >
 >
 (1)
 >
 (2)
 >
 (3)

Using numerical methods, I cannot calculate the derivative of sd with respect to r at r = 2 and t = 4.

Oliveira.

 >

Dirichlet conditions and Neumann values...

Does anyone know how to enter in the pdsolve function Dirichlet conditions and Neumann values?

Oliveira.

dsolve returns signular solution which can be obta...

Why Maple returns -1/x as singular solution below when this solution can be obtained from the general solution when constant of integration is zero?

 > restart;
 > ode:=2*y(x)+2*x*y(x)^2+(2*x+2*x^2*y(x))*diff(y(x),x) = 0; dsolve(ode,singsol=false);

 > sol:=[dsolve(ode,singsol=essential)];

 > subs(_C1=0,sol)

 >

Loop two outputs...

I've made this proc and it has 2 outputs

*

How do I fix that?

strange latex result. Why an extra 1 shows up?...

We all know that Maple's Latex is not the best of Maple to say the least.

But this one is really strange. Maple prints a 1 for no apparant reason in the latex which makes it ugly.

I wonder if Maplesoft still maintains its Latex conversion code at all?  So one can at least hope may be one day all of this will get fixed? Year after year, and Maple's Latex still not changed.

If Mapesoft do not intend to make any changes in its Latex conversion software at all, it will be good if an official statement is made in this regards so that at least customers know.

 > sol:=dsolve((x-a)*(x-b)*diff(y(x),x)+k*(y(x)-a)*(y(x)-b) = 0,y(x)): sol:=subs(_C1=C[1],sol);

 > latex(sol)

y \left( x \right) ={1 \left(  \left( x-b \right) ^{-k} \left( x-a
\right) ^{k}a{{\rm e}^{akC_{{1}}-bkC_{{1}}}}- \left( x-b \right) ^{-k
} \left( x-a \right) ^{k}b{{\rm e}^{akC_{{1}}-bkC_{{1}}}}+b \left( {
\frac {-x+b}{-x+a}} \right) ^{-k}{{\rm e}^{akC_{{1}}-bkC_{{1}}}}-b

\right)  \left( -1+ \left( {\frac {-x+b}{-x+a}} \right) ^{-k}{{\rm e}
^{akC_{{1}}-bkC_{{1}}}} \right) ^{-1}}

 >

why odetest fail to verify solution when using C[1...

Why odetest sometimes fail to give 0  from odetest when simply using C[1] instead of _C1 as constant of integration?

I do not remember now if I asked about this before now. But for me as a user, this does not look right. I like to use C[1] instead of _C1 as the constant of integration as it looks better in Latex. I had no idea it will make a difference to odetest what the constant of integration symbol used is.

Is this a known issue? Do you consider this a bug? Maple 2019.1 on windows 10.

 > restart;
 > ode:= x^2*diff(y(x),x)+x*y(x)+sqrt(y(x)) = 0;

 > sol_1:=sqrt(y(x))=1/x+_C1/sqrt(x); odetest(sol_1,ode)

 > sol_2:=subs(_C1=C[1],sol_1); odetest(sol_2,ode); #why this now fails??

 > sol_3:=subs(C[1]=_C1,sol_2); odetest(sol_3,ode)

A=B but Maple does not want to simplify arctanh(A...

A=B  but not able to simplify arctanh(A)-arctanh(B)  to be zero.  Why? Is there a workaround? Using Maple 2019.1

 > restart;
 > A:=((y*sqrt(3) + 3)*sqrt(3))/(6*sqrt(y^2 + 1)); B:=(y + sqrt(3))/(2*sqrt(y^2 + 1)); simplify(A-B)

 > simplify(arctanh(A)-arctanh(B))

 > simplify(arctanh(A)-arctanh(B),trig)

 > simplify(arctanh(A)-arctanh(B)) assuming positive

 > simplify(arctanh(A)-arctanh(B),trig) assuming positive

 > plot(arctanh(A),y=-Pi..Pi)

 > plot(arctanh(B),y=-Pi..Pi)

 >

Compare to Mathematica:

How can i make assumptions for all subscript-varia...

Hello,

In this simple example in Maple 2019:

k[1] := 16;

solve(k[2]^2 = k[1], {k[2]}, useassumptions) assuming 0 < k;

The result of the solve is

{k[2] = 4}, {k[2] = -4}

which is not what I'm expecting.

I'm expecting

{k[2] = 4}

I'm trying to assume that both the variable named k[1] and the variable named k[2] are greater than 0 by assuming that 0 < k.
My goal is a way to make the same assumption for variables with the same base but different subscripts.
How can I do this?

Combining and converting units...

 >

Note: To enter units, I used the unit key (blue) in the Units palette.

When I use the combine and simplify functions to manipulate temperature units, Maple returns a wrong answer, as we can see below;

 >
 (1)
 >
 (2)

The quantity 20 Celsius is not converted correctly. However, for other types of dimensions, the combine and simplify functions work correctly.

 >
 (3)
 >
 (4)
 >

Another strange thing happens with the convert function, with temperature units. When we use the same source unit, the convert function deletes the unit, leaving only the quantity.

 >
 (5)

I do not understand why this occurs.

 >
 (6)

I do not understand why this occurs.

Oliveira

lprint output (format)...

If I copy the output of a lprint command e.g.
lprint(<1,2;3,4>);

Matrix(2,2,{(1, 1) = 1, (1, 2) = 2, (2, 1) = 3, (2, 2) = 4},datatype = anything
,storage = rectangular,order = Fortran_order,shape = [])

and paste it in a new execution group, I obtain a "Line print output" instead of a Maple (1D) input,
so it cannot be executed.
(Of course, it's possible to paste in Notepad to remove the format, but it's annoying.)

This happens in Maple 2019,  Windows 64, Worksheet mode, 1D input;  but not all the time (sometimes it is as it used to be).
Do you see the same behavior?

When i type sqrt(25), i get 1.81847767202745*10^(-...

Yeah, i have tried evalf[10](sqrt(25)).
How can i get a simple number as answer? I'm loving the software but i just wished i could type in:
int(sin(x), x = 0 .. pi)
and get 2 instead of (2.739493386*10^(-116) + (2.739493386*10^(-116))*I)*pi.
Also, when i type evalf[50](pi), i wish to get all the 50 digits, but i just get \pi :/.

why rewriting the ODE changes it classification by...

This for me seems strange. A first oder ODE of degree one of form   A+B*y'=0. Asking Maple for its classification, it gives exact AND rational.  After rewriting it in the form y'=f(x,y) , now Maple says it is rational only. It no longer recognize it as exact.

How could this be possible? It is the same ODE, just rewrittng differently.  Maple 2019.1 on windows.

 > restart;
 > Typesetting:-Settings(typesetprime=true):
 > ode:=2*x/y(x)-y(x)/(x^2+y(x)^2)+(-x^2/y(x)^2+x/(x^2+y(x)^2))*diff(y(x),x) = 0; DEtools:-odeadvisor(ode);

 > dsolve(ode)

 > new_ode:=convert(ode,D); #so I can solve for y' RHS:=solve(new_ode,D(y)(x));