novaEntradaDados.m

%% Variaveis e funções
ca = 10;
cb = 5;
cinj = 10;
cini = 15;
Lx = 9;
alfa = 0.000008;
tempo = 1000;
divEspaco = 25;
divTempo = 80;
u = 0.001;
deltaX = Lx / divEspaco; % Número de divisões espaciais.
deltaT = tempo / divTempo;  % Número de divisões temporais.
deltaTMAX = (1)/(((2*alfa)/(deltaX^2)) + (u/deltaX));
aux = 0; % Definindo aux = 0 globalmente.
vetorEspaco = zeros(1,divEspaco); % Inicializando vetor da malha espacial.
oldVet = zeros(1,divEspaco); % Inicializando vetor da concentração antiga.
newVet = zeros(1,divEspaco); % Inicializando vetor da concentração nova.


%Preenchendo vetores com as concentrações
for i = 1:divEspaco
    oldVet(i)=cb;
end
for z = 1:divEspaco/2
    oldVet(z)=ca;
end

for i = 1:divEspaco
    newVet(i)=cb;
end
for z = 1:divEspaco/2
    newVet(z)=ca;
end
%% Preenchendo vetores pro plot da condição inicial
vetorCini = zeros(1, tempo);
vetorTempo = zeros(1,tempo);
tamVetorTempo = length(vetorTempo);

for i=1:tempo
    vetorTempo(i) = i;

end
for j=1:tamVetorTempo/2
    vetorCini(j) = cini;
end
for k=(tamVetorTempo/2):tamVetorTempo
    vetorCini(k) = 10;
end




% Gráfico Condição Inicial

% plot(vetorTempo,vetorCini,'r');
% title("Troca da condição inicial");
% xlabel("Lx");
% ylabel("Condição inicial"); 



%%
for i = 2:divEspaco
    vetorEspaco(i) = vetorEspaco(i-1) + deltaX; % Preenchendo a malha começando do i=2 por causa da condição de contorno
end

% Preencher o vetor com as condições iniciais, no while ele vai substituir
% os valores que não sejam o indice 1.
% for i = 1:divEspaco
%     oldVet(i)=cini;
%     newVet(i)=cini;
% end

%oldVet(1) = cini; % Condição inicial no indice 1.
%newVet(1) = cini; % Condição inicial no indice 1.