Note: For $y^2 \equiv x^3 + ax + b$, you want $$m\equiv(3x^2+a)*(2y)^{-1}$$ which is a shortcut that only makes sense once you understand how elliptic curves work.
"Why does it say y=1 and 4 for x=1?"
Like you said, plug in & solve: $$y^2 \equiv 1^3 + 2(1) + 3 \equiv 6 \equiv 1\space(mod\space5)$$
Ask yourself "Is there a square such that it is $\equiv 1 \space(mod\space 5)$?"
$$1^2\equiv1 \space \space and\space\space 4^2\equiv16\equiv1\space (mod \space5).$$
This means (1,1) and (1,4) are points on the discrete version (mod 5) of the curve.
If we weren't working (mod 5), plugging in 1 would give us (1,2) and (1,-2), which are points on the continuous version of the elliptic curve.
Below on the left is the continuous version of an arbitrary elliptic curve; on the right is what happens when you start working mod a number.
Example: Find 2P on the curve $E: y^2 \equiv x^3 + 2x + 3 \space(mod\space5)$ where P = (1,3)
2P = (1,3) + (1,3). When adding two points on an elliptic curve, there are two cases:
Case 1: The points are not the same: Find the line that includes both of the points on the curve. This line will have a third point on the curve. Find that point and reflect about the x-axis to have the final point.
Case 2: The points are the same: Find the line tangent to our point on the curve. This line will intersect the curve at one more point. Find that point and reflect about the x-axis to have the final point.
Now what is the line tangent to (1,3) on E?
Use implicit differentiation (simple calculus):
$$dy/dx = (3x^2 + 2)/(2y)$$.
Now at our point (1,3): $$dy/dx = 5/6 = 5*6^{-1} (mod \space5)$$
It is easy to see (and find) that $6^{-1} \equiv 1 \space(mod\space 5)$
So m = 5 is the slope for our line tangent to (1,3): $$L: y-3 = 5(x - 1)$$
So $$y = 5x - 2.$$
Now where else do E and L intersect?
We have
Plug L into E: $$(5x-2)^2 \equiv x^3 + 2x + 3$$
Expand: $$25x^2 - 20x + 4 \equiv x^3 + 2x + 3$$
There is a useful algebraic property for cubic polynomials that is fun and not hard to prove: The summation of the polynomial's roots equal the negative of the $x^2$ term's coefficient. That is,
$$\sum roots = -(coeff.\space of\space x^2\space term)$$
So in rearranging the polynomial further, we care only about the $x^2$ terms:
$$0 \equiv x^3 - 25x^2 ...$$
We already added (1,3) to itself, so it must be that 1 is a root of multiplicity 2.
$$1 + 1 + r_3 = -(-25) = 25$$ implies the last root is 23. Recall we are working (mod 5) though, so our last root is 23 $\equiv 3$
Now plug $r_3$ into our linear equation: $$y = 5(3) - 2 = 13 \equiv 3\space(mod\space5)$$
Reflect
Lastly, reflect our point (3,3) about the x-axis to get (3,-3). Recall again we're working (mod 5), so add the modulus to -3 until is it in the range $[0,4]$.
$\therefore$ 2P = (3,2)