решение СОДУ методом Рунге-Кутты

[kraftig]

Всем привет. Есть система, которую надо решить с определенным количеством шагов и их количеством.
Вопросы такие:
1) Какой пакет поможет это сделать максимально простым способом? Может примерчик какой дадите?
2) Есть ли у кого библиотека или свой алгоритм на C/C++, который можно использовать для решения системы.
В обоих случаях метод должен быть 4-5 порядка.

\begin{cases} \dot z_1=z_2\\ \dot z_2=z_3\\ \dot z_3=-\frac{65z_3+14z_2+z_1-5}{100} \end{cases}

[cool.phenon]

Есть такой пакет. Например, MATLAB.