The problem of simulating complex systems, such as production lines, industrial plants, food processing, etc., today represents an opportunity that brings with it the great advantage of limiting design costs. However, nowadays the designer, after defining and implementing the mathematical models of the studied process, may need to rebuild the whole simulation framework because he needs to modify the model of even just one subsystem. It is for this reason that in this paper, a new framework for the use of Individual Subsystem Models (ISM) for the modelling and simulation of interconnected systems has been studied and implemented. Furthermore, the study of the state of the art has revealed the lack of efficient and sufficiently general numerical algorithms, but, at the same time, it is simple to use to solve the algebraic-differential equations deriving from the ISM simulation. The proposed new approach follows the paradigm of co-simulation methods, including graph theory methods, to solve the general ISM simply and efficiently. In this approach, each subsystem is required to have its own representation independently of the other subsystems. In this way, it is always possible to replace any subsystem whenever an updated representation becomes available, making maintenance and evolution of the entire ISM flexible. Our framework calls each subsystem separately in an optimal (suboptimal) order based on the structure of the graph. Each calculated output is transferred to the input of the next subsystem in the chosen. The general procedure has been validated in the context of Linear and Time-Invariant ISMs: in these hypotheses, the stability conditions have been calculated and numerical tests have been performed which show the effectiveness of the proposed approach.