Maple 2016 Questions and Posts

These are Posts and Questions associated with the product, Maple 2016

I have this H(alpha,beta,eta) con ‡uent hypergeometric function, which I can easily adopt in Mathematica like this

Hypergeometric1F1[( (\[Alpha] + S \[Beta]))/\[Beta]^2, 
 1 + ((\[Alpha] + S \[Beta]))/\[Beta]^2, -\[Alpha] ((
   E^(-\[Beta] \[Eta]))   )/\[Beta]^2]

But when I try to input this in maple using the built-in `hypergeom` it gives an error.

hypergeom((S*lambda2+alpha)/lambda2^2, 1+(S*lambda2+alpha)/lambda2^2, -alpha*exp(-lambda2*eta)/lambda2^2)

Any suggestions?

 

Hi all,

I have this equation that I can not get all solutions symbolically:

restart:

eq1 := cos(lambda*ln(r1))*cos(lambda*ln(r2))+sin(lambda*ln(r1))*sin(lambda*ln(r2))-1 = 0:

solve(eq1, lambda, allsolutions) assuming r1>0, r2>0, r2>r1

when r1:=1: r2:=2: I get the solution

2*Pi*Z/ln(2)

when r1:=1.1: and r2:=2.1: # no solutions

How to get symbolique solution

Thanks

 

 

I'm trying to approximate the solution of an IVP using Euler's method in the InitialValueProblem command and I keep getting this error (see attached worksheet). Can someone explain why? Thanks!
 

restart

with(Student[NumericalAnalysis]):

InitialValueProblem((D(y))(t) = t*y(t)+1/y(t), y(0) = 3, t = 2, method = euler, numsteps = 5, output = solution)

Error, (in dsolve/numeric) array output cannot be obtained for problems containing global variables

 

``


 

Download ivp_error.mw

I have noticed the second item in Description in

 

http://www.maplesoft.com/support/help/Maple/view.aspx?path=convert/binary

 

The item is: 

“The binary number is returned as a base 10 number consisting of the digits 1 and 0 only.”

 

Does the existence of this item mean the following behavior will not be corrected ? ("10^(-6)" in the following should be "2^(-6)" .)

******************************************************

 

When trying to perform the following:

p1 := proc (x, y) if x^2+y^2 <= 1 then x*y-y^2 else 0 end if end proc;
plot3d(p1, x = -1 .. 1, y = -1 .. 1);
Error, (in plot3d) expected ranges but received x = -1 .. 1 and y = -1 .. 1
 

I get this strange error message. To the best of my knowledge x and y ARE provided as ranges. What am I missing/not understanding?

If I omit the ranges in plot3d Maple returns a correct plot, but the default range (-10 .. 10) does not display sufficient details

I'd like to make a graph like the below. So I know that using the display command and putting graphs in matrix form produces something similar (display(<graph1|graph2|graph3>), however it doesn't allow me to export it as one graph. Any ideas?

 

 

Hello folks, i recently had really confusing times when i plotted the function (x^(5/3)-5*x^(2/3)) in maple 2016. I know how this graph looks like as i have it printed in my undergraduate textbook. But in maple i'm not getting function(above given) values plotted for negative numbers. When i evaluate the function for x=-1 i'm getting some complex numbers, when i'm supposed to have simply -6. Could anybody tell me what's going on?plot(x^(5/3)-5*x^(2/3))

Respected member!
Please help me to find the solution of attached problem, I am a new user so pleaes forgive any mistakes.
 

``


``


``

NULL

NULL

restart

R := 2.0

2.0

(1)

ODEforNum := r^3*((D@@4)(F))(r)+r^2*R*((D@@3)(F))(r)*(F(r)-2/R)+R*((D(F))(r)-r*((D@@2)(F))(r))*(r*(D(F))(r)+3*(F(r)-1/R)) = 0:

numsol := dsolve({BCSforNum, ODEforNum}, numeric, output = listprocedure)

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

``


 

Download mplprimes.mw

Hello.

I have a Pde solution in from of the sum.

pde := diff(u(x, t), t) = diff(u(x, t), x$2)

symbolic := pdsolve([pde, u(x, 0) = 1, u(0, t) = 0, u(1, t) = 0])

symbolic := u(x, t) = Sum(-(2*((-1)^_Z9-1))*sin(_Z9*Pi*x)*exp(-Pi^2*_Z9^2*t)/(Pi*_Z9), _Z9 = 1 .. infinity)

 

I tried a subs or eval command dosen't work.

 

Thanks.

pdex1.mw
 

restart

pde := diff(u(x, t), t) = diff(u(x, t), `$`(x, 2)):

ics := [u(x, 0) = 1, u(0, t) = 0, u(1, t) = 0]:

pds := pdsolve(pde, ics, numeric, time = t, range = 0 .. 1, spacestep = 1/4024, timestep = 1/4024):

symbolic := pdsolve([pde, u(x, 0) = 1, u(0, t) = 0, u(1, t) = 0])

u(x, t) = Sum(-2*((-1)^_Z9-1)*sin(_Z9*Pi*x)*exp(-Pi^2*_Z9^2*t)/(Pi*_Z9), _Z9 = 1 .. infinity)

(1)

eval(rhs(symbolic), `~`[_Z9] = n)

Sum(-2*((-1)^_Z9-1)*sin(_Z9*Pi*x)*exp(-Pi^2*_Z9^2*t)/(Pi*_Z9), _Z9 = 1 .. infinity)

(2)

subs(`~`[_Z9] = n, rhs(symbolic))

Sum(-2*((-1)^_Z9-1)*sin(_Z9*Pi*x)*exp(-Pi^2*_Z9^2*t)/(Pi*_Z9), _Z9 = 1 .. infinity)

(3)

subs[eval](`~`[_Z9] = n, rhs(symbolic))

Sum(-2*((-1)^_Z9-1)*sin(_Z9*Pi*x)*exp(-Pi^2*_Z9^2*t)/(Pi*_Z9), _Z9 = 1 .. infinity)

(4)

``


 

Download pdex1.mw

 

The piecewise plot below displays a sphere truncated by the plane z = 2 - y.

f := proc (x, y, z) options operator, arrow; piecewise(z <= 2-y, x^2+y^2+z^2-16, z-2+y) end proc; implicitplot3d(f, -4 .. 4, -4 .. 4, -4 .. 4, style = surface, numpoints = 50000);

The 3 transforms below when executed in display(T(sphere([0, 0, 0], 4, numpoints = 50000)), scaling = constrained) display the truncated sphere differently:

1) the truncating plane only partly conforms to the boundary of the truncated sphere

2) the truncated sphere is correct provided that the else condition coordinate is in the truncating plane and truncated sphere

3) the truncated sphere is correct but hollow

 

1) T := transform(proc (x, y, z) options operator, arrow; `if`(z <= 2-y, [x, y, z], [x, y, 2-y]) end proc)

2) T := transform(proc (x, y, z) options operator, arrow; `if`(z <= 2-y, [x, y, z], [0, 0, 2]) end proc)

3) T := transform(proc (x, y, z) options operator, arrow; `if`(z <= 2-y, [x, y, z], [`&+-`(sqrt(16-y^2-z^2)), y, 2-y]) end proc)

Please explain the different behavior of the three transforms.           

Good evening sir.

 

I request your valuable support with regard to the above cited query.

 

 

With thanks & regards.

 

Mr.M.Anand

Associate Professor in Mathematics

I need hepl.  I work with the physics paсkage and I set:

with(Physics);
with(StandardModel);
Setup(mathematicalnotation = true)

Coordinates(X)

Define(p)

Setup(su2)

CompactDisplay(p(X))

U := 1+i*(1/f)*p[a](X)*Psigma[a]

H := v*U

DD[mu] := (d_[mu]+2*i*nu_5*KroneckerDelta[mu, 0])*Psigma[0]

And next:

Trace(DD[mu]*H*DD[mu]*H) or simplify(DD[mu]*H*DD[mu]*H)

Maple speaks:

Error, (in Physics:-Trace) invalid input: `union` received Physics:-d_[mu] = F1, which is not valid for its 1st argument

What I am doing wrong? (f, v, nu_5 is constans; a=1,2,3 and mu=0,1,2,3, Psigma[0] is unit matrix 2x2)

 

And if I write:

simplify(Trace(Psigma[a]*Psigma[a]))

Maple doesn't understand that it equals six. 

 

Let us consider 

Student[Precalculus]:-LimitTutor(sqrt(x), x = 2);

One expects a nice illustration of the result sqrt(2). But instead of that one reads "f(x) approaches 1.41 as x approaches 2". This is simply clueless and forms a wrong understanding of limits. It should also be noticed that all the entries (left, 2-sided, and right) produce the same animation. The same issue with other limits I tried, e.g.

Student[Precalculus]:-LimitTutor(sqrt(x), x = 1);

. I think this command should be completely rewritten or excluded from Maple. 

Hi

I have a set P which is filled with 2D elements e.g.

P:={[1,1],[1,2],[1,3]...}


Later I'm wanting to add a new element e.g. {[1,10]}
At the moment I have

P union {[10, 10]}


If I were to print this I get what I want, but I want to update P.
When I try and do

 P:= P union {[10,10]}

I get an error "Error, illegal use of an object as a name"

What am I doing wrong and how can I solve this? Much appreciated!

Hi,

I'm just doing a small project/HW for school where I am trying to create a histogram of the appearance of the digits of Pi. My problem is that I lose connection with the Maple Kernel whenever I go over 1000 digits. Why? Is this because my function scales badly? Any help would be appreciated.

Thanks in advance.

d:=1500;
Digits:=d;
                              1500
                              1500
valpi:=evalf(Pi):

dg:=proc(n,k)
floor(n/10^(ceil(log10(n)-k))) mod 10;
end:

with(Statistics):
L:=[seq(dg(valpi,i),i=1..d)]:
Histogram(L)

 

First 37 38 39 40 41 42 43 Last Page 39 of 60