%u(1,0) Initial distribution of heat through the rod from 0<= x <=1
function res = f(x)
    res = sin(pi*x);
    %res = 0 ;
end