This paper describes the algorithm of a program used to simulate three dimensional models of molecules. In addition to open ended molecules the program also enables simulation of structures with constraints in the form of cyclic regions or fixed location of particular atoms. Several molecules can be handled in a single run and each molecule can have any number of contraints. Further, any number of conformations can be obtained for each constrained region. The program can be used for research in several areas of molecular biology, e.g., structure determination, conformational analysis and topographic comparisons.