:-op(300,xfy,^).
d(X,X,1):-!.
d(C,X,0):-atomic(C).
d(N*U,X,V):-atomic(N),d(U,X,V).
d(X^N,X,N*X^M):-number(N),M is N-1.
d(U+V,X,H+T):-d(U,X,H),d(V,X,T).


simp(E,E):-atomic(E),!.
simp(E,F):-
    E=..[Op,La,Ra],
    simp(La,Lau), simp(Ra,Rau),
    s(Op,Lau,Rau,F).

s(*,0,_,0).
s(*,_,0,0).
s(+,0,Ra,Ra).
s(+,La,0,La).
s(*,1,Ra,Ra).
s(*,La,1,La).

s(*,N,M*Ra, ):-number(N),number(M),


s(+,La,Ra,La+Ra).
s(*,La,Ra,La*Ra).
s(^,La,Ra,La^Ra).