The paper presents a proposal for an in-house mathematical model of a liquid flat-plate solar collector. The proposed model is a one-dimensional distributed parameter model enabling simulations of the collector operation under arbitrarily variable boundary conditions. The model is based on the solution of energy balance equations for all components of the collector. The formulated differential equations are solved iteratively using an implicit difference scheme. In order to obtain a stable numerical solution, it is necessary to use appropriate steps of the time and spatial division. These were found by comparing the results obtained from the model with the results of the analytical solution available in the literature for the transient state, which constitutes the novelty of the present study. The accuracy of the results obtained from the model was verified experimentally by comparing the measured and the calculated history of the fluid temperature at the outlet of the collector. The calculation of the collector time constant is proposed in the paper as an example of the model practical application. The results of the time constant calculation were compared with the values obtained experimentally on the test stand. This is another novelty of the presented research. The analysed collector instantaneous efficiency was then calculated for selected outdoor conditions. The presented mathematical model can also be used to verify the correctness of the collector operation. Comparing, on an ongoing basis, the measured and the calculated values of the fluid temperature at the collector outlet, conclusions can be drawn about the process of solar glass fouling or glycol gelling. The model's simplicity and the low computational demands enable such comparisons in an online mode.