g_0t(y_, t), size(u(:, 1)));
%           u(1, :) = reshape(g_0t(x_, t), size(u(1, :)));
%           u(end, :) = zeros(size(u(end, :)));
%           u(:, end) = zeros(size(u(:, end)));
%           u_t = u + A * u + u * B;
%           if (any(id == slices_id))
%               plotSlice_4_2_3(u_t, X_, Y_, dt, id);
%           end
%           u = u_t;
%       end
%   end
%
%%%
%% Code: Value of Slices
%%%
%   
%   function plotSlice_4_2_3(u, X_, Y_, dt, id)
%       figure; hold on;
%       title(sprintf('t = %.