How can I use Maple to solve a system of nonlinear equations symbolically and display the steps in the solution?

Here are two systems over the reals:

The following results indicate that both and are satisfiable 

QuantifierElimination:-QuantifierEliminate(exists([s,q,r],sys__1));
                              true

QuantifierElimination:-QuantifierEliminate(exists([s,q,r],sys__1));
                              true

but RealDomain:-solve simply returns an empty list (that is, no solution exists) in both cases

RealDomain:-solve(sys__1,[q,s,r]); # ⟹ sys1 cannot be satisfied
                               []

RealDomain:-solve(sys__2,[q,s,r]); # ⟹ sys2 cannot be satisfied
                               []

As discussed in the previous problem, in contrast to using QuantifierElimination:-QuantifierEliminate, the use of RealDomain:-solve is unsafe. Nevertheless, the above output suggests that even the much-more-sophisticated QuantifierElimination:-QuantifierEliminate is still not always reliable (since the correct returnedvalue appears to be in lieu of ). So, what is the right command to handle polynomial systems over real domains in Maple? 

I have learned that the SPECTRA.mla library can solve SDP problems. I have tried to download and use it, but I am still missing the FGb module

How can I resolve this issue on a Windows environment, beause I don't see install file for Windows, only MacOS and Linux:

• FGb page: https://www-polsys.lip6.fr/~jcf/FGb/FGb/index.html
• Spectra page: https://homepages.laas.fr/henrion/software/spectra/

I have this problem:

(x'+x''+epsilon/x)(x'-(x')^3/3!)(x''-(x'')^3/3!)(epsilon/x -(epsilon/x)^3/3!)=1

x(0)=epsilon , x'(epsilon)=epsilon+epsilon^{-1}

0<epsilon<<1

Now I want to try the ansatz: x(epsilon,t)=x_{-1}(t)/epsilon+x_0(t)+x_1(t)epsilon.

Can you suggest me a way to implement this ansatz (i.e finding x_{-1}(t),x_0(t),x_1(t)).

Appreciate your time!

Hi Dear,

I hope everyone is fine here. In the attached file, I have generated a square matrix "Q" using two-dimensional polynomials. The dimension of the square matrix "Q" depends on M1 and M2 parameters. In my simulation, sometimes I need this matrix of 1000 by 1000 dimensions. Using the attached method, it took a lot of time to compute two-dimensional polynomials and then to compute the general square matrix "Q." I wanted to write this matrix using proc (procedures). Maybe by using this way, I don't need to compute the polynomials, and it took less time to compute the square matrix "Q." I know how to generate a matrix using proc when its dimension depends on one parameter. However, here, the dimension of matrix "Q" depends on two parameters, M1 and M2. So, I am a little bit confused about how to adjust them in proc. Please see the attached file and share your useful ideas. 

help.mw

Thanks in advance

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

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