Математическая модель процессов внутри земли

[nkabanoff]

Доброго времени суток! Работаю над проектом по математической физике. Столкнулся с некоторыми трудностями. Во-первых, вычисление производных функций не числено, а аналитически(формульно). Из этого вытекает вопрос с аналитическим нахождением оператора лапласа произвольной функции. Возможно ли написание такого кода на си, может быть есть уже написанные программки? Во-вторых может быть у кого есть программа для решения общих задач Штурма-Лиувилля написанная на си. Понятно, что все это можно сделать в Wolfram Matematica, но мне хотелось бы написать подобную программу на си.