function GraphingSubroutine(x,t,u,Nt,jump)

% Solution for different times (animated)
for n=1:jump:Nt
    plot(x(:,n),u(:,n));
    hold on;
    pause(0.1);
%    pause;
end

figure;
% Solution in 3-D
for n=1:jump:Nt
    plot3(x(:,n),t(:,n),u(:,n));
    hold on;
    pause(0.2);
%    pause;
end

% Solution over xt-plane
figure;
surf(x,t,u);