Skalieren entlang einer Achse musst Du definieren wie das gehen soll. Deine Achse ist 3D (x,y,z), Deine Matrix ist die einer Skalierung 2D in homogenen Koordinaten:
Skalierung (Vergrößerung oder Verkleinerung)
Beim Skalieren eines Objektes um den Faktor s um den Ursprung (0,0) wird ein Punkt (x,y) auf (x ́,y ́) = (s·x, s·y) abgebildet.
3D-Skalierung in homogenen Koordinaten
\(\left(\begin{array}{r}x'\\y'\\z'\\1\\\end{array}\right)= \left(\begin{array}{rrrr}s_x&0&0&0\\0&s_y&0&0\\0&0&s_z&0\\0&0&0&1\\\end{array}\right) \cdot \left(\begin{array}{r}x\\y\\z\\1\\\end{array}\right) \)