Every applied mathematician has used separation of variables. For a given boundary value problem (BVP) in two dimensions, the starting point of this powerful method is the separation of the given PDE into two ODEs. If the spectral analysis of either of these ODEs yields an appropriate transform pair, i.e., a transform consistent with the given boundary conditions, then the given BVP can be reduced to a BVP for an ODE. For simple BVPs it is straightforward to choose an appropriate transform and hence the spectral analysis can be avoided. In spite of its enormous applicability, this method has certain limitations. In particular, it requires the given domain, PDE, and boundary conditions to be separable, and also may not be applicable if the BVP is non-self-adjoint. Furthermore, it expresses the solution as either an integral or a series, neither of which are uniformly convergent on the boundary of the domain (for nonvanishing boundary conditions), which renders such expressions unsuitable for numerical computations. This paper describes a recently introduced transform method that can be applied to certain nonseparable and non-selfadjoint problems. Furthermore, this method expresses the solution as an integral in the complex plane that is uniformly convergent on the boundary of the domain. The starting point of the method is to write the PDE as a one-parameter family of equations formulated in a divergence form, and this allows one to consider the variables together. In this sense, the method is based on the "synthesis" as opposed to the "separation" of variables. The new method has already been applied to a plethora of BVPs and furthermore has led to the development of certain novel numerical techniques. However, a large number of related analytical and numerical questions remain open. This paper illustrates the method by applying it to two particular non-self-adjoint BVPs: one for the linearized KdV equation formulated on the half-line, and the other for the Helmholtz equation in the exterior of the disc (the latter is non-self-adjoint due to the radiation condition). The former problem played a crucial role in the development of the new method, whereas the latter problem was instrumental in the full development of the classical transform method. Although the new method can now be presented using only classical techniques, it actually originated in the theory of certain nonlinear PDEs called integrable, whose crucial feature is the existence of a Lax pair formulation. It is shown here that Lax pairs provide the generalization of the divergence formulation from a separable linear to an integrable nonlinear PDE.