function make_file(filename)
mfile = [filename '.m'];
[fid, message] = fopen(mfile, 'wt');
fprintf(fid, ['function[CV reactions rev_reactions init species start stop increment] = ' filename '()\n']);
fprintf(fid, '%s----Modify the following code to change the model-----\n','%');

fprintf(fid, '\n\n');

fprintf(fid, 'reactions = {''B + P --> I1''      [ 0.1] ;\n');    
fprintf(fid, '             ''I1 + P --> I2''     [0.1]  ;\n'); 
fprintf(fid, '             ''I2 + P --> I3''     [0.1]  ; \n');
fprintf(fid, '             ''I1 --> L''          [0.06667]  ;\n');              
fprintf(fid, '             ''I1 --> K''          [ 0.6]  ;\n');              
fprintf(fid, '             ''I2 --> L''          [0.06667]  ;\n');              
fprintf(fid, '             ''I2 --> K''          [0.15556]  ;\n');              
fprintf(fid, '             ''I3 --> L''          [0.06667 ]  ;\n');              
fprintf(fid, '             ''I3 --> K''          [0.06667]  ;\n');              
fprintf(fid, '             ''K --> 100P''        [0.025];\n'); 
fprintf(fid, '             ''B + FA --> 2B''         [0.05];\n'); 
fprintf(fid, '             ''L + FA --> 2L''          [0.05];\n'); 
fprintf(fid, '             ''FT --> FA''         [0.05]};\n\n');


fprintf(fid, 'rev_reactions = {};\n\n');

fprintf(fid, 'species = {''B'';\n');
fprintf(fid, '             ''P'';\n');
fprintf(fid, '             ''I1'';\n');
fprintf(fid, '             ''I2'';\n');
fprintf(fid, '             ''I3'';\n');          
fprintf(fid, '             ''L'';\n');
fprintf(fid, '             ''K'';\n');
fprintf(fid, '             ''FA'';\n');
fprintf(fid, '             ''FT''};\n\n');


fprintf(fid, '%sInitialValues\n','%');
fprintf(fid, '%sB =\n','%');
fprintf(fid, 'init(:,:,1) = [10000 10000;\n');
fprintf(fid, '               10000 10000];\n\n');

fprintf(fid, '%sP =\n','%');
fprintf(fid, 'init(:,:,2) = [100 0;\n');
fprintf(fid, '               0 0];\n\n');

fprintf(fid, '%sI1 =\n','%');
fprintf(fid, 'init(:,:,3) = [0 0;\n');
fprintf(fid, '               0 0];\n\n');

fprintf(fid, '%sI2 =\n','%');
fprintf(fid, 'init(:,:,4) = [0 0;\n');
fprintf(fid, '               0 0];\n\n');
                
fprintf(fid, '%sI3 =\n','%');
fprintf(fid, 'init(:,:,5) = [0 0;\n');
fprintf(fid, '               0 0];\n\n');
           
fprintf(fid, '%sL =\n','%');
fprintf(fid, 'init(:,:,6) = [0 0;\n');
fprintf(fid, '               0 0];\n\n');
           
fprintf(fid, '%sK =\n','%');
fprintf(fid, 'init(:,:,7) = [0 0;\n');
fprintf(fid, '               0 0];\n\n');
           
fprintf(fid, '%sFA =\n','%');
fprintf(fid, 'init(:,:,8) = [0 0;\n');
fprintf(fid, '               0 0];\n\n');
           
fprintf(fid, '%sFT =\n','%');
fprintf(fid, 'init(:,:,9) = [1e+6 1e+6;\n');
fprintf(fid, '               1e+6 1e+6];\n\n');
           
fprintf(fid, '%sConstant Values\n','%');
fprintf(fid, 'CV = {''B'' ''FA'' ''FT''};\n\n');
           
fprintf(fid, 'start = 0;\n');
fprintf(fid, 'stop = 120;\n');
fprintf(fid, 'increment = 1;\n');

fclose(fid);
