%Define function f79 for exercise 7.9
function f=f79(u)
   
   
   alpha=0.15;
   beta=0.9;
   mu0=0.1; 
   g=9.81;
   f=[u(2) -mu0*(1+alpha*u(1)^beta)*g];
