power(_, 0, 1). % if instead _ is X, 'warning singleton variable [X]'
power(X, 1, X).
power(X, K, N) :-
  K > 1,
  K1 is K - 1,
  power(X, K1, N1),
  N is X * N1.