Prolog. Реализация раскрытия геометрических формул?

Question

[FVoyager]

Процедура должна раскрывать формулы синусов и косинусов суммы

sin(A+B) = sin A cos B + cos A sin B,
cos(A+B) = cos A cos B – sin A sin B.
Пример работы
?- simplify(sin(x+cos(3+x)),X).
X = sin(x)*cos(cos(3)*cos(x)-sin(3)*sin(x))+
cos(x)*sin(cos(3)*cos(x)-sin(3)*sin(x))

Основной вопрос заключается в том, как именно передавать выражение, и функции в ней, т.е, что - то вроде
func(expr), if func == sin ? sin(expr) : cos(expr)