Some problems with c, q, and Ca variables. Implemented 4K explicity Runge-Kutte ODE solver