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