#ftest;

ftest(s_,eq_):= block(f(x)=s,expand(expand(expand(eq))));
ftest(s_,eq_=r_):= block(f(x)=s,expand(expand(expand(eq-r))));