clear all; clf;
%Datos conocidos
t0 = 0;
M0 = 10;
h1 = 0.1;
k = 1.24e-4;
lim = 0.08*M0;

%Metodo de Euler para h = 0.1
    %%Iteracion hasta obtener 0.08*M0
    t1(1) = t0; M1(1) = M0;
    i = 1;
    while M1(i) > lim
        M1(i+1) = M1(i)+h1*(-k*M1(i));
        t1(i+1) = t1(i)+h1;
        i = i+1;
    end

    %%Representacion grafica
    plot(t1,M1,'r')
    grid on
    title('Evolución de la cantidad de Carbono-14 (Método de Euler con paso h = 0.1)')
    xlabel('Tiempo   t    (años)')
    ylabel('Cantidad de Carbono-14   M(t)    (kg)')
    resp = sprintf('La edad de los restos arquelógicos para un paso h = 0.1 es: %d años',t1(end));
    disp(resp)