Resolução de sistemas de EDOs por RK4
vetor solução da variável dependente
valor da variável independente
raiz procurada
instante inicial
tamanho do passo
condições iniciais de y
[ vet_sol, v_indep ] = cpc_rk4( fun, 0, 10, 1, [0,0] );