Maple 2023 Questions and Posts

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

In the end of numerics - Mathematica Stack Exchange,  claimed that "I don't know if it also has generalized precision tracking". (Another answer additionally claimed that Mma chooses "significance arithmetic" to track the precision.) So, the question is simple: Does Maple always (automatically) keep track of the accuracy of inexact numbers through a computation?

I beleave this must be wrong - plea

NULL

int(int(x^2+1, x), x)

(1/12)*(x^2+3)^2``

(1)

int(x^2+1, x)

(1/3)*x^3+x

(2)

int((1/3)*x^3+x, x)

(1/12)*(x^2+3)^2

(3)

NULL

NULL


Download Bug_in_Maple_2023.mw

se see attachment !

Any workaround for this serious error in Maple 2023? Is this new bug? The problem is can't be cought so whole program crashes. Is this related to the new changes made in 2023 to simplification of of trig functions?

restart;

711624

interface(version);

`Standard Worksheet Interface, Maple 2023.0, Windows 10, March 6 2023 Build ID 1689885`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1463. The version installed in this computer is 1462 created 2023, June 10, 2:26 hours Pacific Time, found in the directory C:\Users\Owner\maple\toolbox\2023\Physics Updates\lib\`

expr:=-1/7-(-1/7*I*7^(5/7)*exp(2/7*I*Pi)*sin(1/7*Pi)-1/7*cos(1/7*Pi)*7^(5/7)*exp(2/7*I*Pi))^(7/2);
simplify(expr);

-1/7-(-((1/7)*I)*7^(5/7)*exp(((2/7)*I)*Pi)*sin((1/7)*Pi)-(1/7)*cos((1/7)*Pi)*7^(5/7)*exp(((2/7)*I)*Pi))^(7/2)

Error, (in trig/normal/sincosargs) too many levels of recursion

 

Download simplify_error_june_13_2023.mw

Update

I just tried it in 2022, and it completes instantly with no error.

interface(version)

`Standard Worksheet Interface, Maple 2022.2, Windows 10, October 23 2022 Build ID 1657361`

expr:=-1/7-(-1/7*I*7^(5/7)*exp(2/7*I*Pi)*sin(1/7*Pi)-1/7*cos(1/7*Pi)*7^(5/7)*exp(2/7*I*Pi))^(7/2);
simplify(expr);

-1/7-(-((1/7)*I)*7^(5/7)*exp(((2/7)*I)*Pi)*sin((1/7)*Pi)-(1/7)*cos((1/7)*Pi)*7^(5/7)*exp(((2/7)*I)*Pi))^(7/2)

-1/7-(1/7)*(-(-1)^(3/7))^(1/2)*(-1)^(2/7)

 

Download simplify_ok_2022.mw

Hi community!

In the mw file attached, I have worked a numerical solution to get the standard parameters of minimum and maximum of the distribution Triangular(a,b,c) given 2 quantiles (and its values) and c.  The code works and was simple.  Now, how do I solve this more elegantly to have the reparameterization directly inputed instead of indirectly by solving a & b first (like the attached file does)? Thanks people,altTriangular2.mw

From help it says

But in the following it works on list type. Why? Should have not given an error?

158332

restart;

158332

interface(warnlevel=4);
kernelopts('assertlevel'=2):

3

T:=[ [1,2],[3,4],[6,6]];
whattype(T);
type(T,Matrix);
LinearAlgebra:-RowDimension(T);

[[1, 2], [3, 4], [6, 6]]

list

false

3

T:=Matrix([[1,2],[3,4],[6,6]]);
whattype(T);
type(T,Matrix);
LinearAlgebra:-RowDimension(T);

T := Matrix(3, 2, {(1, 1) = 1, (1, 2) = 2, (2, 1) = 3, (2, 2) = 4, (3, 1) = 6, (3, 2) = 6})

Matrix

true

3

 

Download row_dim_june_13_2023.mw

Where do you think label="dontexpand" is coming from in the following? I never seen this before. Is this a bug?

When using PDEtools:-Solve no such error shows up. Worksheet below

998948

interface(version);

`Standard Worksheet Interface, Maple 2023.0, Windows 10, March 6 2023 Build ID 1689885`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1462 and is the same as the version installed in this computer, created 2023, June 10, 2:26 hours Pacific Time.`

restart;

998948

sol:=sum((-2*_R^2+1)/(6*_R^2+4*_R-1)*ln(y(x)/x-_R),_R = RootOf(2*_Z^3+2*_Z^2-_Z+1))-ln(x)-c[1] = 0;
solve(sol,y(x))

sum((-2*_R^2+1)*ln(y(x)/x-_R)/(6*_R^2+4*_R-1), _R = RootOf(2*_Z^3+2*_Z^2-_Z+1))-ln(x)-c[1] = 0

Error, (in solve) cannot solve expressions with sum((-2*_R^2+1)/(6*_R^2+4*_R-1)*ln(y(x)/x-_R),_R = RootOf(2*_Z^3+2*_Z^2-_Z+1,label = "dontexpand")) for y(x)

PDEtools:-Solve(sol,y(x))

y(x) = exp(-c[1])+_R*x

 

Download dontexpand_june_10_2023.mw

Maple has a myriad of kernel functions for doing different kinds of symbolic replacements to whole expressions: subs, eval, algsubs, applyrule, `simplify/siderels`, `simpl/eval`,  Physics:-Substitute, MTM:-subs, MmaTranslator:-Mma:-ReplaceRepeated, PDEtools:-dsubs, liesymm:-wsubs, student:-powsubs, etc. But if I need to apply transformation rules over and over again until the result no longer changes within  iterations in a singular clean built-in command (so, without explicit while / until / MmaTranslator:-Mma:-FixedPoint) elegantly, only four can be called: eval['recurse']algsubs, applyrule, MmaTranslator['Mma']['ReplaceRepeated']. The difference between applyrule and algsubs has been elucidated in this help page:

applyrule … does not do mathematical transformations as algsubs does.

However, I cannot find any explanation for the potential distinctions between eval['recurse'] and MmaTranslator['Mma']['ReplaceRepeated'] in their documentation. Does anyone know?

Besides, is it possible to use an "operator" form (like ) instead of a "functional" form (like ) to perform (single or repeated) substitutions in Maple?

How can one set their own symbols for DynamicSystems continuoustimevar?  If I try to use DynamicSystems:-SystemOptions('continuoustimevar'=x): or DynamicSystems:-SystemOptions('continuoustimevar'=y): Maple complains that these letters are already assigned. Why can not one have their own choice of which letters to use?  Is there a workaround?


 

858904

interface(version);

`Standard Worksheet Interface, Maple 2023.0, Windows 10, March 6 2023 Build ID 1689885`

restart;

858904

DynamicSystems:-SystemOptions('continuoustimevar'=t):

restart;

858904

DynamicSystems:-SystemOptions('continuoustimevar'=x):

Error, (in DynamicSystems:-SystemOptions) cannot assign x to continuoustimevar, already assigned to statevariable

restart;

858904

DynamicSystems:-SystemOptions('continuoustimevar'=y):

Error, (in DynamicSystems:-SystemOptions) cannot assign y to continuoustimevar, already assigned to outputvariable

restart;

858904

DynamicSystems:-SystemOptions('continuoustimevar'=XXXXXX):

 


 

Download dynamics_june_9_2023.mw

 

Can one evaluate dirac spinor products using the Standard Model package? How far can I take the evaluation in the package? See maplesheet with commentary.

Evaluating_spinor_products_with_the_standard_model_package.mw

When adding singsol=all to this ode in the first example, dsolve returns nothing. But it should have returned the general solution if it can not find singular solution. This is what it does always as can be seen from the second example given below.

Why did dsolve return nothing in the first example? Is this a bug?

319824

interface(version);

`Standard Worksheet Interface, Maple 2023.0, Windows 10, March 6 2023 Build ID 1689885`

restart;

319824

ode:=y(x)=x*diff(y(x),x)+ a*diff(y(x),x)/(sqrt(1+diff(y(x),x)^2));
sol_1:=dsolve(ode,y(x));
sol_2:=dsolve(ode,y(x),singsol=all);

y(x) = x*(diff(y(x), x))+a*(diff(y(x), x))/(1+(diff(y(x), x))^2)^(1/2)

y(x) = x*c__1+a*c__1/(c__1^2+1)^(1/2)

"sol_2 := "

ode:=diff(y(x),x)^2-(1+2*x*y(x))*diff(y(x),x)+2*x*y(x) = 0;
sol_1:=dsolve(ode,y(x));
sol_2:=dsolve(ode,y(x), singsol=all)

(diff(y(x), x))^2-(1+2*y(x)*x)*(diff(y(x), x))+2*y(x)*x = 0

y(x) = c__1*exp(x^2), y(x) = x+c__1

y(x) = c__1*exp(x^2), y(x) = x+c__1

 

Download missing_sol.mw

A chordal graph is one in which all cycles of four or more vertices have a chord, which is an edge that is not part of the cycle but connects two vertices of the cycle. A perfect elimination ordering in a graph is an ordering of the vertices of the graph such that, for each vertex v, v and the neighbors of v that occur after v in the order form a clique.  A graph is chordal if and only if it has a perfect elimination ordering.

I  use IsChordal  to test whether the lexicographic product of  two graphs g1,g2 is a chordal graph. It returned true and provided a perfect elimination sequence 1, 2, ..., 30. However, vertices of s  are "1:1", "1:2", ..., "10:3", rather than using Arabic numerals. Therefore, it is difficult for me to extract useful information from the perfect elimination sequence.

with(GraphTheory):
g1:=PathGraph(10):
g2:=CycleGraph(3):

s:=LexicographicProduct(g1,g2):

IsChordal(s,eliminationordering=true)

true, [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]

(1)

DrawGraph(s)

 

Vertices(s)

["1:1", "1:2", "1:3", "2:1", "2:2", "2:3", "3:1", "3:2", "3:3", "4:1", "4:2", "4:3", "5:1", "5:2", "5:3", "6:1", "6:2", "6:3", "7:1", "7:2", "7:3", "8:1", "8:2", "8:3", "9:1", "9:2", "9:3", "10:1", "10:2", "10:3"]

(2)

Download ischordgraph.mw

Hello community, I am new in this forum and sorry if the following seems rudimentary:

I am replicating a distribution function using Statistics[Distribution] as defined here:

Vose Software's Beta Subjective

I have tested numerically and the function works for this initial values: Min := 3; Mlikely := 8; Avg := 9; Max := 18;

But failed to complete with these values (which I tested works in another software) Min := 1000; Mlikely := 1400; Avg := 1500; Max := 2100;

What puzzles me is the inconsistency of handling the functions (be aware that the server kernel could be slow at times). I hope I don't have to give up and continue my custom distribution project in another software.

Here is the document: 

NULL

The Beta Subjective Distribution

 

 

restart

 

with(Statistics)

NULL

Min := 3; Mlikely := 8; Avg := 9; Max := 18

18

(1)

NULL

Mid := (Min+Max)*(1/2)

21/2

(2)

alpha := 2*(Avg-Min)*(Mid-Mlikely)/((Avg-Mlikely)*(Max-Min))

2

(3)

beta := alpha*(Max-Avg)/(Avg-Min)

3

(4)

NULLNULLNULLNULL

NULL

f := simplify(piecewise(Min <= x and x <= Max, (x-Min)^(alpha-1)*(Max-x)^(beta-1)/(Beta(alpha, beta)*(Max-Min)^(alpha+beta-1)), 0))

piecewise(x < 3, 0, x <= 18, (4/16875)*(x-3)*(-18+x)^2, 18 < x, 0)

(5)

NULL

NULL

MD := Distribution(PDF = unapply(f, x), Conditions = [`and`(Min < Mlikely and Mlikely < Max and Min < Avg and Avg < Max and Mlikely <> Avg, piecewise(Avg < Mlikely, Mlikely > Mid, Mid > Mlikely))])

_m2352873740384

(6)

NULL

NULL

NULL

NULL

NULL

X := RandomVariable(MD)

_R

(7)

NULL

evalf(Mean(X))NULL

9.

(8)

NULL

NULL

evalf(Mode(X))

{8.}

(9)

evalf(Quantile(X, .25))

HFloat(6.645331256340534)

(10)

NULL

A := Sample(X, 10^5)

Mean(A)

HFloat(8.989513601786443)

(11)

Mode(A)

HFloat(8.399525182283146)

(12)

Variance(A)

HFloat(9.029923294780689)

(13)

Skewness(A)

HFloat(0.29817612117709397)

(14)

Kurtosis(A)

HFloat(2.3706624578773154)

(15)

Quantile(A, .25)

HFloat(6.63427086578336)

(16)

NULL

Download Distribution-Beta-Subjective.mw

Hello there, 

This is another issue, associated with substitution. In the following Maple expressions, I tried to substitutte the denominator of 'eq_K1_m4' in order to make it as 'eq_K1_m4_desired', but did not get any success (yet). 

Therefore, would you have a look at this issue to see if the intended goal can be achieved?

restart;

with(LinearAlgebra):

with(DynamicSystems):

interface(imaginaryunit=j):

eq_K1_m4 := K__1 = E__q0*(R__T*E__B*sin(delta) + X__Td*E__B*cos(delta))/(L__aqs*L__l + L__aqs*X__E + L__aqs*L__ads_p + L__l^2 + 2*L__l*X__E + L__l*L__ads_p + R__E^2 + 2*R__E*R__a + R__a^2 + X__E^2 + X__E*L__ads_p) + (X__q - X__dp)*i__q0*(X__Tq*E__B*sin(delta) - R__T*E__B*cos(delta))/(L__aqs*L__l + L__aqs*X__E + L__aqs*L__ads_p + L__l^2 + 2*L__l*X__E + L__l*L__ads_p + R__E^2 + 2*R__E*R__a + R__a^2 + X__E^2 + X__E*L__ads_p);

K__1 = E__q0*(R__T*E__B*sin(delta)+X__Td*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)+(X__q-X__dp)*i__q0*(X__Tq*E__B*sin(delta)-R__T*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)

(1)

eq_K1_m4_desired := K__1 = E__q0*(R__T*E__B*sin(delta) + X__Td*E__B*cos(delta))/Dx + (X__q - X__dp)*i__q0*(X__Tq*E__B*sin(delta) - R__T*E__B*cos(delta))/Dx;

K__1 = E__q0*(R__T*E__B*sin(delta)+X__Td*E__B*cos(delta))/Dx+(X__q-X__dp)*i__q0*(X__Tq*E__B*sin(delta)-R__T*E__B*cos(delta))/Dx

(2)

eq_Dx := Dx = L__aqs*L__l + L__aqs*X__E + L__aqs*L__ads_p + L__l^2 + 2*L__l*X__E + L__l*L__ads_p + R__E^2 + 2*R__E*R__a + R__a^2 + X__E^2 + X__E*L__ads_p;

Dx = L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p

(3)

denom(op(1, rhs(eq_K1_m4))) - rhs(eq_Dx); # checking to see if the denominator expression is the same as the expression of Dx

0

(4)

denom(op(2, rhs(eq_K1_m4))) - rhs(eq_Dx); # checking to see if the denominator expression is the same as the expression of Dx

0

(5)

# 1

map2(applyrule, eq_Dx, eq_K1_m4); # did not work.

K__1 = E__q0*(R__T*E__B*sin(delta)+X__Td*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)+(X__q-X__dp)*i__q0*(X__Tq*E__B*sin(delta)-R__T*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)

(6)

# 2

subs(eq_Dx, eq_K1_m4); # did not work.

K__1 = E__q0*(R__T*E__B*sin(delta)+X__Td*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)+(X__q-X__dp)*i__q0*(X__Tq*E__B*sin(delta)-R__T*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)

(7)

# 3

simplify(eq_K1_m4, {Dx = L__aqs*L__l + L__aqs*X__E + L__aqs*L__ads_p + L__l^2 + 2*L__l*X__E + L__l*L__ads_p + R__E^2 + 2*R__E*R__a + R__a^2 + X__E^2 + X__E*L__ads_p}, [Dx]); # did not work.

K__1 = -E__B*((-R__T*(-X__q+X__dp)*i__q0-E__q0*X__Td)*cos(delta)+sin(delta)*(X__Tq*(-X__q+X__dp)*i__q0-R__T*E__q0))/(L__l^2+(L__aqs+L__ads_p+2*X__E)*L__l+X__E^2+(L__aqs+L__ads_p)*X__E+L__aqs*L__ads_p+(R__E+R__a)^2)

(8)

# 4

algsubs(eq_Dx, eq_K1_m4); # did not work.

K__1 = E__q0*(R__T*E__B*sin(delta)+X__Td*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)+(X__q-X__dp)*i__q0*(X__Tq*E__B*sin(delta)-R__T*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)

(9)

# 5

applyrule(L__aqs*L__l + L__aqs*X__E + L__aqs*L__ads_p + L__l^2 + 2*L__l*X__E + L__l*L__ads_p + R__E^2 + 2*R__E*R__a + R__a^2 + X__E^2 + X__E*L__ads_p = Dx, eq_K1_m4); # did not work.

K__1 = E__q0*(R__T*E__B*sin(delta)+X__Td*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)+(X__q-X__dp)*i__q0*(X__Tq*E__B*sin(delta)-R__T*E__B*cos(delta))/(L__aqs*L__l+L__aqs*X__E+L__aqs*L__ads_p+L__l^2+2*L__l*X__E+L__l*L__ads_p+R__E^2+2*R__E*R__a+R__a^2+X__E^2+X__E*L__ads_p)

(10)

 

Download Q20230606.mw

I am surprised Maple could not do this integral. Any ideas how to make it do it?

restart;
integrand:=(ln(-x + sqrt(x^2 - 1))*x + sqrt(x^2 - 1))*x*sqrt(-x^2 + 1)/sqrt(x^2 - 1);
int(integrand,x)

For reference, another software gives

ClearAll[x]
integrand = (Log[-x + Sqrt[x^2 - 1]]*x + Sqrt[x^2 - 1])*x*Sqrt[-x^2 + 1]/Sqrt[x^2 - 1]
Integrate[integrand, x]

Which I verified to be correct.

Any workarounds? I tried assuming x>0 for example. Not sure what else to try now. I am sure with some pre-processing Maple can do it. 

Maple 2023

I was unsure if this should be a post or a question, or how to word it properly as I feel the attached file is the best way to see my issue. 

 However, when doing a computation involing the D_ specifically D_[mu](F[~mu,~nu]) where is a tensor, Maple is putting in strange "variables" into the output of my expression. I set the spacetime variables myself to tau,x,y,z and in the output there is combinations which appear as x(tau),y(tau) which has completely baffled me. 

Any ideas/help would be greatly appreciated. 

EinsteinPert.mw

EDIT: The problem is resolved when the Physics[Vectors] package is removed and replaced by VectorCalculus for the Lapacian. Still curious as why the Physics[Vectors] package causes this issue. (I of course could just brute force the Laplacian but I thought I would use the power of Maple for this).

First 21 22 23 24 25 26 27 Last Page 23 of 32