How to write code for optimization with constraints of 2nd order partial difference equation?

im using finite difference method to solve the partial differential eqations. i have to miximize the function

Answer 1

[Matt J]

This is good reading,

https://www.mathworks.com/help/optim/ug/optimizing-a-simulation-or-ordinary-differential-equation.html