Automated control systems can be designed for oxygen therapy to adjust the fraction of inspired oxygen (FiO2) to meet the percentage of total hemoglobin bonded with oxygen (SaO2). The goal of this summer research was to establish a computer model of neonatal oxygen transport systems to be later used as a patient model in the development of automated oxygenation control systems. The model was previously reported in the literature and was reproduced in this work in Simulink, a graphical programming environment in MATLAB. The simulated model consisted of a respiratory sub-model and a cardiovascular submodel linked together using an oxygen dissociation curve, a nonlinear curve relating the oxygen measurement in the arterial blood with SaO2. Simulation results were in agreement with those reported in the literature.