%% Valeur maximale de z et intervalle de valeurs de z
zmax= 10; zint= [0 zmax];
% Parametre de tir
alpha= 0.2;

%% Appel a ode45 pour resolution de l'equation differentielle
% 1er argument = fonction qui definit la derivee,
% 2d argument = intervalle des valeurs de z,
% 3eme argument = vecteur des conditions initiales
[z,F]= ode45(@G , zint , [ 0 ; 0 ; alpha ]);

% Recuperation des valeurs tabulees de la fonction et de sa derivee
f= F(:,1); fprime= F(:,2);

% Graphe de f'(z)
plot(z,fprime); xlabel("z"); ylabel("fprime");

%% Fonction G qui definit l'equation differentielle
function res= G(z, F)
  res= [ F(2) ; F(3) ; ... ];
end
