restart;
ls:=["a", "b", "c", "d", "e"]:
selectremove(t->ListTools:-Search(t,ls)::odd, ls);

The easiest way to solve the problem is to show your graphs in the  x  and  y ranges (using the  view  option), where the differences are clearly visible:

restart:
with(plots):
f1:=x->0.95*x;
f2:=x->0.99*x;
f3:=x->0.991*x;
p1:=plot(f1(x),x=0..10,color=red):
p2:=plot(f2(x),x=0..10,color=blue):
p3:=plot(f3(x),x=0..10,color=green):
P1:=display(p1,p2,p3, caption="Plot 1", captionfont=[times,16]); 
P2:=display(p1,p2,p3, view=[9.5..10,9.2..9.7],caption="Plot 2", captionfont=[times,16]);
display(<P1 | P2>);

Download plotting1.mw

Edit.

Maple finds the answer through elliptic functions. Unfortunately, the result is quite cumbersome:

restart;
y:=Int((k/2*x^2+x+c)/sqrt(lambda^2*(1+k*x)^4 - (k/2*x^2+x+c)^2), x);
value(y):
simplify(%);

Download Int.mw

Edit. In my opinion, the value of this overall result is rather dubious. With some parameters and variable  x  values, you can get  0 or negative numbers under the root. It is much more reliable to first set the values of the parameters, and only then calculate the integral.

Use  add  instead of  sum  when summing a specific finite number of terms. The  sum  command is used if at least one of the summation limits is a symbol or the infinity:

restart;
seq(add(x^r*(add((-1)^s*binomial(r, s)*pochhammer((a+s+1)*d, n), s = 0 .. r))/factorial(r), r = 0 .. n), n = 0 .. 4);

It is easy to prove that your second graph is empty. Probably you rewrote something wrong. I modified the second equation a bit to make the graph look like your second graph in the picture:

Edit. Let us prove that the equation  abs(z-I) - abs(z+1)  =  2  has no solutions for complex numbers  z=x+I*y . This is most easily done using the fact that   abs(z1 - z2)  is the distance between the points  z1  and  z2 . Consider 3 points  z1=I ,  z2=-1 and  z . It follows from the triangle inequality that  abs(z-z1) - abs(z-z2) < abs(z1-z2) = sqrt(2)<2

Download z.mw

If you need to do this more than once, then (as a development of vv's idea) a simple procedure  Sort  can be used to reduce typing:

restart;
S:=[seq(_C||k=c[k], k=0..10)]:
sol:=dsolve(diff(y(x),x) = x+y(x),y(x));

Sort:=s->sort(subs(S,s)):

Sort(sol); 

Use  piecewise  instead of  if .. fi  :

Download if_new.mw

I do not know the reasons for this behavior, but to immediately get the desired position of the constants use the  sort  command:

restart;
ode:=diff(diff(y(x),x),x)+8*diff(y(x),x)+25*y(x) = 1:
sol:=sort(dsolve(ode));

This code in Maple 2018.2

Here is another proof without complex numbers using the well-known rotation matrix by an angle. We take the unit vector  <1,0>  and rotate it first by the angle of  2*x , we get the vector  X . Rotating the same vector  <1,0>  sequentially, first by the angle  x , and then again by x , we get a vector  Y . Obviously  X=Y . At the same time, we proved another formula  cos(2*x)=cos(x)^2-sin(x)^2 .

Download proof.mw

You are probably interested in calculating the area of the region bounded by a red closed curve below (this is not an ellipse of course).
It's easy to do numerically:

Download int.mw

To solve your example, it is enough to know the simplest definitions related to the joint distribution of two random variables  X  and  Y . Let's introduce the notations:  f__XY  is the density function of the joint distribution of  X  and  Y ,  f__X  is the density function of  X ,  f__Y  is the density function of  Y ,  f__X_Y  is the density function of the conditional distribution of  X  given  Y ,   f__Y_X  is the density function of the conditional distribution of  Y  given  X .

Download XY.mw

You have chosen too narrow a range of values for the contours. If we increase it, we get more interesting pictures, for example:

restart;
randomize():
a, b, c_3, c_4, A_5:='rand(-5..5)'()$5;
A := a*(A_5)^3-b; B := 3*a*(A_5)^2-3*b/(A_5); C := (c_3)*(A_5)-(c_4);
u:=3*y^2+2*x^3*B/A+3*x^2*C/A;
plots:-contourplot(u, x=-2.25..2.25, y=-2.25..2.25, axes=boxed, contours=[seq(-10..10,1)],grid=[200,200], coloring=["Blue","Red"], size=[500,500]);

Below are 2 solutions to the problem. In the first one, we randomly set all the necessary parameters. In the second one, we use  Explore  command for this:

MaplePrime did not display the result of the last command. See the attached file for this.

Download contour_new.mw

Addition - another way:

T__e_np := sqrt(r^2+1)/(r^2+sqrt(r^2+1)*r+1);
numer(T__e_np)*sqrt(r^2+1)/expand(denom(T__e_np)*sqrt(r^2+1));
simplify(%);

Download Q20201208_new.mw

Maple simply doesn't know what it should return if  k  is a symbol. Therefore, we must provide for this in the body of the procedure. Below is a possible solution to the problem:
 

Download divide5_new.mw