Yes, you are right. Maple in general weakly solves the examples with parameters. Maple got the solution only for general case, when  m <> -1 / 2 .

Here is the complete solution found with Mathematica (in output used notations for logical  and  and  or ):

Reduce [{-x3 * x5 /m - x3 * x5 - 2 * x1 * x5 + 2 * x1 == 0,
x3 * x5 /m - 2 * x2 * x5 + 2 * x2 ==
0,-х3 + 2 * x5 - x3 * x5 + 1 - 3 * x1 - 3 * x2 + x4 ==
0, -2 * x4 - 2 * x4 * x5 == 0,-x3 * x5 ^ 2 / (x1 * m) - 4 * х5 ^ 2 + 4 * x5 ==
0}, {x1, x2, x3, x4, x5}] / / TraditionalForm

It is strange that in yourcritical lines parameter  m  is missing. Also, when writing functions you missed the multiplication sign. I got the following result:

f := (x1, x2, x3, x4, x5)->-x3*x5/m-x3*x5-2*x2*x5+2*x1:

g := (x1, x2, x3, x4, x5)->x3*x5/m-2*x2*x5+2*x2:

h := (x1, x2, x3, x4, x5)->-x3+2*x5-x3*x5+1-3*x1-3*x2+x4:

k := (x1, x2, x3, x4, x5)->-2*x4-2*x4*x5:

l := (x1, x2, x3, x4, x5)->-x3*x5^2/(x1*m)-4*x5^2+4*x5:

solve({f(x1, x2, x3, x4, x5) = 0, g(x1, x2, x3, x4, x5) = 0, h(x1, x2, x3, x4, x5) = 0, k(x1, x2, x3, x4, x5) = 0, l(x1, x2, x3, x4, x5) = 0}, {x1, x2, x3, x4, x5});

L := [[1,2,3,2,3],[4,5,6,2,3],[7,8,9,2,3]]:

map(`[]`@f@op, L);

                         [[f(1,2,3,2,3)], [f(4,5,6,2,3)], [f(7,8,9,2,3)]]

More specific example:

L := [[1,2,3,2,3],[4,5,6,2,3],[7,8,9,2,3]]:

f:=(v,w,x,y,z)->`+`(v,w,x,y,z):

map(`[]`@f@op, L);

                                            [[11], [20], [29]]

restart; with(plots):

f := x->-x^3+4*x^2-3*x+1:

A := animate(spacecurve, [[x, f(x)*cos(phi), f(x)*sin(phi)], x = 0 .. 3, color = red, thickness = 3], phi = 0 .. 2*Pi, frames = 100):

B := animate(plot3d, [[x, f(x)*cos(phi), f(x)*sin(phi)], x = 0 .. 3, phi = 0 .. s, style = wireframe], s = 0 .. 2*Pi, frames = 100):

L1 := animate(spacecurve, [[0, t*cos(phi), t*sin(phi)], t = 0 .. 1, color = red, thickness = 3], phi = 0 .. 2*Pi, frames = 100):

L2 := animate(spacecurve, [[3, t*cos(phi), t*sin(phi)], t = 0 .. 1, color = red, thickness = 3], phi = 0 .. 2*Pi, frames = 100):

S1 := animate(plot3d, [[0, t*cos(phi), t*sin(phi)], t = 0 .. 1, phi = 0 .. s, style = wireframe], s = 0 .. 2*Pi, frames = 100):

S2 := animate(plot3d, [[3, t*cos(phi), t*sin(phi)], t = 0 .. 1, phi = 0 .. s, style = wireframe], s = 0 .. 2*Pi, frames = 100):

display(A, B, L1, L2, S1, S2, axes = normal, orientation = [50, 60], labels = [x, y, z]);

I removed the square brackets where they are not needed, and replaced them  by the parentheses (where required). Also left only two equations: eq1 and eq3.

restart:

with(plots): 

mb:=765 : mp:=587 : Ib:=76.3*10^3 : Ip:=7.3*10^3 : l:=0.92 : d:=10: F:=0.042: omega:=0.56 :

eq1:=(mb+mp)*diff(x(t),t$2)+mp*l*cos(alpha(t))*diff(alpha(t),t$2)+mp*l*diff(alpha(t),t)^2*sin(alpha(t))-F*sin(omega*t)=0;

eq3:=mp*(d+l*cos(alpha(t)))*diff(x(t),t$2)+(Ip+mp*l^2+mp*d*l*cos(alpha(t)))*diff(alpha(t),t$2)-mp*sin(alpha(t))*(-l*d*alpha(t)^2)+mp*l*9.81*sin(alpha(t))=0;

CI:= x(0)=0,alpha(0)=0,D(x)(0)=0,D(alpha)(0)=0;

solution:=dsolve([eq1, eq3, CI],numeric);

odeplot(solution,[[t,x(t)], [t,alpha(t)]], t=0..10, thickness=2);

Herd size that satisfies all the constraints, very large, and Maple can count it only approximately. Here are calculations using your link above:

w := 300426607914281713365*sqrt(609)+84129507677858393258*sqrt(7766):

n := evalf((w^4658-1/w^4658)^2/(4657*79072));

L := [7460514, 7358060, 10366482, 4149387, 4893246, 3515820, 7206360, 5439213]:

N := `+`(op(L));

N*n;

Brian, you can solve your problem by  solve  command, if  d  will be equal the least common multiple of all denominators.

restart;

unprotect(D):

L:=[B,D,W,Y,b,d,w,y]:

solve({W=5/6*B+Y, B=9/20*D+Y, D=13/42*W+Y, w=7/12*(B+b),

  b=9/20*(D+d), d=11/30*(Y+y), y=13/42*(W+w)});

assign(%):

d := ilcm(op(map(denom, L)));

L;

You can simply write:

all_plots:=[plot(sin(10*x+0.2),x=0..1, thickness=10, color=blue), plot(1-sin(10*x),x=0..1, thickness=10, color=red), plot(sin(10*x),x=0..1, thickness=10, color=green)]: 

plots:-display( all_plots );

 Done in M 12 classic.

Idea of Markiyan in above link was used.

f := piecewise(0 <= x and x <= 2, x^2, 2<x and x<3, undefined, 3 <= x and x <= 4, x^2-3, 4<x and x<5, undefined, 5 <= x and x <= 6, x^2-10):

plot(f, x = 0 .. 6, tickmarks = [default, [seq(i = i, i = 1 .. 4), 5 = "...", seq(j = j+3, j = 6 .. 13), 14 = "...", seq(k = k+10, k = 15 .. 26)]], thickness=2, discont = true, title = "Plot of  x^2");

dsolve(-(diff(lambda(s), s))-2*(diff(lambda(s), s, s))-(diff(lambda(s), s, s, s)) = 0):

subs({_C1=C[1], _C2=C[2], _C3=C[3]}, %);

The first result is also incorrect. Should be:  undefined if  -100<n<-10,  -ln|n+10|+ln|n+100| otherwise, if we consider  ln(0)=-infinity

Unfortunately, Maple weakly solves examples with parameters. Here are another two examples where the answer is incorrect:

solve(a*x = b, x);

solve(sqrt(x-a) = x, x);

Mathematica correctly solves both of these examples:

Vectors can be constructed by  plots[arrow]  command.

Example:

A := Matrix(3, 3, [1, 2, 3, 5, 1, 6, 6, 3, 9]);

V1, V2, V3:=A[..,1], A[..,2], A[..,3];

P:=plots[arrow]({V1, V2, V3}, color=red, width=[0.1, relative=false],scaling=constrained, axes=normal, orientation=[45,75]):

T:=plots[textplot3d]([[1, 5, 6, "V1"], [2, 1, 3, "V2"], [3, 6, 9, "V3"]], color=black):

plots[display](P, T, view=-1..10);

Somewhat easier to just use plot command:

f:=x->x^3+3*x^2+1:

 plot([seq([x, f(x)], x=-10..10)]);

Your procedure works correctly.

Example:

Digits := 20:

hw1(g, 0., 9):

seq(x[i], i = 0 .. 10);

 Explanation:  by default, the output of a procedure call is the last calculation.

Your spheres can be plotted by the code in single line without loops. For clarity, I reduced the number of spheres to 5:

plots[implicitplot3d]([seq(x^2+y^2+z^2 = (2*i)^2, i = 1 .. 5)], x = -10 .. 10, y = -10 .. 0, z = -10 .. 10, style = surface, numpoints = 10000, view = [-11 .. 11, -11 .. 7, -11 .. 11], scaling = constrained, axes = normal, orientation = [50, 75]);