The system of ODEs for the characteristic lines is
\frac{dx}{ds}(r,s)=x^2
\frac{dt}{ds}(r,s)=1
\frac{du}{ds}(r,s)=t.
The first ODE can be integrated by separation of variables:
\frac{dx}{ds}(r,s)=x^2
∫\frac{dx}{x^2}=∫ds
-\frac{1}{x(r,s)}=s+g_1(r)
x(r,s)=-\frac{1}{s+g_1(r)}.
The second ODE is straightforward:
\frac{dt}{ds}(r,s)=1
t=s+g_2(r).
The third ODE can be solved using the value of t just computed:
\frac{du}{ds}(r,s)=t=s+g_2(r)
u(r,s)=(s^2)/2+s*g_2(r)+g_3(r).
Note that the constants of integration are functions of r since we are integrating w.r.t. s. The general solution is:
x(r,s)=-\frac{1}{s+g_1(r)}.
t=s+g_2(r)
u(r,s)=(s^2)/2+s*g_2(r)+g_3(r).
The initial conditions are:
x(r,0)=r
t(r,0)=0
u(r,0)=f(r)
which we substitute into the general solution
x(r,0)=-\frac{1}{g_1(r)}=r
t(r,0)=g_2(r)=0
u(r,0)=g_3(r)=f(r)
getting the value for the constants of integration
g_1(r)=-\frac{1}{r}
g_2(r)=0
g_3(r)=f(r).
that we can substitute into u to get it as a function of s,r:
u(r,s)=s^2/2+f(r).
In order to get u as a function of x,t, we need to invert x(r,s),t(r,s) first
t=s
x=-\frac{1}{t-1/r}=\frac{r}{1-rt}
1-rt=\frac{r}{x}
1=r(t+\frac{1}{x})
r=\frac{x}{1+tx}
and then substitute them into the above expression for u:
u(x,t)=t^2/2+f(\frac{x}{1+tx}).