Inicio do calculo para rungekutta

function rk4()