@Preben Alsholm 

Hi thank you one more time for helping
I do not know syntax of maple so well to do what i intend :
PDE2.mw

Looking from the top:
I have chosen u1,u2 to be my control variables. u1 should take values 0.1 or 1.3; u2 should 10 or 30;

doty is an equation that describe how my state y(t) changes. y should be in between 22 and 26 (22<y<26) so there should be some inequality constraints.

parameter Tout should be a function of time - sinusoid with offset preferably with amplitude 10 and offset 10 so it resmbles daily temperature change.

going back to Bellman equation dotW11. At each iteration dotW11 has to be calculated for all combinations(4) off u1 and u2, their min and max values. And max(dotW11(u1,u2)) will decide which combination is the best.

And this will hopefully change according to changes of state y and Tout. If Tout>y or Tout<y.
I have problem with grasping the iterations in maple because there is no visualisation in form of for loop or any similar.

I would appreciate help
Thank you

@Preben Alsholm 

Thank you
That helped a lot