Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Throughout the course of a year, Maple users create wildly varying applications on all sorts of subjects. To mark the end of 2018, I thought I’d share some of the 2018 submissions to the Maple Application Center that I personally found particularly interesting.

Solving the 15-puzzle, by Curtis Bright. You know those puzzles where you have to move the pieces around inside a square to put them in order, and there’s only one free space to move into?  I’m not good at those puzzles, but it turns out Maple can help. This is one of collection of new, varied applications using Maple’s SAT solvers (if you want to solve the world’s hardest Sudoku, Maple’s SAT solvers can help with that, too).

Romeo y Julieta: Un clasico de las historias de amor... y de las ecuaciones diferenciales [Romeo and Juliet: A classic story of love..and differential equations], by Ranferi Gutierrez. This one made me laugh (and even more so once I put some of it in google translate, which is more than enough to let you appreciate the application even if you don’t speak Spanish). What’s not to like about modeling a high drama love story using DEs?

Prediction of malignant/benign of breast mass with DNN classifier, by Sophie Tan. Machine learning can save lives.

Hybrid Image of a Cat and a Dog, by Samir Khan. Signal processing can be more fun that I realized. This is one of those crazy optical illusions where the picture changes depending on how far away you are.

Beyond the 8 Queens Problem, by Yury Zavarovsky. In true mathematical fashion, why have 8 queens when you can have n?  (If you are interested in this problem, you can also see a different solution that uses SAT solvers.)

Gödel's Universe, by Frank Wang.  Can’t say I understood much of it, but it involves Gödel, Einstein, and Hawking, so I don’t need to understand it to find it interesting.

i have two ode equations, when try to solve, i face this error
Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system
any comment or help to solve this? ( or any simplifying comments to solve) . thnx in advance

 

restart

A := 4263736263.736263:B := 11868131.86813187:C := 600146.5201465202:hr := 0.3e7:muz := 0.2632654643708247e-5:P:=1000:

1000

(1)

ode1 := diff(u(r), r, r)+(diff(u(r), r))/r-u(r)/r^2-B*(diff(w(r), r, r, r)+(diff(w(r), r, r))/r-(diff(w(r), r))/r^2)/A = 0;

diff(diff(u(r), r), r)+(diff(u(r), r))/r-u(r)/r^2-0.2783505155e-2*(diff(diff(diff(w(r), r), r), r))-0.2783505155e-2*(diff(diff(w(r), r), r))/r+0.2783505155e-2*(diff(w(r), r))/r^2 = 0

(2)

ode2 := C*(diff(w(r), r, r, r, r))+2*C*(diff(w(r), r, r, r))/r-(C/r^2+muz*hr^2)*(diff(w(r), r, r))+(C/r^3-muz*hr^2/r)*(diff(w(r), r))+muz*hr^2*w(r)/r^2-B*(diff(u(r), r, r, r))-2*B*(diff(u(r), r, r))/r+B*(diff(u(r), r))/r^2-B*u(r)/r^3 = P;

600146.5201465202*(diff(diff(diff(diff(w(r), r), r), r), r))+1200293.040*(diff(diff(diff(w(r), r), r), r))/r-(600146.5201465202/r^2+23693891.80)*(diff(diff(w(r), r), r))+(600146.5201465202/r^3-23693891.80/r)*(diff(w(r), r))+23693891.80*w(r)/r^2-11868131.86813187*(diff(diff(diff(u(r), r), r), r))-23736263.74*(diff(diff(u(r), r), r))/r+11868131.86813187*(diff(u(r), r))/r^2-11868131.86813187*u(r)/r^3 = 1000

(3)

ics := w(.2) = 0, (D(w))(.2) = 0, w(1) = 0, (D(w))(1) = 0, u(.2) = 0, u(1) = 0;

w(.2) = 0, (D(w))(.2) = 0, w(1) = 0, (D(w))(1) = 0, u(.2) = 0, u(1) = 0

(4)

dsolve({ode1,ode2} union {ics},{w(r),u(r)},numeric)

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

 

 

 

Download zzv.mw

Hi!

I want to plot the approximation of a surface by polynomials. The surface is given by (x,y,f(x,y)) where f(x,y) is given by the following expression

proc (x) options operator, arrow; (sum(i*cos((i+1)*(-2+4*x[1])+i), i = 1 .. 5))*(sum(i*cos((i+1)*(-2+4*x[2])+i), i = 1 .. 5)) end proc

with both variables varying in the interval [0,1]. Then, by using the Bernstein polynomials of two variables (see, for instance, this paper for details  https://www.sciencedirect.com/science/article/pii/0021904589900956), the graph of the resulting (plot3d) surface (x,y,p(x,y))  it is not even like to the original surfaces.

Please, see this PDF of what I have done:  plots.pdf

Some idea or suggestion?

Thanks!

numeric.mws

 

Dear Colleagues,
Please give advice - How correct in Maple can be solve the heat nonlinear equation system with perfect thermal contact conditions?

Using procedure dsolve (numeric) Maple cannot determine type of problem (initial or boundary).

****
initial/boundary conditions specified at too many points 1
Error, (in dsolve/numeric/type_check) cannot determine type of problem (initial or boundary), as conditions are given at 4 points
****
 

In attachment Maple file - numeric.

Maybe anyone can get a simple exaple of the maple code for numeric solution of this  task?

 

 


 

restart;

k1:=5; k2:=7;k3:=8; Q:=4;p:=0.001;l:=0.1; d:=0.6;

k1 := 5

k2 := 7

k3 := 8

Q := 4

p := 0.1e-2

l := .1

d := .6

(1)

Eq1:=(k1^2)*diff(T1(x),x$2)=0;

 

Eq1 := 25*(diff(T1(x), `$`(x, 2))) = 0

(2)

Eq2:=(k2^2)*diff(T2(x),x$2)=-Q*Dirac(x)+T2(x)^4;

Eq2 := 49*(diff(T2(x), `$`(x, 2))) = -4*Dirac(x)+T2(x)^4

(3)

Eq3:=(k3^2)*diff(T3(x),x$2)=0;

Eq3 := 64*(diff(T3(x), `$`(x, 2))) = 0

(4)

con1:=T1(0)=T2(0);con2:=T3(p)=T2(p); con3:=k1*(D(T1)(0))=Q-k2*(D(T2)(0)); con4:=k2*(D(T2)(p))=k3*(D(T3)(p)); con5:= T3(l)=300; con6:=T1(-d)=300;

con1 := T1(0) = T2(0)

con2 := T3(0.1e-2) = T2(0.1e-2)

con3 := 5*(D(T1))(0) = 4-7*(D(T2))(0)

con4 := 7*(D(T2))(0.1e-2) = 8*(D(T3))(0.1e-2)

con5 := T3(.1) = 300

con6 := T1(-.6) = 300

(5)

nsol:=dsolve({Eq1, Eq2, Eq3,con1,con2,con3,con4,con5,con6},{T1(x),T2(x),T3(x)}, numeric,bvp);

 

 

initial/boundary conditions specified at too many points 1

Error, (in dsolve/numeric/type_check) cannot determine type of problem (initial or boundary), as conditions are given at 4 points

 

 

 

 


 

Download numeric.mws

I am brand new to Maple Cloud and Maple Player.

I have uploaded two worksheets to the cloud, and my wife has just installed Maple Player on her laptop.

In Maple Player, the second worksheet shows the shareable symbol but the first doesn't even thought I uploaded both in the same way by clicking on the upload symbol in the Maple Cloud palette. Why is the first worksheet not shareable?

When my wife displays the second worksheet she is able to move its sliders but they do not change the display as they do when I move the sliders within Maple2016. How can she change the display?

Here is a link to the second worksheet:

Cassinian_oval.mw

Say I have 2 matrices, in which the elements themselves are vectors.

I'm looking for a way to perform matrix multiplication on these so that rather than having the first element as x11y11+x12y21+x13y31

It would be x11.y11+x12.y21+x13.y31 where . is the dot product on the elements of each matrix.

I know I could write a procedure to do this manually but I was wondering if there's any pre-made operations (or modifiers on the Multiply operation) to do this.

How do I plot the optimal control functions in an optimal control problem ?

Hello,

I have been working on Maxima and minima, I am able to extract the eigen values for the expression.
Based on following conditions I am able to find out the critical point is maxima or minima or saddle or inconsistant

If all the eigenvalues are positive, the point is a minimum.
If all the eigenvalues are all negative, it's a maximum.
If some eigenvalues are positive, some are negative, and none are zero, then it's a saddle point.
If any eigenvalues are zero, the test is inconclusive

 

I want to return all the critical points and their extrema.
just for example : For one perticular function I got a Eigen values as which I can find using sign function.

EigenValues := [[-.381966011250105+0.*I, -2.61803398874989+0.*I], [.414213562373095+0.*I, -2.41421356237309+0.*I]]
signDetails := [seq([seq(sign(EigenValues[i][j]), j = 1 .. nops(EigenValues[i]))], i = 1 .. nops(EigenValues))] #

 signDetails :=[[-1, -1], [1, -1]]

Now if I have a 0 in a list. Sign function returns 1 for 0, which is incorrect. How can I handle such conditions

 

if I have

 

EigenValues := [[-.381966011250105+0.*I, -2.61803398874989+0.*I], [.414213562373095+0.*I, -2.41421356237309+0.*I], [0, 2]]

I would like to have output [[-1, -1], [1, -1], [0, 1]],

I would like to know how is it possible return output based on above list

in this case my return shouble something like this [maxima, saddle, inconclusive].

Thank you

Here is a simple procedure that works fine if entered using 1D Maple input
> Q:=proc(x)
sin(x)
end proc;
but if you use 2D math input
> q:=proc(x)
sin(x);

  end proc;

Error, unterminated procedure
    Typesetting:-mambiguous(qAssignTypesetting:-mambiguous(

      procApplyFunction(x) sinApplyFunction(x),

      Typesetting:-merror("unterminated procedure")))
Error, unable to parse
    Typesetting:-mambiguous(  Typesetting:-mambiguous(end,

      Typesetting:-merror("unable to parse")) procsemi)

Ouch! But to confuse things further the following procedures may be entered using 2D math and work fine:
>H := proc (x) x^2*sin(x) end proc;
>K := proc (x) sin(x^2) end proc;
Doesn't make any sense to me. Perhaps 2D math is not ready for prime time?

 

So I have got the following integral:

Int((2*z+1)/(z-5)(z-1)^2,z)

around the square with corners 2, 2i, -2, -2i oriented counter-clockwise.

Do I need to tell Maple that z is complex? Do I need to manually parametrize z? Is Maple aware of Cauchy's theorem?

What is the quickest way of evaluating this sort of integrals.

I'm still a maple beginner and looking for a workable solution to the following problem. I have already solved these "by hand" and would now like to control them.

The following equation should be solved:

diff(u(x, t), t, t)-c^2*(diff(u(x, t), x, x))

The approach function would be for example:

u := a(x)*sin(k1*(-c*t+x))

It would be great if maple could give me the solution again in this example:

solution = -c^2*(d^2*a*sin(k1*(-c*t+x))/dx^2+2*d*a*k1*cos(k1*(-c*t+x))/dx)

How exactly do I do this best?
Thanks a lot!
Frank

Hi,

How to simulate 1000 times the Dice1 experiment and put the results in a list to search for the percentils of A=[A1,..,A1000]

( See attachment)

Thanks

QuestionSimulation.mw

how can solve ordinary differenatial equation in maple?

how can use laplce transformation for this equation in maple?

U0 := proc (x) options operator, arrow; cosh(sqrt(2)*x)-1 end proc;
sys_ode := diff(Uc(x), x, x)-(2*(1+U0(x)))*Uc(x) = 0;



Download a.mw

 

 

 

Hi,

Persay I have a sinusodal function or dampened sinusodial; it is understand that the domain is from negative infinity to infinity.  Say I wish to find the zeros or a certain y value of a sinusodial on a restricted domain - say from A to B. How would I tell maple -syntax wise- to solve on that restricted domain? Using the solve feature on maple, it only yields 1 value that satisfies the condition I give it .

 

Thank you all for the support; you al l are allowing to me develop a better understanding for maple.

say I have the functions:

y=e^x*sin(x)

y=e

i wish to find all the x values that satisfy that system on the x domain of 0 to 100


Overview of the Physics Updates

 

One of the problems pointed out several times about the Physics package documentation is that the information is scattered. There are the help pages for each Physics command, then there is that page on Physics conventions, one other with Examples in different areas of physics, one "What's new in Physics" page at each release with illustrations only shown there. Then there are a number of Mapleprimes post describing the Physics project and showing how to use the package to tackle different problems. We seldomly find the information we are looking for fast enough.

 

This post thus organizes and presents all those elusive links in one place. All the hyperlinks below are alive from within a Maple worksheet. A link to this page is also appearing in all the Physics help pages in the future Maple release. Comments on practical ways to improve this presentation of information are welcome.

Description

 

As part of its commitment to providing the best possible environment for algebraic computations in Physics, Maplesoft launched, during 2014, a Maple Physics: Research and Development website. That enabled users to ask questions, provide feedback and download updated versions of the Physics package, around the clock.

The "Physics Updates" include improvements, fixes, and the latest new developments, in the areas of Physics, Differential Equations and Mathematical Functions. Since Maple 2018, you can install/uninstall the "Physics Updates" directly from the MapleCloud .

Maplesoft incorporated the results of this accelerated exchange with people around the world into the successive versions of Maple. Below there are two sections

• 

The Updates of Physics, as  an organized collection of links per Maple release, where you can find a description with examples of the subjects developed in the Physics package, from 2012 till 2019.

• 

The Mapleprimes Physics posts, containing the most important posts describing the Physics project and showing the use of the package to tackle problems in General Relativity and Quantum Mechanics.

The update of Physics in Maple 2018 and back to Maple 16 (2012)

 

 

• 

Physics Updates during 2018

a. 

Tensor product of Quantum States using Dirac's Bra-Ket Notation

b. 

Coherent States in Quantum Mechanics

c. 

The Zassenhaus formula and the algebra of the Pauli matrices

d. 

Multivariable Taylor series of expressions involving anticommutative (Grassmannian) variables

e. 

New SortProducts command

f. 

A Complete Guide for Tensor computations using Physics

 

• 

Physics Maple 2018 updates

g. 

Automatic handling of collision of tensor indices in products

h. 

User defined algebraic differential operators

i. 

The Physics:-Cactus package for Numerical Relativity

j. 

Automatic setting of the EnergyMomentumTensor for metrics of the database of solutions to Einstein's equations

k. 

Minimize the number of tensor components according to its symmetries, relabel, redefine or count the number of independent tensor components

l. 

New functionality and display for inert names and inert tensors

m. 

Automatic setting of Dirac, Paul and Gell-Mann algebras

n. 

Simplification of products of Dirac matrices

o. 

New Physics:-Library commands to perform matrix operations in expressions involving spinors with omitted indices

p. 

Miscellaneous improvements

 

• 

Physics Maple 2017 updates

q. 

General Relativity: classification of solutions to Einstein's equations and the Tetrads package

r. 

The 3D metric and the ThreePlusOne (3 + 1) new Physics subpackage

s. 

Tensors in Special and General Relativity

t. 

The StandardModel new Physics subpackage

 

• 

Physics Maple 2016 updates

u. 

Completion of the Database of Solutions to Einstein's Equations

v. 

Operatorial Algebraic Expressions Involving the Differential Operators d_[mu], D_[mu] and Nabla

w. 

Factorization of Expressions Involving Noncommutative Operators

x. 

Tensors in Special and General Relativity

y. 

Vectors Package

z. 

New Physics:-Library commands

aa. 

Redesigned Functionality and Miscellaneous

 

• 

Physics Maple 2015 updates

ab. 

Simplification

ac. 

Tensors

ad. 

Tetrads in General Relativity

ae. 

More Metrics in the Database of Solutions to Einstein's Equations

af. 

Commutators, AntiCommutators, and Dirac notation in quantum mechanics

ag. 

New Assume command and new enhanced Mode: automaticsimplification

ah. 

Vectors Package

ai. 

New Physics:-Library commands

aj. 

Miscellaneous

 

• 

Physics Maple 18 updates

ak. 

Simplification

al. 

4-Vectors, Substituting Tensors

am. 

Functional Differentiation

an. 

More Metrics in the Database of Solutions to Einstein's Equations

ao. 

Commutators, AntiCommutators

ap. 

Expand and Combine

aq. 

New Enhanced Modes in Physics Setup

ar. 

Dagger

as. 

Vectors Package

at. 

New Physics:-Library commands

au. 

Miscellaneous

 

• 

Physics Maple 17 updates

av. 

Tensors and Relativity: ExteriorDerivative, Geodesics, KillingVectors, LieDerivative, LieBracket, Antisymmetrize and Symmetrize

aw. 

Dirac matrices, commutators, anticommutators, and algebras

ax. 

Vector Analysis

ay. 

A new Library of programming commands for Physics

 

• 

Physics Maple 16 updates

az. 

Tensors in Special and General Relativity: contravariant indices and new commands for all the General Relativity tensors

ba. 

New commands for working with expressions involving anticommutative variables and functions: Gtaylor, ToFieldComponents, ToSuperfields

bb. 

Vector Analysis: geometrical coordinates with funcional dependency

Mapleprimes Physics posts

 

 

1. 

The Physics project at Maplesoft

2. 

Mini-Course: Computer Algebra for Physicists

3. 

A Complete Guide for Tensor computations using Physics

4. 

Perimeter Institute-2015, Computer Algebra in Theoretical Physics (I)

5. 

IOP-2016, Computer Algebra in Theoretical Physics (II)

6. 

ACA-2017, Computer Algebra in Theoretical Physics (III) 

 

• 

General Relativity

 

7. 

General Relativity using Computer Algebra

8. 

Exact solutions to Einstein's equations 

9. 

Classification of solutions to Einstein's equations and the ThreePlusOne (3 + 1) package 

10. 

Tetrads and Weyl scalars in canonical form 

11. 

Equivalence problem in General Relativity 

12. 

Automatic handling of collision of tensor indices in products 

13. 

Minimize the number of tensor components according to its symmetries

• 

Quantum Mechanics

 

14. 

Quantum Commutation Rules Basics 

15. 

Quantum Mechanics: Schrödinger vs Heisenberg picture 

16. 

Quantization of the Lorentz Force 

17. 

Magnetic traps in cold-atom physics 

18. 

The hidden SO(4) symmetry of the hydrogen atom

19. 

(I) Ground state of a quantum system of identical boson particles 

20. 

(II) The Gross-Pitaevskii equation and Bogoliubov spectrum 

21. 

(III) The Landau criterion for Superfluidity 

22. 

Simplification of products of Dirac matrices

23. 

Algebra of Dirac matrices with an identity matrix on the right-hand side

24. 

Factorization with non-commutative variables

25. 

Tensor Products of Quantum State Spaces 

26. 

Coherent States in Quantum Mechanics 

27. 

The Zassenhaus formula and the Pauli matrices 

 

• 

Physics package generic functionality

 

28. 

Automatic simplification and a new Assume (as in "extended assuming")

29. 

Wirtinger derivatives and multi-index summation

See Also

 

Conventions used in the Physics package , Physics , Physics examples , A Complete Guide for Tensor computations using Physics


 

Download Physics-Updates.mw
 

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

First 752 753 754 755 756 757 758 Last Page 754 of 2249