MaplePrimes Questions

Drear freinds,

I want to simplify f (a long experssion) in the form of f2. How to determine M1^2 and M0^2?

f1.mw

Dear all,

I try to run the test.java example but still fails with the error:

Error loading libraries: java.lang.UnsatisfiedLinkError: no jopenmaple in java.library.path: /Library/Frameworks/Maple.framework/Versions/2024/java:/Library/Frameworks/Maple.framework/Versions/2024/java:/Library/Frameworks/Maple.framework/Versions/2024/java

the file jopenmaple.jar is located in 

/Library/Frameworks/Maple.framework/Versions/2024/java

the maple part of my .bash_profile :

# set Maple Home

export MAPLEDIR="/Library/Frameworks/Maple.framework/Versions/2024"

export DYLD_LIBRARY_PATH="$DYLD_LIBRARY_PATH:<BINDIR>"

export MAPLE="<MAPLEDIR>"

export DYLD_LIBRARY_PATH="/Library/Frameworks/Maple.framework/Versions/2024/bin.APPLE_ARM64_MACOS"

export CLASSPATH="/Library/Frameworks/Maple.framework/Versions/2024/JAVA"

The eclipse part:

JDK 22 and Compiler 22

Classpath Variables in Eclipse Java Build Path Settings:

DYLD_LIBRARY_PATH /Library/Frameworks/Maple.framework/Versions/2024/bin.APPLE_ARM64_MACOS

JAVAC /Library/Frameworks/Maple.framework/Versions/2024/java/Maple.jar

MAPLEDIR /Library/Frameworks/Maple.framework/Versions/2024

So I don't know what to do.

I want a plot of the function & the approx. calculus of integral:

E0 := evalf(int(T2, x = x0 .. x0 + 1.542976947*10^(-13))); it doesn't compute in in normal time...Audion.mw

in the program:

restart;
a := -1.44670357887361*10^(-7);
b := -1.049267156*10^(-9);
c := 1.890440485*10^(-12);
d := -6.233924848*10^(-16);
Ea := 0.00762014687*t + a*t^2 + b*t^3 + c*t^4 + d*t^5;
E1 := diff(Ea, t);
E2 := subs(t = 435, Ea);
E3 := subs(t = 528, Ea);
E4 := subs(t = 2860, Ea);
 

It seems that Maple needs more help than necessary:

restart:

kernelopts(version)

`Maple 2015.2, APPLE UNIVERSAL OSX, Dec 20 2015, Build ID 1097895`

(1)

expr:= A+B*limit(f(x), x=+infinity);
eval(expr, limit(f(x), x=+infinity)=1)

A+B*(limit(f(x), x = infinity))

 

A+B

(2)

expr:= A+B*limit(2*f(x), x=+infinity);

eval(expr, limit(f(x), x=+infinity)=1);     # Shouldn't this return A+2*B
eval(expr, limit(2*f(x), x=+infinity)=2);   # Can I avoid doing this?

A+B*(limit(2*f(x), x = infinity))

 

A+B*(limit(2*f(x), x = infinity))

 

A+2*B

(3)

expr:= A+B*limit(f(x)^2, x=+infinity);

eval(expr, limit(f(x), x=+infinity)=1);     # Shouldn't this return A+B
eval(expr, limit(f(x)^2, x=+infinity)=1);   # Can I avoid doing this?

A+B*(limit(f(x)^2, x = infinity))

 

A+B*(limit(f(x)^2, x = infinity))

 

A+B

(4)

expr:= A+B*limit(2*f(x)^2, x=+infinity);

eval(expr, limit(f(x)^2, x=+infinity)=1);    # Shouldn't this return A+2*B
eval(expr, limit(2*f(x)^2, x=+infinity)=2);  # Can I avoid doing this?

A+B*(limit(2*f(x)^2, x = infinity))

 

A+B*(limit(2*f(x)^2, x = infinity))

 

A+2*B

(5)
 

 

Download limits.mw

Why don't the commands labelled "Shouldn't this return.." do the job?

TIA

Hi everyone...

How can the coefficients of below nonlinear equation be obtained?

c1*diff(f(x), x$4) + c2*diff(f(x), x$3)/x +c3*diff(f(x), x$2)/x^2 + c4*diff(f(x), x)/x^3 +c5 *f(x)=0

I want to calculate separately, c1,c2,c3,c4,c5 

tnx...

How can I substitute equation (5) into equation (4) to obtain the final form [see image in BE.mw]?

Additionally, how can I use subscript notation for derivatives, such as expressing diff(f(tau,sigma), tau)  as  f[tau]​?

By using the 'declare' command in Maple, one can utilize subscript notation. However, in this case, the function must be written in terms of τ and σ on each step.

In the final plot in this worksheet, the electric field is displayed using blue headless vectors.

Instead, how can the gradient of the dipole's potential field be displayed in field lines?

Dipole_fields.mw

Hi, I've been trying to get help from the Maple website, but for months now, the display hasn't worked properly. The examples section shows "comma" and "colon" texts in a way that makes it impossible to use the site effectively. Has anyone else reported this issue? It's quite frustrating.

What is the correct idiom in Maple to check that one number is greater than another or not?

For example if n=sqrt(3) and m=3 ?

Maple does not like to compare sqrt(3) with other number. It says 

   Error, cannot determine if this expression is true or false: 3 <= 3^(1/2)

I am doing this in code, so solution has to be such that it works for all cases of n and m, without the ability to look at screen and decide.

I found that using is works without the need to convert to float.  Should one then use is(n>=m) instead of evalb(n>=m) always?  

I can also always apply evalf() on each side before. But this seems like an overkill to me.

restart;

n:=sqrt(3);
m:=3;

3^(1/2)

3

if evalb(n>=m) then
   "yes";
else
   "no";
fi;

Error, cannot determine if this expression is true or false: 3 <= 3^(1/2)

if n>=m then
   "yes";
else
   "no";
fi;

Error, cannot determine if this expression is true or false: 3 <= 3^(1/2)

if evalb(evalf(n)>=m) then
   "yes";
else
   "no";
fi;

"no"

if is(n>=m) then
   "yes";
else
   "no";
fi;

"no"

 

 

Download checking_one_number_larger_than_another.mw

For reference, another software does compare these two numbers as is without the need to convert them to floating point;

But it looks the design of Maple in this aspect is different. Since the principal root of sqrt(3) is always the positive one, I did not think it will cause a problem.

how fixed this for ode test

restart

with(PDEtools)

with(Physics)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

``

pde := -I*(diff(U(xi), xi))*gamma*k*mu+I*gamma*(diff(U(xi), xi))*sigma*w+(diff(diff(U(xi), xi), xi))*gamma*k*w+U(xi)*gamma*mu*sigma+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*w-U(xi)*sigma^2-U(xi)*mu

-I*gamma*(diff(U(xi), xi))*k*mu+I*gamma*(diff(U(xi), xi))*sigma*w+gamma*(diff(diff(U(xi), xi), xi))*k*w+gamma*U(xi)*mu*sigma+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*w-U(xi)*sigma^2-U(xi)*mu

(2)

case1 := [mu = -(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1), A[0] = 0, A[1] = -RootOf(_Z^2*alpha+gamma*k*w+k^2), B[1] = RootOf(_Z^2*alpha+gamma*k*w+k^2), w = (gamma*k*mu-2*k*sigma)/(gamma*sigma-1)]

[mu = -(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1), A[0] = 0, A[1] = -RootOf(_Z^2*alpha+gamma*k*w+k^2), B[1] = RootOf(_Z^2*alpha+gamma*k*w+k^2), w = (gamma*k*mu-2*k*sigma)/(gamma*sigma-1)]

(3)

G1 := U(xi) = 2*RootOf(_Z^2*alpha+gamma*k*w+k^2)/sinh(2*xi)

U(xi) = 2*RootOf(_Z^2*alpha+gamma*k*w+k^2)/sinh(2*xi)

(4)

pde1 := subs(case1, pde)

I*gamma*(diff(U(xi), xi))*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)+I*gamma*(diff(U(xi), xi))*sigma*(gamma*k*mu-2*k*sigma)/(gamma*sigma-1)+gamma*(diff(diff(U(xi), xi), xi))*k*(gamma*k*mu-2*k*sigma)/(gamma*sigma-1)-gamma*U(xi)*(4*gamma*k*w+4*k^2-sigma^2)*sigma/(gamma*sigma-1)+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*(gamma*k*mu-2*k*sigma)/(gamma*sigma-1)-U(xi)*sigma^2+U(xi)*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)

(5)

pde2 := subs(case1, pde1)

I*gamma*(diff(U(xi), xi))*k*(4*gamma*(gamma*k*mu-2*k*sigma)*k/(gamma*sigma-1)+4*k^2-sigma^2)/(gamma*sigma-1)+I*gamma*(diff(U(xi), xi))*sigma*(-gamma*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)-2*k*sigma)/(gamma*sigma-1)+gamma*(diff(diff(U(xi), xi), xi))*k*(-gamma*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)-2*k*sigma)/(gamma*sigma-1)-gamma*U(xi)*(4*gamma*(gamma*k*mu-2*k*sigma)*k/(gamma*sigma-1)+4*k^2-sigma^2)*sigma/(gamma*sigma-1)+(2*I)*(diff(U(xi), xi))*k*sigma+2*alpha*U(xi)^3+(diff(diff(U(xi), xi), xi))*k^2-I*(diff(U(xi), xi))*(-gamma*k*(4*gamma*k*w+4*k^2-sigma^2)/(gamma*sigma-1)-2*k*sigma)/(gamma*sigma-1)-U(xi)*sigma^2+U(xi)*(4*gamma*(gamma*k*mu-2*k*sigma)*k/(gamma*sigma-1)+4*k^2-sigma^2)/(gamma*sigma-1)

(6)

odetest(G1, pde2)

 

NULL

Download test_sol_for_PDE1.mw

A task that was famous at the time is worth remembering:

If for whole numbers x and y the number N = (x^2+y^2)/(1+x*y) is a positive whole number, then it is also a square number.

It can be proven that the converse is also true. Therefore, here is the task:

If N is a square number, then the Diophantine equation has solutions. Solutions must be calculated for N = 9, 49 and 729.

Hi,

I want to display the solutions in the form of intervals or empty sets, using mathematical symbols. Any suggestions?

ThanksQEnsembleAffichage.mw

Could anyone help me to convert a code written in Mathematica to a Maple worksheet? I have PDF only. Could any one have a look on Mpale sheet and PDF....

Mathematica__to_Maple.mw

Mathematica_file.pdf

An animation only allows one time parameter.

I have an setup in which I have a time parameter + another parameter or more(e.g., offset, scale, whatever).

I would like some way to easily traverse phase space and see the result and animate with time at those specific parameters(ideally animate along any path in the phase space).

In my case what this means is that there is a 2D rectangle in which represents the phase space. E.g., time x offset. This effects the animation by setting the current parameters. If the animation is animated in time it moves across the vertical. Basically a "2d slider". Alternatively have actually 2(or more) independent sliders.

Currently I have to manually set the parameter then re-execute the animation to see the new animation.

e.g.,

  animation(plot, [f(t,offset)], t=0..1,offset=1..5)

aa1.mw      aa2.mw 

In my view, these two processes are similar, but they yield different outcomes. What accounts for this discrepancy?

First 13 14 15 16 17 18 19 Last Page 15 of 2374