function [sol_ana]=teste(boundary_condition,peclet,x,t) % An analytical solution of the diffusioncon

% Parâmetros
U0 = 0;       % Condição de contorno em x=0 (ppm ou mol/m3)
UL = 0;       % Condição inicial (ppm ou mol/m3)
phi = 1;    % Porosidade
vel = 1;      % Velocidade (m/s)
D = 1;      % Coeficiente de difusão (m2/s)
L = 1;        % Comprimento característico (m)
x = 0:(1/100):L; % Malha espacial (m)
t = 0:(1/100):1; % Tempo específico (s)

% Chamada da função
solucao = teste1(U0, UL, phi, vel, D, L, x, 1);

% Plot
plot(x, solucao);
xlabel('x');
ylabel('C');
title(sprintf('Solução analítica no tempo %.2f', 1));