My question regards what I assume would be called the general method of characteristics for third-order equations. I am looking for PDE for a function of two variables $x$ and $y$ having a characteristic curve in the form $$ y^2=x^3+\alpha x+\beta~~. $$ This type of curve is called an elliptic curve. $\alpha$ and $\beta$ are two constants and the elliptic discriminant is such that the curve is non-singular (obviously). I will make an example using the method of characteristics for first-order equations, and then I will say more precisely what I'm looking for once I have laid out the context. Given a function $u(x,y)$ and $$ \partial_x u+\alpha\,\partial_y u=0~~, $$ Following what I understand to be the normal algorithmic procedure for solving such PDEs, I will introduce $$ \dfrac{dy}{ds}=\alpha\qquad\qquad \dfrac{dx}{ds}=1\qquad\qquad \dfrac{dz}{ds}=0. $$ Integration yields $$ y(s)=as+c_1\qquad\qquad x(s)=s+c_2 \qquad\qquad z(s)=c_3~~. $$ Eliminating $s$, we find the characteristic curves for the PDE are given by the following lines in $\mathbb{R}^3$: $$ y-\alpha x=y_0\qquad\text{and}\qquad z=k~~. $$ Here, $y_0$ and $k$ are two constants. It follows that $z(x,y)$ is constant when $y-\alpha x$ is constant. Setting $u=z$, we find $u(x,y)=f(y-\alpha x)$ is a solution to the stated PDE.
Now, my question: I want to find a PDE for $u(x,y)$ whose characteristic curve in the plane of those variables (meaning the analogue of $y-\alpha x=y_0$) is an elliptic curve.