I am writing a C program to numerically integrate a differential equation. This differential equation corresponds to dynamics of a system. I am using rkf45 to integrate it. But the problem with rkf45 is that it requires Jacobian of the equations. I want to get rid of jacobian. ode45 in MATLAB does the same thing w/o using Jacobian. Is there any way to get rid of jacobian in C ???