Consider $C^\infty$ functions. For functions of one variable, we know that in general, at $x=c$ if the lowest order nonzero derivative is an even order say $2k$, then $f(c)$ is an extremum, depending on the sign of $f^{(2k)}(c)$. For example, if $f'(c)=f''(c)=f'''(c)=0$ and $f''''(c)>0$, then it is a minimum; if $f'(c)=f''(c)=0$, $f'''(c)\neq 0$, then it is not a local extremum.
Now for multivariable functions. We know the theory of Hessian matrix and eigenvalues. What if some of the second derivatives are zero? For example, consider $f(x,y)$ at $(a,b)$: $f_x=f_y=f_{xx}=f_{xy}=0$. Apparently if $f_{xxx}\neq 0$ the function cannot have a local extremum there, as the trend along the $x$-axis is monotone. Suppose $f_{xxxx},f_{yy}$ do not vanish. Should we examine the eigenvalues of \begin{equation} H=\begin{bmatrix} f_{xxxx} & f_{xxy}\\ f_{xxy} & f_{yy}\end{bmatrix}? \end{equation} What is the general theory about it? Any reference on this? thanks!