Matlab backup 4

% Parâmetros
u0 = 1;           % Condição de contorno em x=0
phi = 1;        % Porosidade
V = 1;            % Velocidade
D = 0;            % Coeficiente de difusão
alpha = -1;
L = 1;           % Comprimento do domínio (0,L)
T = 1;            % Tempo final

dx = 0.001;
x = 0:dx:L;  % Malha espacial

dt = 0.01;
t = 0:dt:T;  % Malha temporal

% Solução Analítica
f = @(x, t) (x <= (V/phi) * t) .* exp((alpha/phi) * x);

tempos = [0.2, 0.5, 0.7];
% Plotando a solução analítica
hold on;
for j = 1:length(t)
    if t(j) == tempos(1) || t(j) == tempos(2) || t(j) == tempos(3)
        plot(x, f(x, t(j)), '--', 'Color', [0 0 1]); % Plota em azul
    end
end

hold off;
xlabel('x');
ylabel('u');
title('Comparação da Solução Numérica com Analítica');
grid on;