MaplePrimes Questions

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

I'm struggling a bit to understand the results of the following program.

Apparently variables are not being changed after the first run. Otherwise I can't explain the results of the last line.

GeometrySegment.mw

How to find Z(t) matrix(array) =[ [z[1,0](t),z[1,1](t)],[z[2,0](t),z[2,1](t)]]

restart;

with(IntegrationTools):with(Physics):

 

 

 

h1 := (m,n)->(n+1/2)*KroneckerDelta[n,m];

proc (m, n) options operator, arrow; Physics:-`*`(n+Physics:-`*`(1, Physics:-`^`(2, -1)), Physics:-KroneckerDelta[n, m]) end proc

(1)

h2 := (mu,nu,m2,l)->-(nu*Pi/l)^2/(2*m2)*KroneckerDelta[mu,nu];

proc (mu, nu, m2, l) options operator, arrow; Physics:-`*`(-1, Physics:-`*`(Physics:-`^`(Physics:-`*`(nu, Pi, Physics:-`^`(l, -1)), 2), Physics:-`^`(Physics:-`*`(2, m2), -1), Physics:-KroneckerDelta[mu, nu])) end proc

(2)

v1 := (m,n,m1)->sqrt(min(n,m)!/max(n,m)!)*(2*m1)^(-abs(n-m)/2)*exp(-1/(4*m1))*LaguerreL(min(n,m),abs(n-m),-1/(2*m1));

proc (m, n, m1) options operator, arrow; Physics:-`*`(sqrt(Physics:-`*`(factorial(min(n, m)), Physics:-`^`(factorial(max(n, m)), -1))), Physics:-`^`(Physics:-`*`(2, m1), Physics:-`*`(-1, Physics:-`*`(abs(n-m), Physics:-`^`(2, -1)))), exp(Physics:-`*`(-1, Physics:-`*`(1, Physics:-`^`(Physics:-`*`(4, m1), -1)))), LaguerreL(min(n, m), abs(n-m), Physics:-`*`(-1, Physics:-`*`(1, Physics:-`^`(Physics:-`*`(2, m1), -1))))) end proc

(3)

v2 := (mu,nu,l)->4*Pi^2*l*mu*nu*(exp(l/2)-(-1)^(mu+nu)*exp(-l/2))/((Pi*(mu+nu))^2+l^2)/((Pi*(mu-nu))^2+l^2);

proc (mu, nu, l) options operator, arrow; Physics:-`*`(4, Physics:-`^`(Pi, 2), l, mu, nu, exp(Physics:-`*`(l, Physics:-`^`(2, -1)))-Physics:-`*`(Physics:-`^`(-1, mu+nu), exp(Physics:-`*`(-1, Physics:-`*`(l, Physics:-`^`(2, -1))))), Physics:-`^`(Physics:-`^`(Physics:-`*`(Pi, mu+nu), 2)+Physics:-`^`(l, 2), -1), Physics:-`^`(Physics:-`^`(Physics:-`*`(Pi, mu-nu), 2)+Physics:-`^`(l, 2), -1)) end proc

(4)

h:=(m,n,mu,nu,m1,m2,l)->evalf(h1(m,n)+h2(mu,nu,m2,l)+v1(m,n,m1)+v2(mu,nu,l));

proc (m, n, mu, nu, m1, m2, l) options operator, arrow; evalf(h1(m, n)+h2(mu, nu, m2, l)+v1(m, n, m1)+v2(mu, nu, l)) end proc

(5)

m1:=1:m2:=1:l:=1:

H:= (m,n,mu,nu)->h(m,n,mu,nu,m1,m2,l);

proc (m, n, mu, nu) options operator, arrow; h(m, n, mu, nu, m1, m2, l) end proc

(6)

H(1,1,0,0);

2.668201175

(7)

eq1:= diff(z(m,nu,t),t)=-I*Sum(Sum(H(m,n,mu,nu)*z(n,mu,t),n=1..N),mu=0..M);

diff(z(m, nu, t), t) = -I*(Sum(Sum(((n+.5000000000)*Physics:-KroneckerDelta[m, n]-4.934802202*nu^2*Physics:-KroneckerDelta[mu, nu]+.7788007831*(factorial(min(m, n))/factorial(max(m, n)))^(1/2)*2.^(-.5000000000*abs(-1.*n+m))*LaguerreL(min(m, n), abs(-1.*n+m), -.5000000000)+39.47841762*mu*nu*(1.648721271-.6065306597*(-1.)^(mu+nu))/((9.869604404*(mu+nu)^2+1.)*(9.869604404*(mu-1.*nu)^2+1.)))*z(n, mu, t), n = 1 .. N), mu = 0 .. M))

(8)

 

zint := Array([[1,0],[0,1]]);

Matrix(2, 2, {(1, 1) = 1, (1, 2) = 0, (2, 1) = 0, (2, 2) = 1})

(9)

Z:= Array(1..2,1..2);
 for i from 1to 2 do
    for j from 1 to 2 do
        Z[i,j]:= dsolve({eq1,zint[i,j]},numeric,output=listprocedure);
    end do;
end do;

 

Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 0})

 

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

 

 

Download test3.mw

I want to plot phase portrait for 4 dimension. help me

4_equation_phase_portrait.mw

thank youvery much

There are some excel sheets that cannot be read w/ Maple 2023 that were read w/ Maple 2022.
I attached an example.ImportFerrites.maple
Does anyone else experienced this issue or has a solution?

am attaching the worksheet of the problem please help me to solve not able to  compute coupled

error.mw

in this equations...how to plot (a,e)? thank you guys

plot.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

Any idea how I catch that error?

(in geometry:-DefinedAs) wrong type of argument |C:/Users/andreas/AppData/Local/Temp/com.maplesoft/NODELibrary-0dc05fe4-0ffc9cfc/lib/NODEFunctions-13751:1976|

Here's the code:

    try
        a := map(coordinates, DefinedAs(segm))[1];
        b := map(coordinates, DefinedAs(segm))[2];
    catch "wrong type of arguments":
        Alert(cat("Error, (in geometry:-DefinedAs) wrong type of arguments, variable segm: ", whattype(segm), segm), table(), 5)
    end try;

how to write the code for integration of the orthogonal cosine function with nu=mu or with nu<>mu

restart;

 

phi:= (mu,Q2)->sqrt(2/l)*sin(mu*Pi*(Q2+l/2)/l);    # E:=mu->(Pi*mu/sqrt(2)/l)^2, mu=1,2...

proc (mu, Q2) options operator, arrow; sqrt(2/l)*sin(mu*Pi*(Q2+(1/2)*l)/l) end proc

(1)

fh1:=simplify((-1/2/m2*Int(diff(phi(mu,Q2),Q2)*diff(phi(nu,Q2),Q2),Q2=-l/2..l/2)))

-mu*Pi^2*nu*(Int(cos((1/2)*mu*Pi*(2*Q2+l)/l)*cos((1/2)*nu*Pi*(2*Q2+l)/l), Q2 = -(1/2)*l .. (1/2)*l))/(l^3*m2)

(2)

convert(fh1,int) assuming(mu,integer,nu,integer);

0

(3)

fh1_subs := simplify(subs(nu = mu, fh1));

-mu^2*Pi^2*(Int(cos((1/2)*mu*Pi*(2*Q2+l)/l)^2, Q2 = -(1/2)*l .. (1/2)*l))/(l^3*m2)

(4)

convert(fh1_substituted,int)assuming (mu,integer,nu,integer);

-(1/2)*mu^2*Pi^2/(l^2*m2)

(5)

 

Download test1.mw

Can you change f(eta) to upflow curve and theta(eta) to downflow curve.

In my Problem,Boundary Conditions are

theta(infinity) = 0, (D(f))(infinity) = 1 , (Take, eta =infinity)

Flows will be correct for what value is taken for infinity .

I take  eta = 5. and also tried changing ranges  but could't find it.Please Help to fix the curve.

my code is,

SM.mw

Hi,

I am creating a random question about the trigonometric form of a complex number. My correction code does not accept all possible amplitudes for the argument (for example, it does not accept 13Pi/6 for an argument of Pi/6). Any ideas on how to fix the code?

QComplexCis.mw

Is there a viewer in Maple to view the .mw files within the Windows or Maple directory [or alternative OS environment], not unlike a pdf or MS Word viewer?

It would make it a lot easier to check through working files without having to load all of them into the system, particulary for developmental work.

regards, RH

I am trying to find three integer a, b, c so that equation (x-a)^4 + (x-b)^4 = c has four different solutions integer x. I tried

restart;
k := 0;
for a from -10 to 10 do
    for b from -10 to 10 do
for c from -10 to 10 do
if a*b*c <> 0 and igcd(a, b) = 1 then X := [solve((x - a)^4 + (x - b)^4 = c, x)]; if type(X[1], integer) and type(X[2], integer) and type(X[3], integer) and type(X[4], integer) and nops({X[1], X[2], X[3], X[4]}) = 4 then k := k + 1; L[k] := [a, b, c, X[]]; end if; end if; end do; end do;
end do;
L := convert(L, list);
k;

L;

I do not get any solution

First 190 191 192 193 194 195 196 Last Page 192 of 2427