Pick functions are the analytic maps of the upper half plane into itself. Their boundary behavior was studied classically by many household names such as Caratheodory, Julia and Nevanlinna. For example, Nevanlinna showed that asymptotic expansions of a certain type of Pick function at infinity parametrize the solutions to the Hamburger moment problem. We propose a program for studying interpolation and approximation multivariable Pick functions on the boundary using techniques from functional analysis which have roots in systems engineering.