= u(t_slice_id, :);
%       u_slice_num = u_slice_num(:);
%       figure; hold on;
%       plot(x_, u_slice_num, 'r-');
%       legend('numerical');
%       xlabel('x'); ylabel('u');
%   end
%
%%%
% $t = 0.06$
plotSlice_2_6_3( 0.06, dt, u, x_ );
%%%
% $t = 0.1$
plotSlice_2_6_3( 0.1, dt, u, x_ );
%%%
% $t = 0.9$
plotSlice_2_6_3( 0.9, dt, u, x_ );
u_slice_1 = u(round(0.9 / dt) + 1, :);
u_slice_1 = u_slice_1(:);
dt = 0.1; T = round(1 / dt); N = t * T;
r = nu * dt / dx / dx;
s = a * dt / dx / 2;
t_ = 0 : dt : t