%
% MATLAB5 : script de test 
%

function matlab5()

    close all
    Tstart = 0;
    Tend   = 5e5;
    Ustart = [1e8 20];
    n = 5e2;
    [T,U] = forest(Tstart,Tend,Ustart,n);
    
    figure
    subplot(2,1,1)
    plot(T/365,U(:,1),'b')
    title('Wood in the forest ')
    xlabel('Time [years]')
    ylabel('Global wood mass [kg]')
    subplot(2,1,2)
    plot(T/365,U(:,2),'r')
    title('Population')
    xlabel('Time [years]')
    ylabel('Population')

    figure
    plot(U(:,1),U(:,2));
    xlabel('Global wood mass [kg]')
    ylabel('Population')

  
    finalValue = U(end,1);
end