For example, he invented the Kruskal count, a magical effect that has been known to perplex professional magicians because it was based not on sleight of hand but on a mathematical phenomenon.
His Ph.D. dissertation, written under the direction of Richard Courant and Bernard Friedman at New York University, was on the topic "The Bridge Theorem For Minimal Surfaces".
In 1960, Kruskal discovered the full classical spacetime structure of the simplest type of black hole in general relativity.
Kruskal (in parallel with George Szekeres) discovered the maximal analytic continuation of the Schwarzschild solution, which he exhibited elegantly using what are now called Kruskal–Szekeres coordinates.
This led Kruskal to the astonishing discovery that the interior of the black hole looks like a "wormhole" connecting two identical, asymptotically flat universes.
In the 1970s, when the thermal nature of black hole physics was discovered, the wormhole property of the Schwarzschild solution turned out to be an important ingredient.
Kruskal's most widely known work was the discovery in the 1960s of the integrability of certain nonlinear partial differential equations involving functions of one spatial variable as well as time.
This work was partly motivated by the near-recurrence paradox that had been observed in a very early computer simulation[12] of a certain nonlinear lattice by Enrico Fermi, John Pasta, Stanislaw Ulam and Mary Tsingou at Los Alamos in 1955.
Those authors had observed long-time nearly recurrent behavior of a one-dimensional chain of anharmonic oscillators, in contrast to the rapid thermalization that had been expected.
Solitary wave phenomena had been a 19th-century mystery dating back to work by John Scott Russell who, in 1834, observed what we now call a soliton, propagating in a canal, and chased it on horseback.
This was a clue that enabled Kruskal, with Clifford S. Gardner, John M. Greene, and Miura (GGKM),[15] to discover a general technique for exact solution of the KdV equation and understanding of its conservation laws.
This discovery gave the modern basis for understanding of the soliton phenomenon: the solitary wave is recreated in the outgoing state because this is the only way to satisfy all of the conservation laws.
In seminal work preceding AKNS, Zakharov and Shabat discovered an inverse scattering method for the nonlinear Schrödinger equation.
In 1986, Kruskal and Zabusky shared the Howard N. Potts Gold Medal from the Franklin Institute "for contributions to mathematical physics and early creative combinations of analysis and computation, but most especially for seminal work in the properties of solitons".
In awarding the 2006 Steele Prize to Gardner, Greene, Kruskal, and Miura, the American Mathematical Society stated that before their work "there was no general theory for the exact solution of any important class of nonlinear differential equations".
The AMS added, "In applications of mathematics, solitons and their descendants (kinks, anti-kinks, instantons, and breathers) have entered and changed such diverse fields as nonlinear optics, plasma physics, and ocean, atmospheric, and planetary sciences.
Much of his subsequent research was driven by a desire to understand this relationship and to develop new direct and simple methods for studying the Painlevé equations.
In Kruskal's opinion, since this property defines the Painlevé equations, one should be able to start with this, without any additional unnecessary structures, to work out all the required information about their solutions.
The first result was an asymptotic study of the Painlevé equations with Nalini Joshi, unusual at the time in that it did not require the use of associated linear problems.
This question was answered negatively in the full generality, for which Conway et al. had hoped, by Costin, Friedman and Ehrlich in 2015.
Kruskal coined the term asymptotology to describe the "art of dealing with applied mathematical systems in limiting cases".