MaplePrimes Questions

hi...

how can I just find the real roots of this equation in Maple:

eq:=-0.0004*x^2 - 2.7680*10^(-28)*x^12 - 2.1685*10^(-43)*x^18 - 1.3245*10^(-37)*x^16 - 1.6917*10^(-32)*x^14 + 0.7650 + 6.6773*10^(-18)*x^8 - 2.5543*10^(-23)*x^10 - 8.0002*10^(-13)*x^6 + 3.6079*10^(-8)*x^4 = 0:

I want to solve it first with fsolve (with options) command and second other command than fsolve. 

tnx...

I am trying to plot the shearforce of a hollow section circular beam along the hight of the cross section.

I don't seem to be able to solve the equation. I get the message:"Warning, solutions may have been lost"
Can anyone help?

The tricky part is the distance of the "point of grafity of the area above y' to the 'point of grafity of the cross section (y=0)'. This distance is a function of y and is called y_

this is de code: M_V_as_function_of_hight_cc__temporary.mw

restart: with(LinearAlgebra): with(plots, textplot, display): #Digits :=5: evalf(%): with(RealDomain):
Rout:=168.3/2:
Rin :=Rout-12.5:

T:=V*Q/(II*b):
Q:=A_*y_:
bout:=2*sqrt(Rout**2 - y**2):
bin:=piecewise( y<Rin,2*sqrt(Rin**2 - y**2)  ,  y>=Rin,0):
#bin:=Re(2*sqrt(Rin**2 - y**2)):
b:=bout-bin;
A_ :=piecewise(  y<=0,0   ,    y<Rin,evalf( int(b ,  y=y..Rout ) ) ,    y>=Rin,int(bout , y=y..Rout)    ,    y>=Rout,0);
y:=y_: AA_:=A_;
y:='y':
eq1:=2*AA_ = A_:
sol1:=solve(eq1, y_ );
y_ := sol1:
y_;
#plot(y_,y=0..Rout);

Let us say f(x) is a procedure that can take two possible values for x.
How do I store the stats of CodeTools:-Usage(f(1))  and CodeTools:-Usage(f(2)) for ease of comparison?  Hopefully, I am not missing any help statements/examples.

Hi all
I have a simple problem with the following matrix entries. I probably have a problem with the indices. Because the matrix is ​​not calculated correctly. Anyone have suggestion?

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.

First 64 65 66 67 68 69 70 Last Page 66 of 2425