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))));