model HarmonicS2P1Dirichlet.mod;
data  DataProb3_S2_Dirichletp1.dat;
option ipopt_options 'wantsol=2';


write gHarmonicS2P1Dirichlet;
