simple factorization

% factor
:- use_module(library(clpfd)).

/*
prime(2).
prime(3).
prime(5).
prime(7).
prime(11).
*/

composite(N,X,Y) :-
  N #> 0, X #> 0, Y #> 0,
  X  #\= 1, Y  #\= 1,
  X  #\=  N , Y  #\= N,
  fac(N,X,Y).

prime(N) :-
  N #> 0, N #\=1 ,
  \+ composite(N,X,Y).

fac(N,X,Y) :- X * Y #= N.

simple_factor(N,X,Y) :-
  prime(X), prime(Y),
  X * Y #= N.