  We demonstrate the use of the Unified Transform Method or Method of Fokas for
boundary value problems for systems of constant-coefficient linear partial
differential equations. We discuss how the apparent branch singularities
typically appearing in the global relation are removable, allowing the method
to proceed, in essence, as for scalar problems. We illustrate the use of the
method with boundary value problems for the Klein-Gordon equation and the
linearized Fitzhugh-Nagumo system. The case of wave equations is treated
separately in an appendix.
