# Simple A+B <-> C equilibrium # Output is text to stdout with numbers of A, B, and C molecules. # This runs with either 1 or 3 molecule lists; comment and uncomment code below as required. # I found runtime was 7.37 s with 3 lists and 10.47 s with 1 list. graphics none graphic_iter 10 dim 3 species A B C max_mol 10000 boxsize 4 accuracy 10 # For 1 list, use following code: /* molecule_lists onlylist mol_list all onlylist */ # For 3 lists, use following code: molecule_lists Alist Blist Clist mol_list A Alist mol_list B Blist mol_list C Clist difc A 1 difc B 1 difc C 1 color A 1 0 0 color B 0 1 0 color C 0 0 1 time_start 0 time_stop 100 time_step 0.01 boundaries 0 0 100 p boundaries 1 0 100 p boundaries 2 0 100 p cmd i 0 100 10 molcount mol 1000 A u u u mol 1000 B u u u start_reaction 1 max_rxn 1 reactant C revrxn rate revrxn 0.1 product revrxn A + B product_param revrxn x 0.2 end_reaction start_reaction 2 max_rxn 1 reactant A + B fwdrxn rate fwdrxn 100 product fwdrxn C end_reaction end_file