Here's a partially geometric solution. For brevity let's write $x$, $y$, $z$ for the distances $SP$, $SQ$, $SR$, and $b$, $a$, $c$ for the corresponding side lengths, as in the following diagram. Furthermore, let $h$ be the length of altitude $CH$, parallel to the line $z$:

We want to find an $S$ that maximizes $xyz$.
From calculus we know that if $S$ is at a local extremum at we move it an infinitesimal distance in any direction, $xyz$ must not change (to first order). In particular, this must be true if we imagine moving $S$ a small distance $\varepsilon$ to the left (away from point $F$ in the diagram). When we do so, $x$ increases and $y$ decreases, but since we're moving $S$ horizontally, $z$ stays the same.
Because triangles $HAC$ and $PFS$ are similar we can compute the increase in $x$ as $\Delta x = \varepsilon h/b$. Mutatis mutandis we find $\Delta y = -\varepsilon h/a$ (with the opposite sign because an increase in $x$ corresponds to a decreas in $y$).
A bit of simple algebra now gives $ \Delta\, xyz = (x+\Delta x)(y+\Delta y)z - xyz = \frac{\varepsilon hyz}b - \frac{\varepsilon hxz}{a} + \text{terms involving }\varepsilon^2$
Because $\varepsilon$ is small we can ignore the $\varepsilon^2$ term. Since $\Delta xyz$ must be $0$ we then have $ \frac{\varepsilon hyz}b = \frac{\varepsilon hxz}a$ Canceling common factors (all nonzero, since $z>0$ inside the triangle) and cross-multiplying gives us $ \tag{1} ya = xb $ which we recognize as the equation of a line in the $x,y$ coordinate system. (This coordinate system is not rectangular, but equations of lines look the same in skew coordinates). If we can find two different solutions to $(1)$, the line connecting them will be a locus that every local extremum of $xyz$ inside the triangle must lie on.
The corner $C$ is is clearly a solution to $(1)$ because there we have $x=y=0$. Where can we find another one? The midpoint between $A$ and $B$ must be one, because there both sides of $(1)$ happen to be the area of triangle $ABC$.
So every local extremum of $xyz$ must lie on the median from $C$, and repeating the argument with another orientation we can see that it must lie on the other medians too. So the intersection of the medians is the only place in the inside of a triangle that can be a local extremum.
Since there quite obviously must be a global maximum somewhere in the triangle, the intersection of the medians has to be it.