The problem of finding a rational function of a real variable interpolating at given points (with suitable intersection numbers) while being everywhere positive is considered. The positivity requirement is essential, e.g., in modeling the square magnitude of a frequency response. This goal is achieved by resorting to a Hermite-type approach which allows one to obtain positive interpolating functions of minimum complexity by means of an efficient numerical procedure
