function [ f ] = example_7p4_func( t, u )
% Function for example 7.8 of the manuscript.

global N;
f = zeros(N-1,1);
h = 1/N;
f(1) = N^2*(u(2) - 2*u(1));
for i=2:N-2
    f(i) = N^2*(u(i+1)-2*u(i)+u(i-1));
end
f(N-1) = N^2*(t - 2*u(N-1) + u(N-2));
    
%f(1) = 16*u(2) - 32*u(1);
%f(2) = 16*u(3) - 32*u(2) + 16*u(1);
%f(3) = 16*t - 32*u(3) + 16*u(2);

end

