import("bfct")$
II=ann( t*(x-t) );
I=map(subst,II,s,1/2);
/* No Fourier command */
J=[- t*dt - 2*x*dt + x*dx, x*dt - x*dx - 3/2];
generic_bfct(J,[t,x],[dt,dx],[1,0]);
/* s(s+2) */
dp_weyl_gr_main(J | v=[t,x,dt,dx], order=[[t,-1,dt,1]]);
end$