SoloCodigo

hola: estoy haciendo un programa en matlab, el codigo es el siguente

function[C2iterado C3iterado C]=calcc1(D,rop,C1,C4,L,k,K), tic
format long
clc
clf
Citer=1067911;
C2=C1-1;
J1=rop.*D.*(C1-C2)./L;
C3=0;
C=[C1 C2 J1 C3];
deltat=1;
Vp=7.3E-7;
Vas=1.23E-4;
A=5E-3;
for i=1:10000
    C(i+1,1)=C(1,1);
    C(i+1,2)=C(i,2)-1;
    C(i+1,3)=rop.*D.*(C(i+1,1)-C(i+1,2))./L;
    C(i+1,4)=C(i+1,3)./k;
    if C(i,2)<=C(i,4).*K;
        break
    end
end
C
C3iterado=[C(size(C,1),4)];
C2iterado=[C(size(C,1),2)];
x=[1:size(C(:,2))];
subplot(1,3,1),plot(x,C(:,2))
xlabel('i')
ylabel('C2 [mg/m^3]')
subplot(1,3,2),plot(x,C(:,3))
xlabel('i')
ylabel('J1 [mg/(m^2*s)]')
subplot(1,3,3),plot(x,C(:,4))
xlabel('i')
ylabel('C3 [mg/m^3]')
CiEE=[C2iterado C3iterado];
C1paranuevaiteracion=((1-rop.*D.*A.*deltat./(L.*Vp)).*C(1,1)+(rop.*D.*A.*deltat./(L.*Vp)).*C2iterado);
C4paranuevaiteracion=((1-k.*A.*deltat./Vas).*C(1,4)+(k.*A.*deltat./Vas).*C3iterado);

    if C1paranuevaiteracion>=Citer.*0.9
    calcc1(D,rop,C1paranuevaiteracion,C4paranuevaiteracion,L,k,K)
   
    end
toc
 necesito que en el termino definido como CiEE se guarde el resultado de la primera solucion (como ya esta escrito) y a continuacion al hacer el if empiece a guardar los resultados para cada nueva solucion del problema, para obtener una matriz CiEE con la soluciones de todos los calculos que hizo, alguien sabe como hacerlo???

si no se entiende lo q necesito porfavor avisenme

saludos
y gracias de ante mano