function T0=g0(n, Dx);

%supponiamo che ci sia temperatura uniforme pari a 20 gradi

T0=20*ones(n+1,1);