@Markiyan Hirnyk  I do not understand the meaning of your comment.

@Preben Alsholm  Thank you. Very useful comment!

@Carl Love  Can we write this using  InertForm  package:

add(convert(i, symbol), i = 1 .. 20) = `+`($ 1 .. 20);

@Markiyan Hirnyk   This actually works. In fact, the result for a=1 (formula (1))  is the same as yours.

g := abs((x+4)/(x*a+2)) > 1/x:

sol := solve(g, [a, x]):

select(t -> is(eval(t[1], a = 1)), sol);

eval(map(t -> subsop(1 = NULL, t), %), a = 1);

@Preben Alsholm   If you write

restart;
n:=5: #Example
for i to n do B[i]:=LinearAlgebra:-RandomMatrix(n) end do; #Example B's.
A:=LinearAlgebra:-IdentityMatrix(n);
for i from 1 to n do A:=A.B[i] end do:
A;

then get the same result. 

Addition.  The same result can be gotten shorter:

`.`(seq(B[i], i=1..5));

Another simple way to remove the contradiction - use  unapply  instead  arrow to define a function:

f:=unapply(int(x - t, t=0 .. x), x);

f(t);

                  f := x -> 1/2*x^2

                       1/2*t^2

@sarra   In order that the matrix was square, the number of equations should equal the number of unknowns. You have 16 unknowns but only four equations.

What do you mean by a transposition?  Transposition (accoding to Wiki) a permutation which exchanges two elements and keeps all others fixed. In this sense, for your example is enough 2 transpositions: [w,x,y,z]->[y,x,w,z]->[y,x,z,w]

evalf(Int(abs(exp((2*Pi*I)*x)+exp((2*Pi*I)*y)), [x = 0 .. 1, y = 0 .. 1]));

identify(%);

                                                                             4/Pi

@Carl Love 

A := Matrix([[1, 2, 3], [3, 1, 2], [2, 3, 1]]):
B := Matrix([[1, 2, 3], [3, 1, 2], [2, 3, 1]]):

is(A=B);

                             false

@Carl Love   I do not understand why we need to override the matrices. Maple tests  A  and  A1  as the identical matrices:

restart;

with(LinearAlgebra):

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

A1 := Matrix(3, shape = Circulant[[1, 2, 3]]);

Equal(A, A1);

@Muhammad Usman   

applyrule(exp(x::anything)^y::anything=exp(x*y), expand(exp(2*a+b)));

@Markiyan Hirnyk   I do not know the causes of your problems. Here's a screenshot of the playing this animation in Power Point. Below is visible panel to control animation:

@Markiyan Hirnyk   I do not understand what's the problem? Click on the picture  and the bottom panel appears to start the animation.

Here is a sample ppt-file with the animation of the tangent line as the limiting position of the secant. Animation is made in Maple (M 12 classic), saved as a gif-file and inserted to Power Point in the above method.

Animation_tangent_line.pptx