%u(1,t) variaton in heat at endpoint, x=1 function res = g1(t) res = 0; end