We begin with part (1). Let
\(I_f=\int_a^b f\) and
\(I_g=\int_a^b g\text{,}\) and fix
\(\varepsilon \gt 0\text{.}\) By the
Integrability Criterion, there are partitions
\(P_f\) and
\(P_g\) of
\([a,b]\) such that
\begin{equation*}
U(P_f,f)-L(P_f,f) \lt \frac{\varepsilon}{2},
\qquad
U(P_g,g)-L(P_g,g) \lt \frac{\varepsilon}{2}.
\end{equation*}
Let \(P=P_f \cup P_g\text{.}\) Since refinement decreases upper sums and increases lower sums,
\begin{equation*}
U(P,f)-L(P,f) \lt \frac{\varepsilon}{2},
\qquad
U(P,g)-L(P,g) \lt \frac{\varepsilon}{2}.
\end{equation*}
For each subinterval \(I_i\) of \(P\text{,}\) let \(m_i^f,M_i^f\) and \(m_i^g,M_i^g\) be the infima and suprema of \(f\) and \(g\) on \(I_i\text{.}\) Then
\begin{equation*}
m_i^f+m_i^g \le \inf_{x \in I_i}(f(x)+g(x)),
\qquad
\sup_{x \in I_i}(f(x)+g(x)) \le M_i^f+M_i^g.
\end{equation*}
Summing over the partition gives
\begin{equation*}
L(P,f)+L(P,g) \le L(P,f+g) \le U(P,f+g) \le U(P,f)+U(P,g).
\end{equation*}
Hence
\begin{equation*}
U(P,f+g)-L(P,f+g)
\le
\bigl(U(P,f)-L(P,f)\bigr)+\bigl(U(P,g)-L(P,g)\bigr)
\lt \varepsilon.
\end{equation*}
So \(f+g\) is integrable.
Also, because \(I_f\) lies between \(L(P,f)\) and \(U(P,f)\text{,}\) and similarly for \(I_g\text{,}\) we have
\begin{equation*}
I_f+I_g-\varepsilon
\lt
L(P,f)+L(P,g)
\le
\int_a^b (f+g)
\le
U(P,f)+U(P,g)
\lt
I_f+I_g+\varepsilon.
\end{equation*}
Since \(\varepsilon\) is arbitrary, \(\int_a^b (f+g)=I_f+I_g\text{.}\)
For scalar multiples, let \(P=\{x_0,\dots,x_n\}\) be any partition, and let \(m_i,M_i\) be the infima and suprema of \(f\) on \([x_{i-1},x_i]\text{.}\) If \(c \ge 0\text{,}\) then the infimum and supremum of \(cf\) on the same interval are \(cm_i\) and \(cM_i\text{.}\) If \(c \lt 0\text{,}\) they are \(cM_i\) and \(cm_i\text{.}\) In either case,
\begin{equation*}
U(P,cf)-L(P,cf)=|c|\bigl(U(P,f)-L(P,f)\bigr).
\end{equation*}
Hence \(cf\) is integrable. When \(c \ge 0\text{,}\)
\begin{equation*}
L(P,cf)=cL(P,f),
\qquad
U(P,cf)=cU(P,f),
\end{equation*}
and the same squeeze argument as above gives \(\int_a^b (cf)=c\int_a^b f\text{.}\) When \(c \lt 0\text{,}\) the same conclusion follows from the identities
\begin{equation*}
L(P,cf)=cU(P,f),
\qquad
U(P,cf)=cL(P,f).
\end{equation*}
For part (2), if \(f(x) \le g(x)\) on \([a,b]\text{,}\) then on every subinterval of every partition \(P\) the infimum and supremum of \(f\) are at most the corresponding infimum and supremum of \(g\text{.}\) Thus
\begin{equation*}
L(P,f) \le L(P,g)
\qquad\text{and}\qquad
U(P,f) \le U(P,g).
\end{equation*}
Since \(\int_a^b f \le U(P,f)\) and \(L(P,g) \le \int_a^b g\) for every partition \(P\text{,}\) it follows that \(\int_a^b f \le \int_a^b g\text{.}\)
For part (3), let
\(I=\int_a^b f\) and fix
\(\varepsilon \gt 0\text{.}\) By the
Integrability Criterion, there is a partition
\(P\) of
\([a,b]\) such that
\(U(P,f)-L(P,f) \lt \varepsilon/2\text{.}\) If
\(c \notin P\text{,}\) adjoin
\(c\) to obtain a refinement; the gap does not increase. Write
\(P_1\) and
\(P_2\) for the induced partitions of
\([a,c]\) and
\([c,b]\text{.}\) Then
\begin{equation*}
L(P,f)=L(P_1,f)+L(P_2,f),
\qquad
U(P,f)=U(P_1,f)+U(P_2,f),
\end{equation*}
so
\begin{equation*}
\bigl(U(P_1,f)-L(P_1,f)\bigr)+\bigl(U(P_2,f)-L(P_2,f)\bigr)
=
U(P,f)-L(P,f)
\lt
\frac{\varepsilon}{2}.
\end{equation*}
Since each summand is nonnegative, both are less than \(\varepsilon/2\text{.}\) Therefore the restrictions of \(f\) to \([a,c]\) and \([c,b]\) are integrable.
Let \(I_1=\int_a^c f\) and \(I_2=\int_c^b f\text{.}\) Because \(I_1\) and \(I_2\) lie between the corresponding lower and upper sums, we obtain
\begin{gather*}
I_1+I_2-\varepsilon
\lt
L(P_1,f)+L(P_2,f)
=
L(P,f)
\le
I,\\
I
\le
U(P,f)
=
U(P_1,f)+U(P_2,f)
\lt
I_1+I_2+\varepsilon.
\end{gather*}
Hence \(I=I_1+I_2\text{.}\)
Finally, for part (4), the hypothesis gives \(-M \le f(x) \le M\) on \([a,b]\text{.}\) By part (2),
\begin{equation*}
\int_a^b (-M) \le \int_a^b f \le \int_a^b M.
\end{equation*}
Since constant functions are integrable and
\begin{equation*}
\int_a^b (-M)=-M(b-a),
\qquad
\int_a^b M=M(b-a),
\end{equation*}
we conclude that
\begin{equation*}
-M(b-a) \le \int_a^b f \le M(b-a).
\end{equation*}
Therefore \(\left|\int_a^b f\right| \le M(b-a)=M|b-a|\text{.}\)