Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I'm having a symmetric bivariate polynomial P(x,y) with a high degree (larger than 8). I want to know if there exists a root of P(x,y) or not. I may not need a specific root but rather to know if P(x,y) has a root or not. I've been doing some research but I only find the command that helps solve univariate functions or solve systems of bivariate functions. 

In conclusion, my question is: How can I know that a symmetric bivariate polynomial has a root or not using Maple?

Here is the polynomial I want to check

P(x,y)=x^12*y^3 + 10*x^11*y^4 + 43*x^10*y^5 + 105*x^9*y^6 + 161*x^8*y^7 + 161*x^7*y^8 + 105*x^6*y^9 + 43*x^5*y^10 + 10*x^4*y^11 + x^3*y^12 - x^12*y^2 - 8*x^11*y^3 - 17*x^10*y^4 + 8*x^9*y^5 + 82*x^8*y^6 + 128*x^7*y^7 + 82*x^6*y^8 + 8*x^5*y^9 - 17*x^4*y^10 - 8*x^3*y^11 - x^2*y^12 - x^12*y - 4*x^11*y^2 - 62*x^10*y^3 - 341*x^9*y^4 - 902*x^8*y^5 - 1410*x^7*y^6 - 1410*x^6*y^7 - 902*x^5*y^8 - 341*x^4*y^9 - 62*x^3*y^10 - 4*x^2*y^11 - x*y^12 + x^12 - 8*x^11*y - 62*x^10*y^2 - 680*x^9*y^3 - 3169*x^8*y^4 - 7312*x^7*y^5 - 9540*x^6*y^6 - 7312*x^5*y^7 - 3169*x^4*y^8 - 680*x^3*y^9 - 62*x^2*y^10 - 8*x*y^11 + y^12 + 10*x^11 - 17*x^10*y - 341*x^9*y^2 - 3169*x^8*y^3 - 11838*x^7*y^4 - 21793*x^6*y^5 - 21793*x^5*y^6 - 11838*x^4*y^7 - 3169*x^3*y^8 - 341*x^2*y^9 - 17*x*y^10 + 10*y^11 + 43*x^10 + 8*x^9*y - 902*x^8*y^2 - 7312*x^7*y^3 - 21793*x^6*y^4 - 30696*x^5*y^5 - 21793*x^4*y^6 - 7312*x^3*y^7 - 902*x^2*y^8 + 8*x*y^9 + 43*y^10 + 105*x^9 + 82*x^8*y - 1410*x^7*y^2 - 9540*x^6*y^3 - 21793*x^5*y^4 - 21793*x^4*y^5 - 9540*x^3*y^6 - 1410*x^2*y^7 + 82*x*y^8 + 105*y^9 + 161*x^8 + 128*x^7*y - 1410*x^6*y^2 - 7312*x^5*y^3 - 11838*x^4*y^4 - 7312*x^3*y^5 - 1410*x^2*y^6 + 128*x*y^7 + 161*y^8 + 161*x^7 + 82*x^6*y - 902*x^5*y^2 - 3169*x^4*y^3 - 3169*x^3*y^4 - 902*x^2*y^5 + 82*x*y^6 + 161*y^7 + 105*x^6 + 8*x^5*y - 341*x^4*y^2 - 680*x^3*y^3 - 341*x^2*y^4 + 8*x*y^5 + 105*y^6 + 43*x^5 - 17*x^4*y - 62*x^3*y^2 - 62*x^2*y^3 - 17*x*y^4 + 43*y^5 + 10*x^4 - 8*x^3*y - 4*x^2*y^2 - 8*x*y^3 + 10*y^4 + x^3 - x^2*y - x*y^2 + y^3

`debugger/no_output`;

From the Maple 2023 user manual for AudioTools, when I attempt to input:

with(AudioTools):
Tone := f -> Play(Create(x -> evalhf(sin(1/22050*x*Pi*f)), duration = 1)):
Tone(440);

I get the following error message:

Error, (in DocumentTools:-GetProperty) Attempted to retrieve property of unknown component Speaker0.

Clicking the error message in Maple then goes to this page:

https://www.maplesoft.com/support/help/errors/view.aspx?path=Error,%20(in%20DocumentTools:-GetProperty)%20Attempted%20to%20retrieve%20property%20of%20unknown%20component%20Speaker0.

How do I get a function to play as an audio file in Maple 2023?

Is there a way to enclose in parentheses or brackets a highlighted selection of Maple Input?  This is a nice feature on WXMaxima to quickly and accurately enclose any selection of input syntax in parentheses or brackets.  Otherwise it can become tricky to add the opening and closing parentheses/brackets separately.

How do I insert an Input Prompt (">")?

On Maple 2023, each new input is marked on the right-hand side as (1), (2), (3), etc.  This is equivalent on WXMaxima as %i1, %i2, %i3, etc.  On WXMaxima, I can left-click between any two previous inputs to insert a new input, for example %i4 might be between %i1 and %i2 if I added %i4 after %i3.  However, I can't seem to do this on Maple.  Is there a way to go back and add a new input between two previous inputs on Maple 2023?

Likewise, I can't seem to add a new input prompt after the last input prompt if the last input prompt has not executed.

Is there a way to express any number in any number base using Maple?

Start with converting any positive real number (rational or irrational) into any positive real number base (rational or irrational), such as: what is "1/65536" in base "999" or what is "sqrt(2)" in base "exp(1)"?

But also for negative, imaginary, and complex numbers: for example, what is "(1/3)-sqrt(3)*I" in base "-Pi*I"?

Specifically, I would like to know what are the first 1024 digits of "exp(1)" in base 16777216, but a general function to convert any number into any base would be preferred.  Base 16777216 is a fun example since each digit can be expressed as a hexadecimal color (from #000000 to #ffffff) on a computer monitor.  Not sure how to get Maple to color the digits automatically though.

It's unfortunate that the "convert/base" feature on Maple only works for integers and bases that are positive integers.

https://en.wikipedia.org/wiki/golden_ratio_base

https://en.wikipedia.org/wiki/negative_base

https://en.wikipedia.org/wiki/quater-imaginary_base

Hi,

I am trying to animate quadrics, and I don't understand the issue with my fourth animation on the elliptic paraboloid

Thanks for your help

QQEllipticPar.mw

Please help me how to plot the following Chemotaxis System in Maple. Great Thanks

PDE 1= diff(p(x,t),t)-diff(p(x,t)*q(x,t),x)=D*diff(p(x,t),x,x)

PDE 2 = diff(q(x,t),t)-diff(p(x,t),x)=0

with the initial conditions

p(x)=0.5+1/(1+exp(2*(x-30)))

and

p(x)=0.5+1/(1+exp(-2*(x-30)))

and the boundary conditions

p(0, t) = 1.5, p(1, t) = 0.5

and

q(0, t) = 0.5, q(1, t) = 1.5

Certainly, ½ cannot be a root of the following equations: 
 

interface(version)

restart;

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

(1)

solvefor[x]((x-RealDomain:-`^`(2, -1))*RealDomain:-`^`(x-2*RealDomain:-`^`(3, -1), -1)*(1-sqrt(1-x*(x-2*RealDomain:-`^`(3, -1))*RealDomain:-`^`(RealDomain:-`^`(x-RealDomain:-`^`(2, -1), 2), -1))) = 0)

Warning, solvefor is deprecated. Please use solve command.

 

[x = 0, x = 1/2]

(2)

`~`[limit]((x-1/2)*(1-sqrt(1-x*(x-2/3)/(x-1/2)^2))/(x-2/3), [x = 0, x = 1/2])

[0, undefined]

(3)

solvefor[x]((x-RealDomain:-`^`(2, -1))*RealDomain:-`^`(x-2*RealDomain:-`^`(3, -1), -1)*(1-sqrt(1+3*x*(x-2*RealDomain:-`^`(3, -1))*RealDomain:-`^`(5*RealDomain:-`^`(x-RealDomain:-`^`(2, -1), 2), -1))) = 0)

Warning, solvefor is deprecated. Please use solve command.

 

[x = 0, x = 1/2]

(4)

`~`[limit]((x-1/2)*(1-sqrt(1+3*x*(x-2/3)/(5*(x-1/2)^2)))/(x-2/3), [x = 0, x = 1/2])

[0, undefined]

(5)

``


 

Download solvefor_BUG.mw

But why can't Maple's solvefor rigorously verify (or at least try to check) the solution by itself?
Please note that this issue is irrelevant to the alleged deprecated command. You may reproduce these via :-solve, Degrees:-solveRealDomain:-solvePDEtools:-Solve, etc.

Say I have 20 userdefined functions defined by me say

Now i required function which picks one of the userdefined from the above 20 functions randomly

ADM-1.mw

I need a help for solving this non linear equation by Adomian Decomposition Method.How to find A0,A1,A2... and u1,u2,u3 and a series.I am getting this error,Please help me.

restart

PDEtools[declare](prime = x);

`derivatives with respect to`*x*`of functions of one variable will now be displayed with '`

(1)

equ1 := u[tt] = 1/2*(u[xx]+u[yy])+u^2:

ICS; u(x, y, 0) = 1, u[t](x, y, 0) = e^(x+y), lambda = 0, u[0] = 1+t*e^(x+y), F(u[0]) = u[0]^2

for n from 0 to 5 do A[n] := d^n*[F*(sum(lambda^i*u[i], i = 0 .. n))]/(factorial(n)*`dλ`^n); u[n+1] = (1/2)*int[diff(sum(u[n], n = 0 .. 5), x, x)+diff(sum(u[n], n = 0 .. 5), y, y), t = 0 .. t, t = 0 .. t, dt*dt]+int[A[n], t = 0 .. t, t = 0 .. t, dt*dt] end do

``

Download ADM-1.mw

How do I display Unicode characters in Maple 2023?  From the user manual, it seems like this should be supported, but all I see are tofu characters (empty white boxes) and/or question marks.  I have the correct Windows fonts installed for Unicode support and the characters display correctly in Microsoft Notepad.

For example, in sexagesimal (base-60) cuneiform, 1 = 𒐕 (𒐕), 2 = 𒐖 (𒐖), 3 = 𒐗 (𒐗), etc.  If the cuneiform characters display on your web browser, then you should have the correct fonts installed, but they don't display on Maple.  Using the special characters as numerical constants, then you should be able to say in Maple:

𒐕 := 1;
𒐖 := 2;
𒐗 := 3;
𒐕+𒐖;

I tried adding "with(XMLTools):" and "with(StringTools):" but that didn't help, nor did "DecodeEntities("𒐕");".

Even more problematic than simply not displaying the characters, but Maple does not seem to be able to save these characters as constants.  All three characters save as the value "3" instead of 1, 2, and 3, so that the sum of 𒐕 (1) and 𒐖 (2) equals 6 instead of 3.  Something else wrong is that I see two tofu boxes and/or two question marks for a single cuneiform character.  Not sure if this is a bug with Maple 2023.

I seem to have the same problem with က (က) and 𐀀 (𐀀) but not with Ā (Ā) so it looks like maybe there is not full Unicode support for Maple 2023?  𐀀 displays as two missing characters which somehow is already saved as the decimal value 3 whereas က displays as a single missing character with no defined value.

https://en.wikipedia.org/wiki/Babylonian_cuneiform_numerals

https://en.wikipedia.org/wiki/Cuneiform_Numbers_and_Punctuation#Block

I could not figure how to fix my code. Could you help me?

Thanks in advance,

LinearAnalysis.mw

restart;

:with(DynamicSystems):with(DEtools):with(plots):with(plottools):with(DETools):

wiht(LinearAlgebra): with(VectorCalculus):

 

 

e1:=0.5: e2:=0.2:

sigma1:=e1*alpha: sigma2:=e2*delta:

g:=x/(1+beta*x^2);
f:=(theta*x-1)*(1-x)*(1+beta*x^2)-y;
h:=alpha*g-z-sigma1;
j:=delta*y-sigma2;
F:=g*f;
H:=y*h;
G:=z*j;

x/(1+.4*x^2)

 

(5*x-1)*(1-x)*(1+.4*x^2)-y

 

x/(1+.4*x^2)-z-.5

 

y-.2

 

x*((5*x-1)*(1-x)*(1+.4*x^2)-y)/(1+.4*x^2)

 

y*(x/(1+.4*x^2)-z-.5)

 

z*(y-.2)

(1)

 

theta:=5;alpha:=1;delta:=1;beta:=0.4;
Ffunc:=unapply(F,x,y); Hfunc:=unapply(H,x,y,z);Gfunc:=unapply(G,y,z);
J1:=Jacobian([Ffunc(x,y),Hfunc(x,y,z),Gfunc(y,z)],[x,y,z]);
C:= CharacteristicPolynomial(J1,lambda);
R:= RouthTable(C,lambda);

theta := 5

 

alpha := 1

 

delta := 1

 

beta := .4

 

Ffunc := proc (x, y) options operator, arrow; x*((5*x-1)*(1-x)*(1+.4*x^2)-y)/(1+.4*x^2) end proc

 

Hfunc := proc (x, y, z) options operator, arrow; y*(x/(1+.4*x^2)-z-.5) end proc

 

Gfunc := proc (y, z) options operator, arrow; z*(y-.2) end proc

 

Matrix(3, 3, {(1, 1) = ((5*x-1)*(1-x)*(1+.4*x^2)-y)/(1+.4*x^2)-.8*x^2*((5*x-1)*(1-x)*(1+.4*x^2)-y)/(1+.4*x^2)^2+x*((5-5*x)*(1+.4*x^2)-(5*x-1)*(1+.4*x^2)+(4.0*x-.8)*(1-x)*x)/(1+.4*x^2), (1, 2) = -x/(1+.4*x^2), (1, 3) = 0, (2, 1) = y*(1/(1+.4*x^2)-.8*x^2/(1+.4*x^2)^2), (2, 2) = x/(1+.4*x^2)-z-.5, (2, 3) = -y, (3, 1) = 0, (3, 2) = z, (3, 3) = y-.2})

 

Error, (in Verify:-CommonExports) system object is not a module

 

Matrix(%id = 36893491070243237028)

(2)

E0:= [0,0,0];
JE0:=subs(x=E0[1],y=E0[2],z=E0[3],J1);
CE0:=CharacteristicPolynomial(JE0,lambda);
RE0:= RouthTable(CE0,lambda);
RouthTable(CE0,lambda,'stablecondition'=true);

E0 := [0, 0, 0]

 

Matrix(3, 3, {(1, 1) = -1., (1, 2) = -0., (1, 3) = 0, (2, 1) = 0., (2, 2) = -.5, (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = -.2})

 

Error, (in Verify:-CommonExports) system object is not a module

 

Matrix(%id = 36893491070243233172)

 

true

(3)

E1:= [1,0,0];
JE1:=subs(x=E1[1],y=E1[2],z=E1[3],J1);
CE1:=CharacteristicPolynomial(JE1,lambda);
RE1:= RouthTable(CE1,lambda);
RouthTable(CE1,lambda,'stablecondition'=true);

E1 := [1, 0, 0]

 

Matrix(3, 3, {(1, 1) = -4.000000000, (1, 2) = -.7142857143, (1, 3) = 0, (2, 1) = 0., (2, 2) = .2142857143, (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = -.2})

 

Error, (in Verify:-CommonExports) system object is not a module

 

Matrix(%id = 36893491070243221500)

 

true

(4)

E2:= [1/theta,0,0];
JE2:=subs(x=E2[1],y=E2[2],z=E2[3],J1);
CE2:=CharacteristicPolynomial(JE2,lambda);
RE2:= RouthTable(CE2,lambda);
RouthTable(CE2,lambda,'stablecondition'=true);

E2 := [1/5, 0, 0]

 

Matrix(3, 3, {(1, 1) = .8000000000, (1, 2) = -.1968503937, (1, 3) = 0, (2, 1) = 0., (2, 2) = -.3031496063, (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = -.2})

 

Error, (in Verify:-CommonExports) system object is not a module

 

Matrix(%id = 36893491070243218004)

 

true

(5)

 

Download LinearAnalysis.mw

 

I have been using WXMaxima for years and just upgraded to Maple 2023, but I am having a very difficult time adjusting from Maxima to Maple.  Is there a way to make Maple more Maxima-like?

The first problem that I encounter is that "Shift + Enter" and "Enter" should be reversed for people who are used to inputting with Maxima.  I would like to reconfigure Maple so that "Shift + Enter" will be for "Evaluate Expression" and "Enter" will be for "Soft New Line".  This makes inputting multi-line commands (using ":") much easier, so that you "Enter" to add the next line, and then "Shift + Enter" to evaluate.  I thought there was a way to switch these, but I can't find it in the manual.

The next problem is the syntax itself.  "MmaToMaple" will translate Mathematica syntax into Maple syntax, but why is there no equivalent for Maxima?  Maxima/Macsyma is the original computer algebra system (CAS) that inspired both Maple and Mathematica, and being free and open-source, has the widest possible user base (since anyone can use it for free).  Maple needs a way to translate to and from Maxima syntax.  For personal use, I already have a large body of Maxima worksheets saved as plaintext MAC files that I would like to use in Maple.  It will save me a lot of time and trouble if I can import Maxima worksheets directly into Maple.

Is there an easy way to save commonly used Maxima expressions that are not already reserved by Maple, such as "float" (for "evalf") and "%e" (for "exp(1)") so that some Maxima syntax can run directly in Maple?  The only major changes then would be replacing "$" with ":" and ":" with ":=" since ";" is the same in both Maxima and Maple.  Not sure what a "$" means in Maple....

Maple syntax is more similar to Maxima syntax than to Mathematica syntax.  Below is a sample problem that is un-solvable in Maxima (application will crash or run out of memory) but solvable in Maple, to show the differences in syntax.  Unfortunately, it seems that in many cases the Maple syntax is more verbose than the Maxima syntax, requiring more characters to type (such as when defining functions).

/* Maxima syntax */ $
c : 299792458 $
h : 6.62607015E-34 $
k : 1.380649E-23 $
n : exp(1) $
T : 25200 $
P (lmin,lmax) := integrate(2*(c^2)*h/(l^5*(n^(c*h/(k*l*T))-1)),l,lmin,lmax)/integrate(2*(c^2)*h/(l^5*(n^(c*h/(k*l*T))-1)),l,0,inf) $
float(P(380E-9,750E-9));

 

(* Maple syntax *) :
c := 299792458 :
h := 6.62607015E-34 :
k := 1.380649E-23 :
n := exp(1) :
T := 25200 :
P := (lmin,lmax) -> integrate(2*(c^2)*h/(l^5*(n^(c*h/(k*l*T))-1)),l=lmin..lmax)/integrate(2*(c^2)*h/(l^5*(n^(c*h/(k*l*T))-1)),l=0..infinity) :
evalf(P(380E-9,750E-9));

It seems like Maplesoft has gone to a lot of effort to create intercompatibility between Maple and Mathematica, but their real competition is getting customers to upgrade from Maxima (which is completely free and open-source) to Maple (which is non-free and closed-source like Mathematica is).  Making Maple more Maxima-friendly should increase the number of Maple customers by making it easier to upgrade from Maxima to Maple.

Can anyone help me fix a corrupted maple file?

I have tried to fix it myself and searched all over for a solution but can' seem to fix it myself. Please help

Download Styrkelære2regneopgaver.mwStyrkelære2regneopgaver.mwStyrkelære2regneopgaver.mw

First 92 93 94 95 96 97 98 Last Page 94 of 2097