Advanced Search

Maple Questions and Posts

These are Posts and Questions associated with the product, Maple
View unanswered posts
Maple Questions and Posts Feed

Maple Transactions Autumn Issue 2024 is now up

by: rcorless 735 Maple
October 31 2024
8  0

Maple Transactions has just published the Autumn 2024 issue at mapletransactions.org

From the header:

This Autumn Issue contains a "Puzzles" section, with some recherché questions, which we hope you will find to be fun to think about.  The Borwein integral (not the Borwein integral of XKCD fame, another one) set out in that section is, so far as we know, open: we "know" the value of the integral because how could the identity be true for thousands of digits but yet not be really true? Even if there is no proof.  But, Jon and Peter Borwein had this wonderful paper on Strange Series and High Precision Fraud showing examples of just that kind of trickery.  So, we don't know.  Maybe you will be the one to prove it! (Or prove it false.)

We also have some historical papers (one by a student, discussing the work of his great grandfather), and another paper describing what I think is a fun use of Maple not only to compute integrals (and to compute them very rapidly) but which actually required us to make an improvement to a well-known tool in asymptotic evaluation of integrals, namely Watson's Lemma, just to explain why Maple is so successful here.

Finally, we have an important paper on rational interpolation, which tells you how to deal well with interpolation points that are not so well distributed.

Enjoy the issue, and keep your contributions coming.

Maple and MapleSim 2024.2 updates

by: TechnicalSupport 790 Maple 2024 , MapleSim 2024
October 31 2024
3  0

We have just released updates to Maple and MapleSim.

Maple 2024.2 includes ability to tear away tabs into new windows, improvements to scrollable matrices, corrections to PDF export, small improvements throughout the math engine, and more.  We recommend that all Maple 2024 users install this update.

This update also include a fix to the problem with the simplify extension mechanism, as first reported on MaplePrimes. Thanks, as always, for helping us make Maple better.

This update is available through Tools>Check for Updates in Maple, and is also available from the Maple 2024.2 download page, where you can find more details.

At the same time, we have also released an update to MapleSim, which contains a variety of improvements to MapleSim and its add-ons. You can find more information on the MapleSim 2024.2 download page.

Content of exported SVG files is always outside th...

Asked by: mafitch 10
export svg
October 31 2024
2  0

When I export images as SVG the resulting file always has the image much larger than the viewbox resulting in only part of it showing. My normal workflow is to then load in Inkscape and correct the error, however, I would love to not have an extra step. Is there any other fix for this?

How do you convert this system of equations into m...

Asked by: GFY 65
matrix ode differential-equations
October 31 2024
0  3

How do you convert this system of equations into matrix form? The decoupling process is performed.Convert the equations into photographs.fuxian1030.mw

Change of coordinates: Are there other ways to loo...

Asked by: C_R 3437
coords
October 31 2024
1  2

What I do at the moment

u-> changecoords(u,[x,y,z],spherical_physics,[r,theta,phi]):
map(%,[x,y,z])
  [r cos(phi) sin(theta), r sin(phi) sin(theta), r cos(theta)]

Any other (preferably shorter) ways to look up transformations defined in ?coords

How do i define a tensor with an equation of three...

Asked by: Steve255 15
tensor differentialgeometry physics
October 30 2024
3  4

I am studying linear operators that have vanishing Nijenhuis Torsion, that is a (1-1) tensor L whose corresponding (1-2) tensor N given by a tensor equation of L (see attached) is identically zero.

I am new to maple, i have used it to plot vector fields and solve systems of equations in the past but i am unfamiliar with the DifferentialGeometry and Physics packages.

Attached is my best effort at solving this problem directly for a simple 2d case without the use of any packages and i am wondering if it is possible to do it all in one line without having to define tensor components one by one.

nijenhuis_torsion.mw

Solving Nonlinear Systems Symbolically in Maple...

Asked by: repromag 15
solve step-by-step
October 30 2024
0  0

How can I use Maple to solve a system of nonlinear equations symbolically and display the steps in the solution?

Which function is more credible?...

Asked by: sursumCorda 1244
polynomial realdomain
October 30 2024
1  0

Here are two systems over the reals:

sys__1:=And(r*(387*r+52)+2<r*(226*q+121*s)+9*q*(q*(2*q-5)-3*s+2)+6*s,4*q^3+r*(27*r+4)+s^2=q*(q+18*r),q>=0,r>=0):
sys__2:=And((392-1739*q)*r+4*(2-9*q)**2+2151*r**2<75*r*s,4*q**3+r*(27*r+4)+s**2=q*(q+18*r),q>=0,r>=0):

The following results indicate that both and are satisfiable 

QuantifierElimination:-QuantifierEliminate(exists([s,q,r],sys__1));
                              true

QuantifierElimination:-QuantifierEliminate(exists([s,q,r],sys__1));
                              true

but RealDomain:-solve simply returns an empty list (that is, no solution exists) in both cases

RealDomain:-solve(sys__1,[q,s,r]); # ⟹ sys1 cannot be satisfied
                               []

RealDomain:-solve(sys__2,[q,s,r]); # ⟹ sys2 cannot be satisfied
                               []

As discussed in the previous problem, in contrast to using QuantifierElimination:-QuantifierEliminate, the use of RealDomain:-solve is unsafe. Nevertheless, the above output suggests that even the much-more-sophisticated QuantifierElimination:-QuantifierEliminate is still not always reliable (since the correct returnedvalue appears to be in lieu of ). So, what is the right command to handle polynomial systems over real domains in Maple? 

Standard form ode ...

Asked by: janhardo 700
ode differential-equations duplicate-question
October 30 2024
0  8

Can i get this ode in a "standardform"  ?

verg:= (-delta*eta^2 + alpha*eta)*diff(diff(U(xi), xi), xi) - U(xi)*(2*eta*gamma*theta*(delta*eta - alpha)*U(xi)^2 + eta^2*delta*k^2 + (-alpha*k^2 - 2*delta*k)*eta + 2*k*alpha + delta) = 0;

How to simplify the following term?...

Asked by: Alfred_F 335
step-by-step radical simplification
October 30 2024
1  3

Since the puzzle task "A circle is to be disturbed ..." makes no fun, here is a Maple task:
The term to be simplified step by step:
(2+10/(3*sqrt(3)))^(1/3)+(2-10/(3*sqrt(3)))^(1/3)

How can I use SPECTRA.mla in Maple on Windows?...

Asked by: bathudaide 25
package fgb spectra
October 29 2024
0  12

I have learned that the SPECTRA.mla library can solve SDP problems. I have tried to download and use it, but I am still missing the FGb module

> with(SPECTRA);
> M := Matrix(6, 6, [[10, 1, 0, m[1], -m[3], m[2]], [1, -2*m[1]+27, m[3], -27/2, -m[4], -m[5]], [0, m[3], -2*m[2], m[4], m[5], 0], [m[1], -27/2, m[4], 10, 0, m[6]], [-m[3], -m[4], m[5], 0, -2*m[6], 0], [m[2], -m[5], 0, m[6], 0, 1]]);
> SolveLMI(M);
=> Error, (in SPECTRA:-SolveLMI) `FGb` does not evaluate to a module

 

How can I resolve this issue on a Windows environment, beause I don't see install file for Windows, only MacOS and Linux:

Which of the following is a possible value of 𝑥?...

Asked by: yangtheary 290
gcd recreational-mathematics
October 29 2024
0  4

CUDA - Can I use it to speed up symbolic computati...

Asked by: Calin 10
performance linearalgebra cuda
October 28 2024
2  1

We are working to obtain a fully symbolic dynamic model for a robot. 

Using the CPU (I9-12900K) and 128 GB of DDR5 RAM (5600 MHZ) did not compute a 7x7 Inverse of a symbolic matrix

Is is possible to exploit the CUDA functions to compute it on the GPU? I have a NVidia RTX A6000 (48 GB of DDR6 GPU memory)

I tried this: 

CUDA:-Enable(true)

CUDA: -MatrixInverse(D_Q):

But it does not use the GPU to compute this. 

Maybe I'm doing smth wrong. 

Thank you, 

Calin

2F1 hypergeometric...

Asked by: masa13 15
mathematica import hypergeom
October 28 2024
3  8

Hi

Dear friends, I am a relatively new user and I have a problem in entering and calculating the 2F1 hypergeometric function. My question is how to enter this function in Maple in equations so that Maple recognizes it? Because hypergeom ([1], [2], [3]) Maple itself is a 3-element function, while 2F1 hypergeometric ([1], [2], [3], [4]) is a four-element!

Confusing of methods for odetest?...

Asked by: salim-barzani 1575
ode series differential-equations
October 28 2024
1  11

Almost i did 10 method for this ode equation all of them are succes but this one is giving me some confusing and i am looking for  get my answer, the mothod say if we have the auxilary equation if substitute the solution of this auxilary equation in our series solution then substitute in ode equation must be satisfy but it is not satisfy so when he did assumption for the auxilary equation he say it satisfy if we sabstitute this assumption in our series solution!

My question is this how we get thus assumption ? and why finding exact  solution of auxilary equation not satisfy?

> 

restart
> 

with(PDEtools)
> 

with(LinearAlgebra)
> 

with(Physics)
> 

with(SolveTools)
> 

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

 (1)
> 

NULL
> 

Fode := (-delta*eta^2+alpha*eta)*(diff(diff(U(xi), xi), xi))-U(xi)*(2*eta*gamma*theta*(delta*eta-alpha)*U(xi)^2+eta^2*delta*k^2+(-alpha*k^2-2*delta*k)*eta+2*k*alpha+delta) = 0

(-delta*eta^2+alpha*eta)*(diff(diff(U(xi), xi), xi))-U(xi)*(2*gamma*eta*theta*(delta*eta-alpha)*U(xi)^2+eta^2*delta*k^2+(-alpha*k^2-2*delta*k)*eta+2*k*alpha+delta) = 0

 (2)
> 

NULL
> 

F := sum(a[i]*G(xi)^i, i = 0 .. 1)

a[0]+a[1]*G(xi)

 (3)

``

 (4)
> 

D1 := diff(F, xi)

a[1]*(diff(G(xi), xi))

 (5)
> 

NULL
> 

S := (diff(G(xi), xi))^2 = G(xi)^4+A[2]*G(xi)^2+A[1]

(diff(G(xi), xi))^2 = G(xi)^4+A[2]*G(xi)^2+A[1]

 (6)
> 

S1 := diff(G(xi), xi) = sqrt(G(xi)^4+A[2]*G(xi)^2+A[1])

diff(G(xi), xi) = (G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (7)
> 

E1 := subs(S1, D1)

a[1]*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (8)
> 

D2 := diff(E1, xi)

(1/2)*a[1]*(4*G(xi)^3*(diff(G(xi), xi))+2*A[2]*G(xi)*(diff(G(xi), xi)))/(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (9)
> 

E2 := subs(S1, D2)

(1/2)*a[1]*(4*G(xi)^3*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)+2*A[2]*G(xi)*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2))/(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (10)
> 

K := U(xi) = F

U(xi) = a[0]+a[1]*G(xi)

 (11)
> 

K1 := diff(U(xi), xi) = E1

diff(U(xi), xi) = a[1]*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (12)
> 

K2 := diff(U(xi), xi, xi) = E2

diff(diff(U(xi), xi), xi) = (1/2)*a[1]*(4*G(xi)^3*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)+2*A[2]*G(xi)*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2))/(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (13)
> 

L := eval(Fode, {K, K1, K2})

(1/2)*(-delta*eta^2+alpha*eta)*a[1]*(4*G(xi)^3*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)+2*A[2]*G(xi)*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2))/(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)-(a[0]+a[1]*G(xi))*(2*gamma*eta*theta*(delta*eta-alpha)*(a[0]+a[1]*G(xi))^2+eta^2*delta*k^2+(-alpha*k^2-2*delta*k)*eta+2*k*alpha+delta) = 0

 (14)
> 

L1 := normal((1/2)*(-delta*eta^2+alpha*eta)*a[1]*(4*G(xi)^3*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)+2*A[2]*G(xi)*(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2))/(G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)-(a[0]+a[1]*G(xi))*(2*gamma*eta*theta*(delta*eta-alpha)*(a[0]+a[1]*G(xi))^2+eta^2*delta*k^2+(-alpha*k^2-2*delta*k)*eta+2*k*alpha+delta) = 0)
> 

 

> 

collect(L1, {G(xi)})

(-2*delta*eta^2*gamma*theta*a[1]^3+2*alpha*eta*gamma*theta*a[1]^3-2*delta*eta^2*a[1]+2*alpha*eta*a[1])*G(xi)^3+(-6*delta*eta^2*gamma*theta*a[0]*a[1]^2+6*alpha*eta*gamma*theta*a[0]*a[1]^2)*G(xi)^2+(-6*delta*eta^2*gamma*theta*a[0]^2*a[1]+6*alpha*eta*gamma*theta*a[0]^2*a[1]-delta*eta^2*k^2*a[1]+alpha*eta*k^2*a[1]-delta*eta^2*A[2]*a[1]+alpha*eta*A[2]*a[1]+2*delta*eta*k*a[1]-2*alpha*k*a[1]-delta*a[1])*G(xi)-2*gamma*delta*eta^2*theta*a[0]^3+2*gamma*alpha*eta*theta*a[0]^3-delta*eta^2*k^2*a[0]+alpha*eta*k^2*a[0]+2*delta*eta*k*a[0]-2*alpha*k*a[0]-delta*a[0] = 0

 (15)
> 

eq0 := -2*delta*eta^2*gamma*theta*a[0]^3+2*alpha*eta*gamma*theta*a[0]^3-delta*eta^2*k^2*a[0]+alpha*eta*k^2*a[0]+2*delta*eta*k*a[0]-2*alpha*k*a[0]-delta*a[0] = 0
> 

eq1 := -6*delta*eta^2*gamma*theta*a[0]^2*a[1]+6*alpha*eta*gamma*theta*a[0]^2*a[1]-delta*eta^2*k^2*a[1]+alpha*eta*k^2*a[1]-delta*eta^2*A[2]*a[1]+alpha*eta*A[2]*a[1]+2*delta*eta*k*a[1]-2*alpha*k*a[1]-delta*a[1] = 0
> 

eq2 := -6*delta*eta^2*gamma*theta*a[0]*a[1]^2+6*alpha*eta*gamma*theta*a[0]*a[1]^2 = 0
> 

eq3 := -2*delta*eta^2*gamma*theta*a[1]^3+2*alpha*eta*gamma*theta*a[1]^3-2*delta*eta^2*a[1]+2*alpha*eta*a[1] = 0
> 

COEFFS := solve({eq0, eq1, eq2, eq3}, {alpha, eta, a[0], a[1]}, explicit)
> 

case1 := COEFFS[4]

{alpha = delta*(eta^2*k^2+eta^2*A[2]-2*eta*k+1)/(eta*k^2+eta*A[2]-2*k), eta = eta, a[0] = 0, a[1] = 1/(-gamma*theta)^(1/2)}

 (16)
> 

NULL
> 

S

(diff(G(xi), xi))^2 = G(xi)^4+A[2]*G(xi)^2+A[1]

 (17)
> 

S1

diff(G(xi), xi) = (G(xi)^4+A[2]*G(xi)^2+A[1])^(1/2)

 (18)
> 

S2 := dsolve(S, G(xi))

G(xi) = -(1/2)*(-2*A[2]-2*(A[2]^2-4*A[1])^(1/2))^(1/2), G(xi) = (1/2)*(-2*A[2]-2*(A[2]^2-4*A[1])^(1/2))^(1/2), G(xi) = -(1/2)*(2*(A[2]^2-4*A[1])^(1/2)-2*A[2])^(1/2), G(xi) = (1/2)*(2*(A[2]^2-4*A[1])^(1/2)-2*A[2])^(1/2), G(xi) = JacobiSN((1/2)*(2*(A[2]^2-4*A[1])^(1/2)-2*A[2])^(1/2)*xi+c__1, (-2*(A[2]*(A[2]^2-4*A[1])^(1/2)-A[2]^2+2*A[1])*A[1])^(1/2)/(A[2]*(A[2]^2-4*A[1])^(1/2)-A[2]^2+2*A[1]))*A[1]*2^(1/2)/(A[1]*(-A[2]+(A[2]^2-4*A[1])^(1/2)))^(1/2)

 (19)
> 

K

U(xi) = a[0]+a[1]*G(xi)

 (20)
> 

K4 := subs(case1, K)

U(xi) = G(xi)/(-gamma*theta)^(1/2)

 (21)
> 

NULL
> 

K5 := subs(S2, K4)

U(xi) = -(1/2)*(-2*A[2]-2*(A[2]^2-4*A[1])^(1/2))^(1/2)/(-gamma*theta)^(1/2)

 (22)
> 

NULL
> 

F2 := subs(case1, Fode)

(-eta^2*delta+delta*(eta^2*k^2+eta^2*A[2]-2*eta*k+1)*eta/(eta*k^2+eta*A[2]-2*k))*(diff(diff(U(xi), xi), xi))-U(xi)*(2*gamma*eta*theta*(delta*eta-delta*(eta^2*k^2+eta^2*A[2]-2*eta*k+1)/(eta*k^2+eta*A[2]-2*k))*U(xi)^2+eta^2*delta*k^2+(-k^2*delta*(eta^2*k^2+eta^2*A[2]-2*eta*k+1)/(eta*k^2+eta*A[2]-2*k)-2*k*delta)*eta+2*k*delta*(eta^2*k^2+eta^2*A[2]-2*eta*k+1)/(eta*k^2+eta*A[2]-2*k)+delta) = 0

 (23)

``

 (24)
> 

NULL
> 

odetest(K5, F2)

-(1/2)*delta*eta*(A[2]^2-4*A[1])^(1/2)*(-2*(A[2]+(A[2]^2-4*A[1])^(1/2))/gamma)^(1/2)/((eta*k^2+eta*A[2]-2*k)*(-theta)^(1/2))

 (25)


and i hope mapleprimes don't delete this question becuase of this pictures also it help for undrestanding

 

there is other picture for different auxilary equation just  add one multiply term for G(xi)^4 in case anyone needed i will upload

Download odetest.mw

First 61 62 63 64 65 66 67 Last Page 63 of 2218