sol := dsolve(sys union ics, {x(t), y(t)}, numeric, method = rkf45, range = 0 .. 20000, maxfun = 10000000, output = listprocedure, abserr = 1e-12, relerr = 1e-12, stepsize = 0.001);

Download maple_primes_post_partiel_diff.mw

with(Student:-Calculus1)
infolevel[Student[Calculus1]] := 1

To get a idea , you can use Hint  from Rule package
If you are uncertain about which rule to apply next to a problem, ask for a hint.
Hint(Limit(-tanh(sqrt(2)*(a*x+b)),x=0));
Check one-sided limits separately
 

The output from the Hint command can always be used as the index to the Rule command.
Rule[[chain]](Diff(sin(x^2), x)); (see example Maple ) , how it goes further ?
 

I assume ax+ b that all variables are real numbers ?

@GFY 
This simplification for now
 

variable := 15

Leaf count value of the original expression: 6879

Leaf count value of method 15 (simplify(e, symbolic)  ,  e is expression                                                  ): 2967

you can check both expressions if they are equal 
example :

restart;
                      "maple.ini in users"

e1 := -sqrt(-(exp(-2 + 2*x) - 2)*exp(-2 + 2*x))/(exp(-2 + 2*x) - 2);
e2 := 1/sqrt(2*exp(-2*x)*exp(2) - 1);
                                                   (1/2)
               (-(exp(-2 + 2 x) - 2) exp(-2 + 2 x))     
       e1 := - -----------------------------------------
                           exp(-2 + 2 x) - 2            

                                  1              
              e2 := -----------------------------
                                            (1/2)
                    (2 exp(-2 x) exp(2) - 1)     

# Definieer de expressies
simplified_e1 := simplify(e1);
simplified_e2 := simplify(e2);


                                                         (1/2)
                     (-(exp(-2 + 2 x) - 2) exp(-2 + 2 x))     
  simplified_e1 := - -----------------------------------------
                                 exp(-2 + 2 x) - 2            

                                       1             
          simplified_e2 := --------------------------
                                                (1/2)
                           (2 exp(-2 x + 2) - 1)     

# Stel de vergelijking op
eq := simplified_e1 = simplified_e2;


                                                  (1/2)   
              (-(exp(-2 + 2 x) - 2) exp(-2 + 2 x))        
      eq := - ----------------------------------------- = 
                          exp(-2 + 2 x) - 2               

                    1             
        --------------------------
                             (1/2)
        (2 exp(-2 x + 2) - 1)     


# Vermenigvuldig beide zijden van de vergelijking met de noemer van de rechterkant
left_side := lhs(eq):
right_side := rhs(eq):
new_eq := left_side*(2*exp(-2*x + 2) - 1) = right_side*(2*exp(-2*x + 2) - 1);


new_eq := - 

                                      (1/2)                         
  (-(exp(-2 + 2 x) - 2) exp(-2 + 2 x))      (2 exp(-2 x + 2) - 1)   
  --------------------------------------------------------------- = 
                         exp(-2 + 2 x) - 2                          

                       (1/2)
  (2 exp(-2 x + 2) - 1)     


# Vereenvoudig de nieuwe vergelijking
simplified_new_eq := simplify(new_eq);


simplified_new_eq := 

                                                             (1/2)   
  (-2 exp(-2 x + 2) + 1) (-(exp(-2 + 2 x) - 2) exp(-2 + 2 x))        
  ---------------------------------------------------------------- = 
                         exp(-2 + 2 x) - 2                           

                       (1/2)
  (2 exp(-2 x + 2) - 1)     


solutions := solve(simplified_new_eq, x);


                         solutions := x

solutions;
                               x


==========================================================================
                      
another example 
# Definieer de expressies
expr1 := exp(x) + exp(-x) assuming(x::complex);
expr2 := 2*cosh(x) assuming(x::complex);

# Stel de vergelijking op
eq := expr1 = expr2:

# Vereenvoudig de vergelijking
simplified_eq := simplify(eq):

# Los de vergelijking op voor x (om te controleren of ze gelijk zijn)
solutions := solve(simplified_eq, x):

# Output de resultaten
simplified_eq;
solutions;

                   expr1 := exp(x) + exp(-x)

                       expr2 := 2 cosh(x)

                  exp(x) + exp(-x) = 2 cosh(x)

                               x

is(expr1 = expr2) assuming(x::complex);
                              true

..as you can see here :
 Each trigonometric function in terms of each of the other five List of trigonometric identities - Wikipedia

(all this mathematical knowledge here on wiki is included in Maple i think) 

Knowing the unit circle and all the graphs of sin, cos , tan for deducing some formulas and for goniometric equations.

Manual adding table in Maple and then add plots into the cells is probably not useful for you?

Got this from someone else here on forum

Download complex_log_plot.mw

f := x -> x^sin(x):

Inflpts := [fsolve(D(D(f))(x), x=0..16, maxsols=6)];

Q := map(p->[p,f(p)], Inflpts):

T := seq(plot(D(f)(p)*(x-p)+f(p), x=p-1..p+1, color=red), p=Inflpts):

plots:-display(plot(f, 0.0..16.0, color=black), T,

        plots:-pointplot(Q, symbolsize=10, symbol=solidcircle, color=blue),

        map(p->plots:-textplot(evalf[4]([p[1]-sign(D(f)(p[1]))*2/3,p[2]+1,p]),

                               font=[Times,8]),Q), size=[800,400]);

a example fo infliction points, can be adjusted for find min/max  :to give a idea

with(Student[Basics]);
  [ExpandSteps, FactorSteps, LinearSolveSteps, LongDivision, 

    OutputStepsRecord, PracticeSheet, SolveSteps]

-----------------------------------------------------------------

FactorSteps(Y^2*x^3-x^3);  
How about for complex numbers ?

Given a polynome in C : z^4-2.z^3+ 3.z^2-2.z+2