как задать в фортране  производную в общем виде
или можно какую-нибудь библиотеку подключить

на самом деле, решаю следующее
есть выражение вида x(j+1)=1-a(x(j))^2

на каждом шаге надо вычислять производную dx(j+1)/dx(j),
но скажым шагом полином увеличивается 

подскажите, примерный алгоритм 

Это сообщение отредактировал(а) f999t1 - 22.12.2008, 17:11

Вы бы хоть функцию нормально записали, а то возведение в степень там осталось не у дел перед скобкой, да и не ясно, с чего бы это  

Тогда, может, кто чего и подскажет.

по-моему, библиотеки тут не подключишь, можно попробовать вычислять производную либо аналитически, либо численными методами. Коды вычисления численными методами обычно простые.

так с производными боле мене понятно
вот текст



как задать рекурсию?
x(j+1)=1-a( x(j) )^2

т.е в первом ходе вычислил а из уравнения
fa=1-a*(x^2)
во втором ходе надо это сделать из уравнения
fa2=1-a*(fa^2)=1-a+2(a^3)(x^2)-(a^2)(x^4)
и т.д.

Производные на фортране сложно записать аналитически, не парсер же прикручивать Воспользуйтесь конечными разностями, как по определению производной:



А рекурсию (хотя лучше её избегать) делают с помощью объявлений функции типа 
recurcive function f(x) result A........

Это сообщение отредактировал(а) popovda - 27.12.2008, 15:58