Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

As by title, is it possible to check intermediate results found by nodes when using Grid:Map, and stop the execution according to these?

every time i have a small problem 2  why not cancel number 2 in denominator , i don't want see a number with fraction like 1/3  3/4 how i fix this

restart

K := (1/2)*sqrt(2*sqrt(2)*sqrt(lambda*a[5]/a[4])+2*sqrt(-2*a[5]/a[4])*(B[1]*sqrt(-lambda)*sinh(sigma)+B[2]*sqrt(-lambda)*cosh(sigma))/(B[1]*cosh(sigma)+B[2]*sinh(sigma)+mu/lambda)+2*sqrt(-(2*(lambda^2*B[1]^2*a[5]-lambda^2*B[2]^2*a[5]-mu^2*a[5]))/(lambda*a[4]))/(B[1]*cosh(sigma)+B[2]*sinh(sigma)+mu/lambda))*e^(i*psi(x, t))

Warning, if e is meant to be the exponential e, use command/symbol completion or palettes to enter this special symbol, or use the exp function

 

(1/2)*(2*2^(1/2)*(lambda*a[5]/a[4])^(1/2)+2*(-2*a[5]/a[4])^(1/2)*(B[1]*(-lambda)^(1/2)*sinh(sigma)+B[2]*(-lambda)^(1/2)*cosh(sigma))/(B[1]*cosh(sigma)+B[2]*sinh(sigma)+mu/lambda)+2*(-2*(lambda^2*B[1]^2*a[5]-lambda^2*B[2]^2*a[5]-mu^2*a[5])/(lambda*a[4]))^(1/2)/(B[1]*cosh(sigma)+B[2]*sinh(sigma)+mu/lambda))^(1/2)*e^(i*psi(x, t))

(1)

B[1] := 0; mu := 0

0

 

0

(2)

simplify(eval(K))

(1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t))

(3)

simplify((1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t)), 'trig')

(1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t))

(4)

simplify((1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t)))

(1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t))

(5)

simplify((1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t)), 'constant')

(1/2)*2^(3/4)*(((coth(sigma)*(-lambda)^(1/2)*(-a[5]/a[4])^(1/2)+(lambda*a[5]/a[4])^(1/2))*B[2]+csch(sigma)*(lambda*a[5]*B[2]^2/a[4])^(1/2))/B[2])^(1/2)*e^(i*psi(x, t))

(6)
 

 

 

NULL

Download cancelation.mw

Recently a new configuration of a computer, with a better cpu - AMD Ryzen 7 8700G w/ Radeon 780M Graphics, but maple computation is particularly slow, cpu call special less computation is particularly slow, is there any way to improve?

OneFrame := proc(k)
local Courbe, T, a, b, c, t, P, Q, NormM, F, Ell, sol, N1, N2, dr, tx;
a := 11; b := 7; c := sqrt(a^2 - b^2); t := 1/3*Pi;
Ell := x^2/a^2 + y^2/b^2 = 1;
geometry:-point(T, (a^2 - b^2)*cos(t)^3/a, -(a^2 - b^2)*sin(t)^3/b);
Courbe := plots:-implicitplot(Ell, x = -a - 10 .. a + 10, y = -b - 10 .. b + 10, scaling = constrained, color = blue);
NormM := plots:-implicitplot(y - b*sin(t) = a*sin(t)*(x - a*cos(t))/(b*cos(t)), x = -a - 5 .. a + 10, y = -b - 10 .. b + 10, color = orange); geometry:-line(Per, y - b*sin(t) = a*sin(t)*(x - a*cos(t))/(b*cos(t)), [x, y]);
geometry:-point(P, subs(y = 0
, geometry:-Equation(Per), 0));
geometry:-point(Q, 0, subs(x = 0, geometry:-Equation(Per)));
geometry:-point(M, a*cos(t), b*sin(t));
geometry:-point(N1, a*cos(k), b*sin(k));
geometry:-point(F, 2.329411765, -2.567510609);
geometry:-line(L, N1, F);
sol := solve({geometry:-Equation(L), Ell}, {x, y},explicit);
geometry:-point(N2, subs(sol[2], x), subs(sol[2], y));
geometry:-segment(sg, N1, N2);
tx := plots:-textplot([[geometry:-coordinates(M)[], "M"],
[geometry:-coordinates(N1)[], "N1"], [geometry:-coordinates(N2)[], "N2"],
[geometry:-coordinates(P)[], "P"],
[geometry:-coordinates(Q)[], "Q"],
[geometry:-coordinates(F)[], "F point de Frégier"],
[geometry:-coordinates(T)[], "T"]], font = [times, bold, 16], align = [above, left]);
dr := geometry:-draw([sg(color = magenta, linestyle = dash),
Per(color = black), P(color = red, symbol = solidcircle, symbolsize = 12),
Q(color = red, symbol = solidcircle, symbolsize = 12),
M(color = black, symbol = solidcircle, symbolsize = 12),
F(color = red, symbol = solidcircle, symbolsize = 12),
N1(color = black, symbol = solidcircle, symbolsize = 8),
N2(color = black, symbol = solidcircle, symbolsize = 8),
T(color = black, symbol = solidcircle, symbolsize = 8)]);
plots:-display(Courbe, tx, dr, scaling = constrained, axes = none); end proc;

plots:-animate(OneFrame, [k], k = Pi/3 .. Pi, frames = 50);
Error, (in plots/animate) wrong type of arguments
Why this animation does't work ? Thank you very much.
 

 

 

Any ideas, guys?

I am trying to create a ring of difference operators with customly defined actions. I have managed to get the code to work to an extent, but I am having trouble creating the inverse operators. For example, if Sn sends n to n+2, I want to create Sn^(-1) that sends n to n-2. When I tried to input this into the algebra itself, it gave me the following error: "Error, (in Ore_algebra:-skew_algebra) indeterminate n may appear in a single commutation only". Below is the code that creates the difference algebra. Apologies for the bad formatting.


> with(Ore_algebra);

> CreateDifferenceAlgebra := proc (vars::(list(symbol)), offsets::(list(numeric)))


local actions, i, var, diff_op,
shift_offset, algebra; actions := table();
for i to numelems(vars) do
   var := vars[i];
   diff_op := convert(cat("S",\var), symbol);
   shift_offset := offsets[i];
   actions[diff_op] := (proc (values) options operator, arrow; proc (u, order)
      local res; res := u;
      to order do res := subs(values[1] = values[1]+values[3], res*values[2]) end do;
      res
      end proc

   end proc)([var, diff_op, shift_offset]) end do;


algebra := Ore_algebra:-skew_algebra(seq(shift = [convert(cat("S",

vars[i]), symbol), vars[i]], i = 1 .. numelems(vars)),

action = {seq(convert(cat("S", vars[i]), symbol) =
actions[convert(cat("S", vars[i]), symbol)], i = 1 ..
numelems(vars))}); return algebra end proc


;

 

 

I have currently a Maple session running with only one open worksheet.

The worksheet only contains an input line with the name "a" and an output line with the name "a".

The fan of my laptop is running full throttle and the task manager displays 10% CPU usage and a very high power usage (3 mserver.exe running, the process with the cpu load is javaw.exe. Suspending this process shuts down the fan).

After about an hour I decided to ask. While typing this post (in Firefox, probably unrelated) the CPU usage went down and the fan went quite. The total system CPU usage is now down to 2%.

Has anybody seen the same? What could have cause it? It's not the first time I observe this. Anything that I could check before restarting Maple? All on Windows 10 after system restart.

Update:

The thread that consumes cpu is called ucrtbase.dll!configthreadlocate

Starting Maplesim in parallel almost immediately turns off the fan while the displayed cpu load is still high but now variing.

I can't understand why Maple interprets 1 .05  as 1 * 5 = 5 , and 2 .05  as 2 * 5 = 10 . Note the space between 1 and .05
In a different calculation I accidentally inserted a space between 1 and .05, and received a strange answer, and finally narrowed it down to a space.
But now I wonder why it is it is interpreting this way. Also I see that 2  0.05 produces an error. But  2 .05 is treated as 10. There is an implied multiplication? But the multiplication should be 2 * .05 , which is 1 and not 10.

i found thus condition which if we substitute in equation must be equal to zero, i don't know  how i can get zero

test_pde1.mw

I have encountered an issue: eq1 is not satisfied, though eq2 is satisfied for the parametric value (10). I need assistance in finding a way to ensure that both equations are satisfied simultaneously. Please provide guidance or suggest a potential approach for addressing this issue.verf_kk.mw

I have  a big problem in transformation How we can do suh transformation in  type of  procure  without use any hand work for example in physic abs|-| remove the exponential term how the maple remove that term automatically and collect all term and do my transformation this example is really hard one which is must do a lot by hand and mixed them which maybe a week take my time to get results and how i reach the results without spending that time i have a result of this equation and i am try to get but i don't know the results of this person is correct or not but i will share in here,  i did some try i will share in here too if in DEchange add U(xi) it will work and give me the other step but i need something more effective, when q^* is conjugate of q =exp(-ipsi(x,t))U(xi)

NULL

restart

with(PDEtools)

with(Physics)

with(SolveTools)

undeclare(prime)

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

(1)

 

 

tr := {t = tau, x = xi/k+v*tau^alpha/(k*alpha)+theta, u(x, t) = U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta)), u[1](x, t) = U(xi)*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))}

{t = tau, x = xi/k+v*tau^alpha/(k*alpha)+theta, u(x, t) = U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta)), u[1](x, t) = U(xi)*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))}

(2)

pde := I*(I*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)-mu*tau+theta))*w-exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))*v)+a*(diff(u(x, t), `$`(x, 2)))+b*U(xi)^2*u(x, t)+C[1](h[1]*(diff(u(x, t), `$`(x, 4)))+h[2]*(diff(u(x, t), x))^2*u[1](x, t)+h[3]*abs(diff(u(x, t), x))^2*u(x, t)+h[4]*U(xi)^2*(diff(u(x, t), `$`(x, 2)))+h[5]*u(x, t)^2*(diff(u[1](x, t), `$`(x, 2)))+h[6]*U(xi)^4*u(x, t))+I*C[2]*(h[7]*(diff(u(x, t), `$`(x, 4)))+h[8]*U(xi)^2*(diff(u(x, t), x))+h[9]*u(x, t)^2*(diff(u[1](x, t), x))) = 0

I*(I*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)-mu*tau+theta))*w-exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))*v)+a*(diff(diff(u(x, t), x), x))+b*U(xi)^2*u(x, t)+C[1](h[1]*(diff(diff(diff(diff(u(x, t), x), x), x), x))+h[2]*(diff(u(x, t), x))^2*u[1](x, t)+h[3]*abs(diff(u(x, t), x))^2*u(x, t)+h[4]*U(xi)^2*(diff(diff(u(x, t), x), x))+h[5]*u(x, t)^2*(diff(diff(u[1](x, t), x), x))+h[6]*U(xi)^4*u(x, t))+I*C[2]*(h[7]*(diff(diff(diff(diff(u(x, t), x), x), x), x))+h[8]*U(xi)^2*(diff(u(x, t), x))+h[9]*u(x, t)^2*(diff(u[1](x, t), x))) = 0

(3)

``

PDEtools:-dchange(tr, pde, [xi, tau, U, U(xi)])

I*(I*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)-mu*tau+theta))*w-exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))*v)+a*((2*I)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))/k+exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(U(xi), xi), xi))-U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k^2)*k^2+b*U(xi)^3*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))+C[1](h[1]*(-(4*I)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))/k^3-6*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(U(xi), xi), xi))/k^2+(4*I)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(diff(U(xi), xi), xi), xi))/k+exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))+U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k^4)*k^4+h[2]*(exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))+I*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k)^2*k^2*U(xi)*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))+h[3]*abs((exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))+I*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k)*k)^2*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))+h[4]*U(xi)^2*((2*I)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))/k+exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(U(xi), xi), xi))-U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k^2)*k^2+h[5]*U(xi)^2*(exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta)))^2*((diff(diff(U(xi), xi), xi))*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))-(2*I)*(diff(U(xi), xi))*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k-U(xi)*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k^2)*k^2+h[6]*U(xi)^5*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta)))+I*C[2]*(h[7]*(-(4*I)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))/k^3-6*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(U(xi), xi), xi))/k^2+(4*I)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(diff(U(xi), xi), xi), xi))/k+exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))+U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k^4)*k^4+h[8]*U(xi)^2*(exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))*(diff(U(xi), xi))+I*U(xi)*exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k)*k+h[9]*U(xi)^2*(exp(I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta)))^2*((diff(U(xi), xi))*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))-I*U(xi)*exp(-I*(xi/k+v*tau^alpha/(k*alpha)+mu*tau+theta))/k)*k) = 0

(4)
 

NULL


Download find_ODE.mw

How can I check if a name has been used/entered already but was not assigned to a value

The variable palette only lists assigned names.

I tried unames() but this lists all unassigned names. A 'user' option (which filters for user-assigned names) as in anames() does not seem to exist.

One of my failed attempts (in 1D-Math):

restart;
unames():
initial_unames := {%}:
new_name;
{unames()} minus initial_unames; # should ideally return a reduced set containing new_name;
has(%,new_name)

What else can be done? (I am probably overlooking something very simple.)

@Rouben Rostamian  

Dear Sir Professor Rostamian my name is Viorel Popescu from the Polytechnic University of Bucharest if you remember in the summer of 2019 you helped me to solve the equation: rH''(r)+H'(r)+(rk^2-r^2*b^2/R^2)H(r)=0 where k, b, and R are real constants positive number, with condition H(R)=0 and H'(1/R)=R. I appreciate it very much, please I'm in a similarly embarrassing situation to beg you for an answer. I want to find the equation of audion and complete the experiment http://www.michaelvio.byethost8.com/Audion.pdf

My account in Maple Primes is the same michaelvio (35) as the email michaelvio@yahoo.com and also @gmail.com it's an experiment that I want to make for my PhD. Practically I suppose that the energy can be approx. as a series of power of frequency t from I selected severaral terms Ea := 0.00762014687*t + a*t^2 + b*t^3 + c*t^4 + d*t^5 and I guess that satisfies an equation as in the document. The case of photons is beyond my possibility, but a little help from a distinguished Professor as you should cheer me up Audion1.mw

Audion.docx

Please help! 

Hello,

I noticed that the Linearly Implicit Euler method (also known as the Semi-Implicit Euler method) is not available in Maple's built-in ODE solvers. This method is useful for stiff ODEs, where part of the function is treated implicitly (for the linear term) and part is treated explicitly (for the non-linear term).

I know that the Linearly Implicit Euler method is a specialized method that probably does not find enough widespread use to justify its inclusion as a standard feature in Maple, especially given Maple's focus on numerical methods such as Runge-Kutta methods and fully implicit methods for rigid equations.

I’m wondering:

  1. Why isn’t this method included in Maple’s standard set of numerical solvers?
  2. How can I implement this method in my own code in Maple to solve stiff ODEs?

Any guidance or examples of implementation would be greatly appreciated!

Thank you!

Linearly_Implicit_Method.pdf

The usual ODE must be solved:
y´´*(y^3-y)+y´^2 *(y^2+1)=0
"Dangerous places" of the definition domain must be described: Where are the general solution y(x) and its derivatives continuous?

First 11 12 13 14 15 16 17 Last Page 13 of 2215