Go to the directory problem1. Create a file called ConvDiffReacSolver.m and write your program in this file. You also have to edit the files Reaction1.m and Reaction2.m, where you specify the variational formulation.