The object of this study is to develop a mathematical model of the hemodialysis system including the mechanism of solute kinetics, water exchange and also cardiovascular dynamics. The cardiovascular system model used in this study simulates the short-term transient and steady-state hemodynamic responses such as hypotension and disequilibrium syndrome (which are main complications to hemodialysis patients) during hemodialysis. It consists of a 12 lumped-parameter representation of the cardiovascular circulation connected to set-point models of the arterial baroreflexes, a kinetic model (hemodialysis system model) with 3 compartmental body fluids and 2 compartmental solutes. We formulate mathematically this model in terms of an electric analog model. All resistors and most capacitors are assumed to be linear. The control mechanisms are mediated by the information detected from arterial pressoreceptors, and they work on systemic arterial resistance, heart rate, and systemic venous unstressed volume. The hemodialysis model includes the dynamics of urea, creatinine, sodium and potassium in the intracellular and extracellular pools as well as fluid balance equations for the intracellular, interstitial, and plasma volumes. Model parameters are largely based on literature values. We have presented the results on the simulations performed by changing some model parameters with respect to their basal values. In each case, the percentage changes of each compartmental pressure, heart rate (HR), total systemic resistance (TSR), ventricular compliance, zero pressure filling volume and solute concentration profiles are represented during hemodialysis.