## 1540 Reputation

19 years, 144 days

## Corrupted file....

Hope this helps (edited):

restart;
A[1]:=alpha=1-(1/2)/(1-(RootOf(16*_Z*(_Z*(2*_Z*(_Z*(8*_Z*(_Z*(_Z*(_Z*(32*_Z*(8*_Z-33)+1513)-812)-13)+267)-1469)-330)+811)+279)+345,index=2)-1/2)**2);
A[2]:=expr=(1+alpha)*sqrt(1-alpha**2)+(3+4*alpha)/12*sqrt(3-4*alpha**2)+2*(1+alpha)/3*sqrt(2*(1+alpha)*(1-2*alpha))+(1+2*alpha)/6*sqrt(2*((1-alpha)**2-3*alpha**2));
alias(beta=indets(A[1],RootOf)[1]);
F := lambda^2+(8159503/855552)*beta-(236960/1671)*beta^9+(2085340/5013)*beta^8-(13276063/40104)*beta^7+(414451/20052)*beta^6+(10073755/160416)*beta^5-(13077487/160416)*beta^4+(1471997/213888)*beta^3+(5380211/160416)*beta^2+1770613/2566656;
evalf(subs(lambda = expr, evalf(subs(A[1], A[2]), 40), F), 40);
F:=(evala@Norm)(F);
evalf(subs(lambda = expr, evalf(subs(A[1], A[2]), 40), F), 40);

Also, the file you've submitted appears corrupted, so says Maple 2017.
Sorry, my previous posting was based on a typo. This one is corrected.

## Just the plot and a few solutions:...

restart;
f1 := x1 - x1*sin(x1 + 5*x2) - x2*cos(5*x1 - x2);
f2 := x2 - x2*sin(5*x1 - 3*x2) + x1*cos(3*x1 + 5*x2);
ds := {d1 = 3/13, d2 = 2/13};
rgs := subs(ds, [seq(seq([x1 = 2*Pi*i*d1 .. 2*Pi*(i + 1)*d1, x2 = 2*Pi*j*d2 .. 2*Pi*(j + 1)*d2], i=(floor(subs(ds, -1/d1)) .. ceil(subs(ds, 1/d1))-1)), j = (floor(subs(ds, -1/d2)) .. ceil(subs(ds, 1/d2)) - 1))]):
t0 := time():
Digits := 32:
S := NULL:
for r in rgs do
#print(r);
S0 := NULL;
for i from 1 to 20 do
q := fsolve({f1, f2}, {x1, x2}, avoid = {S0}, op(r)):
if type(q, specfunc(anything, fsolve)) then break else S0 := S0, q fi; od:
S := S, S0;
od:
t1 := time():
print(t1 - t0);
print(nops([S]));
C := plots[pointplot](map(subs, [S], [x1, x2]), symbol = solidcircle, color = violet):
Digits := 12:
plots[display](C, seq(plots[implicitplot](L[1], x1 = -2*Pi .. 2*Pi, x2 = -2*Pi .. 2*Pi, grid = [300, 300], color = L[2]), L=[[f1, red], [f2, blue]]), size = [1200, 1200], scaling = constrained);

Also replace color specification in plot C: violet for COLOUR(RGB, .723921568, .673725492, .723921568). It looks better.

Notice that x1=0 and x2=0 cases are not in the solution set. That is probably because these are boundary cases in the relevant sub ranges. Easily solved using solve.

Moderator edit: Here's the plot, showing 190 solutions.

## Use map....

map2(map, convert, A, D);

## Two errors in the code corrected....

with(PolynomialIdeals);
R := PolynomialRing([x, y, z, u, v], 'lex');
P := PolynomialIdeal([x^7*v, y^8*v, z^9*v, u*v]);
decomposition := PrimaryDecomposition(P);
htP := nops([decomposition]);

## Complex vs real...

f := x -> x^3:
g := x -> x^(1/3):

#For complex domain these answers are correct:
f(g(x));
g(f(x));

#For real domain these answers are correct:
A := [x < 0, x >= 0]:
f(g(x)) assuming x, real;
proc(E) global A; local a; local S := proc(E, a) a, simplify(E) assuming a end; 'piecewise'(seq(S(E, a), a = A)); end(g(f(x)));
F2 := %:

plot([Re, Im](F2), x = -2 .. 2, color = [red, blue], axes = framed);
plot([Re, Im](g(f(x))), x = -2 .. 2, color = [red, blue], axes = framed);

## After simplification:...

After simplifying for real t, res appears as:

-abs(t)^3*csgn(2*sqrt(2)+3*ln(t))+t^3;

## Trig simplification....

If You know the final form then you can demonstrate with this:

R := RootOf(4*_Z^2+(4*RootOf(60*_Z^3-60*_Z^2+15*_Z-1)-4)*_Z+4*RootOf(60*_Z^3-60*_Z^2+15*_Z-1)^2-4*RootOf(60*_Z^3-60*_Z^2+15*_Z-1)+1);
final := 1/3+(1/3)*cos((1/3)*arctan(3/4));

Finalform := (1/3)*cos(alpha)+1/3;
(`@`(combine, factors, evala, Norm))(Finalform-R)[2, 1, 1];
subs(solve(%, {alpha}), Finalform);
evalf(Testzero(evalf(%-final)), 20);

#or more precisely:

Finalform := (1/3)*cos(arctan(alpha)/3)+1/3;
(`@`(combine, factors, evala, Norm))(Finalform-R)[2, 1, 1];
subs(solve(%, {alpha}), Finalform);
evalf(Testzero(evalf(%-final)), 20);

## Use assume, not asumme....

Also additional assumptions should use command additionally not assume. Using assume again on the same variable reinitializes all assumptions on it (and related variables I believe).

## For that form of differential problem us...

ode := diff(y(x), x) = sin(y(x))+1;
ic := y(0) = Pi;
Order := 40;
sol := dsolve([ode, ic], y(x), series);
ser1 := gfun[seriestorec](subs(sol, y(x)), a(n));
diffeq1 := gfun[rectodiffeq](ser1[1], a(n), y(x));
dsol1 := dsolve(diffeq1, y(x));

## Maple types....

This should help:

indets(cS(a)(b)(c), specfunc(anything, cS));

Maple does not resolve types in the 0th operand.

## Do not use unapply....

Replace
Int(unapply(eq7,t),-25..25,method=_d01ajc)
with
Int(eq7,t=-25..25,method=_d01ajc)

## ?ifactors...

Instead of ifactor use ifactors. Much easier to work with.

## A minor improvement....

y = -1/2 .. 3, z = -3 .. 3, x = 0 .. 1;
u = 2*x*y, v = 2*x*z;
u = (minimize .. maximize)(2*x*y, x = 0 .. 1, y = -1/2 .. 3);
v = (minimize .. maximize)(2*x*z, x = 0 .. 1, z = -3 .. 3);
u = -1 .. 6, v = -6 .. 6, x = 0 .. 1;
uvEQ := collect(combine(subs(solve({u = 2*x*y, v =  2*x*z}, {y, z}), expand([eq1, eq2]*x^2*300))), [exp, sin, cos]);

Digits := 16;
PLT1:=plots[intersectplot](
op(map(E->map(evalf, surface(E, u = -1 .. 6, v = -6 .. 6, x = 0 .. 1)), uvEQ)),
maxlev = 5, maxtet=100000, grid = [31, 61, 11], thickness = 3, transparency = 0.3):
print(PLT1);

subs(solve({u = 2*x*y, v = 2*x*z}, {y, z}), [y, z, x]);
F1 := unapply(%, [u, v, x]);
plots[display](plottools[transform](F1)(PLT1), view=[ -1/2 .. 3,  -3 .. 3,  0 .. 1]);

## See ?plots[intersectplot]...

If you'd like a plot of the solution you can use:

Digits := 16;
plots[intersectplot](
map(evalf, surface(eq1, y = -1/2 .. 3, z = -3 .. 3, x = 0 .. 1)),
map(evalf, surface(eq2, y = -1/2 .. 3, z = -3 .. 3, x = 0 .. 1)),
maxlev = 5, grid = [31, 61, 11], thickness = 3, transparency = 0.3);

## ArrayTools....

You could use this, if you mean exact match.

ArrayTools:-IsZero(A) and LinearAlgebra:-Dimensions(A)=(4,4);

Also see ?ArrayTools:-IsEqual

 1 2 3 4 5 6 7 Last Page 1 of 21
﻿