Advanced Search

Product Suggestions

Post your suggestions on new features and products.
Create a new Post in this category

Integrate Heaviside bugs

by: vv 13620 Maple
integral heaviside
November 15 2016
3  7

Correct computatiton for

for reasonable expressions f(x,y), g(x,y) would be very useful in double integrals.

For the moment this is not possible. Too many bugs:

int(Heaviside(1-x^2-y^2), x=-infinity..infinity, y=-infinity..infinity); #should be Pi
                           undefined
int(Heaviside(1-x^2-y^2), x=-1..1, y=-1..1); #should be Pi
                               0
int(Heaviside(y-x^2), x=-1..1, y=-1..1); #should be 4/3
                               -2

int(Heaviside(y-x^2), y=-1..1, x=-1..1); #This one is OK!
                              4/3

 

 

 

 

Cannot submit SCRs

by: Markiyan Hirnyk 7809 Maple , MaplePrimes
August 29 2016
3  1

In order to change Maple for the better, I use to submit SCRs. However, as i was kindly
informed by Bryon (a copy of his e-letter on demand), MaplePrimes are under reconstruction and do not
work properly. At least my three messages sent through the Contact button were lost.
I have  unsuccessfully tried to reach beta.maplesoft.com (see the result of ping in the screen screen.29.08.16.docx).
Please, help me!

Defining a quantum operator as a tensor – And...

by: Pascal4QM 478 Maple 2016
tensor commutator quantum
August 24 2016
2  2

Hi,

In the following example I introduce some commutation rules that are standard in Quantum Mechanics. A major feature of the Maple Physics Package, is that it is possible to define tensors as Quantum Operators. This is of great interest because powerful tensor simplification rules can then be used in Quantum Mechanics. For an example, see the commutation rules of the components of the angular momentum operator in ?Physics,Examples. Here, I focus on a possible issue: when destroying all quantum operators, the pre-defined commutation rules still apply, which should not be the case. As shown in the post, this is link to the fact that these operators are also tensors.
 

NULL

 

> 

Physics:-Version()[2]

`2016, August 16, 18:56 hours`

 (1)

NULL

NULL

> 

restart; with(Physics); interface(imaginaryunit = I)

First, set a 3D Euclidian space

> 

Setup(mathematicalnotation = true, dimension = 3, signature = `+`, spacetimeindices = lowercaselatin, quiet)

[dimension = 3, mathematicalnotation = true, signature = `+ + +`, spacetimeindices = lowercaselatin]

 (2)

Define two rank 1 tensors

> 

Define(x[k], p[k])

`Defined objects with tensor properties`

  

{Physics:-Dgamma[a], Physics:-Psigma[a], Physics:-d_[a], Physics:-g_[a, b], p[k], x[k], Physics:-KroneckerDelta[a, b], Physics:-LeviCivita[a, b, c]}

 (3)

Now, further define these tensors as quantum operators and gives the usual commutation rule between position and momentum operators (Quantum Mechanics).

> 

Setup(hermitianoperators = {p, x}, algebrarules = {%Commutator(p[k], p[n]) = 0, %Commutator(x[k], p[l]) = I*`ℏ`*KroneckerDelta[k, l], %Commutator(x[k], x[l]) = 0}, realobjects = {`ℏ`})

[algebrarules = {%Commutator(p[k], p[n]) = 0, %Commutator(x[k], p[l]) = I*`ℏ`*Physics:-KroneckerDelta[k, l], %Commutator(x[k], x[l]) = 0}, hermitianoperators = {p, x}, realobjects = {`ℏ`}]

 (4)

As expected:

> 

(%Commutator = Commutator)(p[a], x[b])

%Commutator(p[a], x[b]) = -I*`ℏ`*Physics:-KroneckerDelta[a, b]

 (5)

Now, reset all the Hermitian operators, so that all quantum operators are destroyed. This is useful if, for instance, one needs to compare some the result with the commutative case.

> 

Setup(redo, hermitianoperators = {})

[hermitianoperators = none]

 (6)

As expected, there are no quantum operators anymore...

> 

Setup(quantumoperators)

[quantumoperators = {}]

 (7)

...so that the following expressions should commute (result should be true)

> 

Library:-Commute(p[a], x[b])

false

 (8)

Result should be 0NULL

> 

Commutator(p[a], x[b])

-I*`ℏ`*Physics:-KroneckerDelta[a, b]

 (9)
> 

p[a], x[b]

p[a], x[b]

 (10)
> 

NULL

NULL

``

NULLNULL

Below is just a copy & paste of the above section. The only difference, is that "Define(x[k], p[k])" has been commented, so that x[k]and p[k] are not a tensor. In that case, everything behaves as expected (but of course, the interesting feature of tensors is not available).

````

NULL

> 

restart; with(Physics); interface(imaginaryunit = I)

First, set a 3D Euclidian space

> 

Physics:-Setup(mathematicalnotation = true, dimension = 3, signature = `+`, spacetimeindices = lowercaselatin, quiet)

[dimension = 3, mathematicalnotation = true, signature = `+ + +`, spacetimeindices = lowercaselatin]

 (11)

#Define two rank 1 tensors

> 

Now, further define these tensors as quantum operators and gives the usual commutation rule between position and momentum operators (Quantum Mechanics)

> 

Physics:-Setup(hermitianoperators = {p, x}, algebrarules = {%Commutator(p[k], p[n]) = 0, %Commutator(x[k], p[l]) = Physics:-`*`(Physics:-`*`(I, `ℏ`), Physics:-KroneckerDelta[k, l]), %Commutator(x[k], x[l]) = 0}, realobjects = {`ℏ`})

[algebrarules = {%Commutator(p[k], p[n]) = 0, %Commutator(x[k], p[l]) = I*`ℏ`*Physics:-KroneckerDelta[k, l], %Commutator(x[k], x[l]) = 0}, hermitianoperators = {p, x}, realobjects = {`ℏ`}]

 (12)

As expected:

> 

(%Commutator = Physics:-Commutator)(p[a], x[b])

%Commutator(p[a], x[b]) = -I*`ℏ`*Physics:-KroneckerDelta[a, b]

 (13)

Now, reset all the Hermitian operators, so that all quantum operators are destroyed.

> 

Physics:-Setup(redo, hermitianoperators = {})

[hermitianoperators = none]

 (14)

As expected, there are no quantum operators anymore...

> 

Physics:-Setup(quantumoperators)

[quantumoperators = {}]

 (15)

...so that the following expressions should commute (result should be true)

> 

Physics:-Library:-Commute(p[a], x[b])

true

 (16)

Result should be 0``

> 

Physics:-Commutator(p[a], x[b])

0

 (17)
> 

p[a], x[b]

p[a], x[b]

 (18)
> 

NULL

``

NULL``

NULL


Download Quantum_operator_as_Tensors_August_23_2016.mw

odd Zeta value series approximation to the value of Ei...

by: Adam Ledger 360 Maple
constants exponential
August 18 2016
0  0


> 

(41403013072+527237505*Zeta(3)^2+8631147672*Zeta(3)+527237505*Zeta(5)^2+sqrt(1714209491440202877184+118155152758998350304*Zeta(3)^2+714711039780356736768*Zeta(3)+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(3)^4+9101329527743676720*Zeta(3)^3+555958773357250050*Zeta(3)^2*Zeta(5)^2+9101329527743676720*Zeta(3)*Zeta(5)^2+277979386678625025*Zeta(5)^4))*(41403013072+527237505*Zeta(5)^2-sqrt(1714209491440202877184+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(5)^4))/((41403013072+527237505*Zeta(3)^2+8631147672*Zeta(3)+527237505*Zeta(5)^2-sqrt(1714209491440202877184+118155152758998350304*Zeta(3)^2+714711039780356736768*Zeta(3)+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(3)^4+9101329527743676720*Zeta(3)^3+555958773357250050*Zeta(3)^2*Zeta(5)^2+9101329527743676720*Zeta(3)*Zeta(5)^2+277979386678625025*Zeta(5)^4))*(41403013072+527237505*Zeta(5)^2+sqrt(1714209491440202877184+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(5)^4)))-3/7-(1/10)*(exp(Pi)+5)/(Ei(1)+1/2)

(41403013072+527237505*Zeta(3)^2+8631147672*Zeta(3)+527237505*Zeta(5)^2+(1714209491440202877184+118155152758998350304*Zeta(3)^2+714711039780356736768*Zeta(3)+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(3)^4+9101329527743676720*Zeta(3)^3+555958773357250050*Zeta(3)^2*Zeta(5)^2+9101329527743676720*Zeta(3)*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2))*(41403013072+527237505*Zeta(5)^2-(1714209491440202877184+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2))/((41403013072+527237505*Zeta(3)^2+8631147672*Zeta(3)+527237505*Zeta(5)^2-(1714209491440202877184+118155152758998350304*Zeta(3)^2+714711039780356736768*Zeta(3)+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(3)^4+9101329527743676720*Zeta(3)^3+555958773357250050*Zeta(3)^2*Zeta(5)^2+9101329527743676720*Zeta(3)*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2))*(41403013072+527237505*Zeta(5)^2+(1714209491440202877184+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2)))-3/7-(1/10)*(exp(Pi)+5)/(Ei(1)+1/2)

 (1)
> 

abs(evalf[1000]((41403013072+527237505*Zeta(3)^2+8631147672*Zeta(3)+527237505*Zeta(5)^2+(1714209491440202877184+118155152758998350304*Zeta(3)^2+714711039780356736768*Zeta(3)+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(3)^4+9101329527743676720*Zeta(3)^3+555958773357250050*Zeta(3)^2*Zeta(5)^2+9101329527743676720*Zeta(3)*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2))*(41403013072+527237505*Zeta(5)^2-(1714209491440202877184+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2))/((41403013072+527237505*Zeta(3)^2+8631147672*Zeta(3)+527237505*Zeta(5)^2-(1714209491440202877184+118155152758998350304*Zeta(3)^2+714711039780356736768*Zeta(3)+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(3)^4+9101329527743676720*Zeta(3)^3+555958773357250050*Zeta(3)^2*Zeta(5)^2+9101329527743676720*Zeta(3)*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2))*(41403013072+527237505*Zeta(5)^2+(1714209491440202877184+30838267512743688864*Zeta(5)^2+277979386678625025*Zeta(5)^4)^(1/2)))-3/7-(1/10)*(exp(Pi)+5)/(Ei(1)+1/2)))

0.4022755651e-7

 (2)
> 

``


Download EXPONENTIAL_INTEGRAL_ODD_ZETA_SERIES_APPROXIMATION.mw

CONGRUENT FUNCTIONS OF THE FRACTIONAL PART OVER...

by: Adam Ledger 360 Maple
August 17 2016
0  1

A more honest and specific version of lemma 3.

CONGRUENT_FUNCTIONS_OF_THE_FRACTIONAL_PART_OVER_Q_LEMMA_4.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/CONGRUENT_FUNCTIONS_OF_THE_FRACTIONAL_PART_OVER_Q_LEMMA_4.mw .

Download CONGRUENT_FUNCTIONS_OF_THE_FRACTIONAL_PART_OVER_Q_LEMMA_4.mw

Help Menu Bugs

by: Samuel Hollis 100 Maple 2016
August 15 2016
2  1

Using a mouse's back/forward buttons or alt+ left arrow/ alt+right arrow did not work in Maple 2015 and they still don't work in Maple 2016. 

Also, if you select anywhere in the search bar, it automatically selects all the words.  I sometimes (most the time) want to change only one word.  

Im using x64 Windows 10 and it also did it when running Windows 8/8.1 .   Does anybody else have this problem?

Can it be fixed please if not?

Make Maple touch friendly!

by: Samuel Hollis 100 Maple 2016 , MapleSim 2016
August 13 2016
1  2

If you try and scroll a page up or down, it just selects the content.  How hard would it be to fix this?

Primitive words in StringTools

by: JamieSimpson 5 Maple
stringtools
July 11 2016
1  1

It would be nice if the StringTools package had a function giving the number of primitive words of length n on an alphabet of size a.  The formula for this is well-known but tedious to code.

Linkage mechanisms

by: one man 1335 Maple 15
application animation
July 10 2016
1  5

        General description of the method of solving underdetermined systems of equations. As a particular application of the idea proposed a universal method  kinematic analysis for all kinds of  spatial and planar link mechanisms with any number degrees of freedom.  The method can be used for powerful CAD linkages.   
         http://www.maplesoft.com/applications/view.aspx?SID=154228

       


      Some examples of a much larger number calculated by the proposed method. Examples gathered here not to look for them on the forum and opportunity to demonstrate the method.  Among the examples, I think, there are very complicated.

https://vk.com/doc242471809_408704758
https://vk.com/doc242471809_408704572
https://vk.com/doc242471809_376439263
https://vk.com/doc242471809_402619761
https://vk.com/doc242471809_402610228
https://vk.com/doc242471809_401188803
https://vk.com/doc242471809_400465891
https://vk.com/doc242471809_400711315
https://vk.com/doc242471809_387358164
https://vk.com/doc242471809_380837279
https://vk.com/doc242471809_379935473
https://vk.com/doc242471809_380217387
https://vk.com/doc242471809_363266817
https://vk.com/doc242471809_353980472
https://vk.com/doc242471809_375452868
https://vk.com/doc242471809_353988163 
https://vk.com/doc242471809_353986884 
https://vk.com/doc242471809_353987119
https://vk.com/doc242471809_324249241
https://vk.com/doc242471809_324102889
https://vk.com/doc242471809_322219275
https://vk.com/doc242471809_437298137
https://vk.com/doc242471809_437308238
https://vk.com/doc242471809_437308241
https://vk.com/doc242471809_437308243
https://vk.com/doc242471809_437308245
https://vk.com/doc242471809_437308246
https://vk.com/doc242471809_437401651
https://vk.com/doc242471809_437664558

 

 

Error and nonsense in eulermac outputs

by: vv 13620 Maple
euler-maclaurin
March 16 2016
1  3

eulermac(1/(n*ln(n)^2),n=2..N,1);  #Error
Error, (in SumTools:-DefiniteSum:-ClosedForm) summand is singular in the interval of summation


eulermac(1/(n*ln(n)^2+1),n=2..N,1);  #nonsense

 

 

Calculating linkage mechanisms

by: one man 1335 Maple 17
animation
March 04 2016
3  12

The method of solving underdetermined systems of equations, and universal method for calculating link mechanisms. It is based on the Draghilev’s method for solving systems of nonlinear equations. 
When calculating link mechanisms we can use geometrical relationships to produce their mathematical models without specifying the “input link”. The new method allows us to specify the “input link”, any link of mechanism.

Example.
Three-bar mechanism.  The system of equations linkages in this mechanism is as follows:

f1 := x1^2+(x2+1)^2+(x3-.5)^2-R^2;
f2 := x1-.5*x2+.5*x3;
f3 := (x1-x4)^2+(x2-x5)^2+(x3-x6)^2-19;
f4 := sin(x4)-x5;
f5 := sin(2*x4)-x6;

Coordinates green point x'i', i = 1..3, the coordinates of red point x'i', i = 4..6.
Set of x0'i', i = 1..6 searched arbitrarily, is the solution of the system of equations and is the initial point for the solution of the ODE system. The solution of ODE system is the solution of system of equations linkages for concrete assembly linkage.
Two texts of the program for one mechanism. In one case, the “input link” is the red-green, other case the “input link” is the green-blue.
After the calculation trajectories of points, we can always find the values of other variables, for example, the angles.
Animation displays the kinematics of the mechanism.
MECAN_3_GR_P_bar.mw 
MECAN_3_Red_P_bar.mw

(if to use another color instead of color = "Niagara Dark Orchid", the version of Maple <17)

Method_Mechan_PDF.pdf






Calculating linkage mechanisms

by: one man 1335 Maple 17
animation
March 04 2016
3  12

The method of solving underdetermined systems of equations, and universal method for calculating link mechanisms. It is based on the Draghilev’s method for solving systems of nonlinear equations. 
When calculating link mechanisms we can use geometrical relationships to produce their mathematical models without specifying the “input link”. The new method allows us to specify the “input link”, any link of mechanism.

Example.
Three-bar mechanism.  The system of equations linkages in this mechanism is as follows:

f1 := x1^2+(x2+1)^2+(x3-.5)^2-R^2;
f2 := x1-.5*x2+.5*x3;
f3 := (x1-x4)^2+(x2-x5)^2+(x3-x6)^2-19;
f4 := sin(x4)-x5;
f5 := sin(2*x4)-x6;

Coordinates green point x'i', i = 1..3, the coordinates of red point x'i', i = 4..6.
Set of x0'i', i = 1..6 searched arbitrarily, is the solution of the system of equations and is the initial point for the solution of the ODE system. The solution of ODE system is the solution of system of equations linkages for concrete assembly linkage.
Two texts of the program for one mechanism. In one case, the “input link” is the red-green, other case the “input link” is the green-blue.
After the calculation trajectories of points, we can always find the values of other variables, for example, the angles.
Animation displays the kinematics of the mechanism.
MECAN_3_GR_P_bar.mw 
MECAN_3_Red_P_bar.mw

(if to use another color instead of color = "Niagara Dark Orchid", the version of Maple <17)

Method_Mechan_PDF.pdf






Apply to be a Maple Beta Tester!

by: Daniel Skoog 1771 Maple
beta-testing
December 22 2015
4  2

As the year draws to a close, we start looking forward to a new year and a new release of Maple. With every new release comes many new features and updates to explore.

We are looking for several new beta testers with a good working knowledge of Maple; We need your input, your ideas, and your experience with our products to help us improve the software and get it ready for general release.

There are many benefits to becoming a beta tester:

  • You’ll get to use the new software before anyone else does.
  • You’ll help us make our software better in ways that work for you.
  • Your suggestions could determine the future direction of the software.
  • You’ll get feedback right from the development team.

If you are interested in becoming a beta tester for the next version of Maple, please email: beta (at) maplesoft.com for more information.

Apply to be a Maple Beta Tester!

by: Daniel Skoog 1771 Maple
beta-testing
December 22 2015
4  2

As the year draws to a close, we start looking forward to a new year and a new release of Maple. With every new release comes many new features and updates to explore.

We are looking for several new beta testers with a good working knowledge of Maple; We need your input, your ideas, and your experience with our products to help us improve the software and get it ready for general release.

There are many benefits to becoming a beta tester:

  • You’ll get to use the new software before anyone else does.
  • You’ll help us make our software better in ways that work for you.
  • Your suggestions could determine the future direction of the software.
  • You’ll get feedback right from the development team.

If you are interested in becoming a beta tester for the next version of Maple, please email: beta (at) maplesoft.com for more information.

New features in Mathematica 10.3.0

by: Markiyan Hirnyk 7809 Maple
mathematica residue
October 16 2015
2  16

Mathematica 10.3.0 was announced yesterday. This is the 6th release of Mathematica 10 during 16 months. I wonder its  MathematicaFunctionData and   FindFormula . At first sight, the former is an analog of FunctionAdvisor of Maple, but the latter isn't any analog. Also compare the outputs of

Residue[Binomial[n,k],{n,-j}]

(-1)^j/(j!*k!*(-j-k)!)

 and

>`assuming`([residue(binomial(n, k), n = -j)], [integer, j > 0]);

                residue(binomial(n, k), n = -j)
Let us wait for Maple 2016.

 

4 5 6 7 8 9 10 Last Page 6 of 24