This work presents a method for the mathematical modelling of dynamics of complex environmental systems and ecological processes called the "response function method". Using time series observations and prior knowledge or information, thisapproach allows the researcher to take all essential features of ecosystem processes, such as complexity, multi-dimensioness, uncertainty and irreducability, into account. The authors demonstrate the potential of the method using several models, and a diskette containing these is included.