import("bfct")$
F=generic_bfct([x^2*dx-x*dx+4*x-1],[x],[dx],[-1]);
print(fctr(F));  /* same output with M2 */
end$