The work is aimed to establish mathematical model for functional dependencies of anthropomethric characteristics on nutritional assessment for adolescent pupils in Croatian boarding schools. Investigated are 1117 female and male pupils from 14 to 19 years of age. Experimentally are determined the following parameters: body mass, body weight, body mass index, obesity index, circumference (waist, hips, chest, arm), skin fold thicknesses ( subscapular, triceps). Proposed are multivariable linear models for relationships between the parameters of nutritional assessment and intake of principal food components. They are related to intake of energy, protein, fat and carbohydrates. The models are determined for each female and male generation, and evolution of the model parameters through generations is determined. In view of recommended diet allowances (RDA) and econometric considerations, proposed are linear programming models for optimisation of diet planning in the Croatian boarding schools.