Although molecular electronics is in its infancy, there are many significant advances in understanding and implementing computational functions based on molecular electronic devices. In order to evaluate a molecular electronic circuit, a model is needed. In existing literature, only a few models have been introduced for molecular electronics, all of which do not support the hysteresis phenomenon. To the best of the authors' knowledge, for the first time, this paper introduces a behavioral model in SPICE which is able to imitate hysteresis behavior in a desired molecular circuit. The proposed model is then modified in order to guarantee continuous behavior in its I/V curves. Some molecular electronic I/V curves from the literature are selected and by utilizing the proposed model, its correct operation are investigated.