:- use_package(clpr).

prod([],[], Result) :-
    Result .=. 0.

prod([X|Xs],[Y|Ys], Result) :-
    Result .=. X * Y + Rest,
    prod(Xs, Ys, Rest).