Табулирование функции в среде Delphi

Условия задачи:
Разработать схему алгоритма и Delphi-проект вычисления таблицы значений функции:
y=ax^4+b*ln⁡(x^2+b^2)
Параметр А принимает N различных значений, начиная с начального значения An и изменяясь с шагом Da, параметр X принимает М различных значений, с шагом Dx, а параметр В численно равен корню уравнения:

√(2.31xcos(πx/2)) - x^3 = 0
Вычисленному на интервале изоляции от c=0 до d=9, с заданной погрешностью Ꜫ=от 10-3 до 10-6.