eyeprolog 1.6.12 → 1.6.13

package/examples/proof/peano-calculus.pl

@@ -0,0 +1,294 @@
+peano_answer(two_plus_three, s(s(s(s(s(z)))))).
+peano_answer(two_times_three, s(s(s(s(s(s(z))))))).
+peano_answer(factorial_four, s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))).
+
+clause(1, padd(var('A'), z, var('A')), true).
+clause(2, padd(var('A'), s(var('B')), s(var('C'))), padd(var('A'), var('B'), var('C'))).
+clause(3, pmul(anonymous(1), z, z), true).
+clause(4,
+       pmul(var('A'), s(var('B')), var('C')),
+       (pmul(var('A'), var('B'), var('D')), padd(var('A'), var('D'), var('C')))).
+clause(5, pfact(var('N'), var('Value')), pfac(var('N'), s(z), var('Value'))).
+clause(6, pfac(z, var('Acc'), var('Acc')), true).
+clause(7,
+       pfac(s(var('N')), var('Acc'), var('Value')),
+       (pmul(var('Acc'), s(var('N')), var('Next')), pfac(var('N'), var('Next'), var('Value')))).
+clause(8, peano_answer(two_plus_three, var('N')), padd(s(s(z)), s(s(s(z))), var('N'))).
+clause(9, peano_answer(two_times_three, var('N')), pmul(s(s(z)), s(s(s(z))), var('N'))).
+clause(10, peano_answer(factorial_four, var('N')), pfact(s(s(s(s(z)))), var('N'))).
+
+step(peano_answer(two_plus_three, s(s(s(s(s(z)))))),
+     rule(8),
+     ['N' = s(s(s(s(s(z)))))],
+     [padd(s(s(z)), s(s(s(z))), s(s(s(s(s(z))))))]).
+step(padd(s(s(z)), s(s(s(z))), s(s(s(s(s(z)))))),
+     rule(2),
+     ['A' = s(s(z)), 'B' = s(s(z)), 'C' = s(s(s(s(z))))],
+     [padd(s(s(z)), s(s(z)), s(s(s(s(z)))))]).
+step(padd(s(s(z)), s(s(z)), s(s(s(s(z))))),
+     rule(2),
+     ['A' = s(s(z)), 'B' = s(z), 'C' = s(s(s(z)))],
+     [padd(s(s(z)), s(z), s(s(s(z))))]).
+step(padd(s(s(z)), s(z), s(s(s(z)))),
+     rule(2),
+     ['A' = s(s(z)), 'B' = z, 'C' = s(s(z))],
+     [padd(s(s(z)), z, s(s(z)))]).
+step(padd(s(s(z)), z, s(s(z))), fact(1), ['A' = s(s(z))], []).
+step(peano_answer(two_times_three, s(s(s(s(s(s(z))))))),
+     rule(9),
+     ['N' = s(s(s(s(s(s(z))))))],
+     [pmul(s(s(z)), s(s(s(z))), s(s(s(s(s(s(z)))))))]).
+step(pmul(s(s(z)), s(s(s(z))), s(s(s(s(s(s(z))))))),
+     rule(4),
+     ['A' = s(s(z)), 'B' = s(s(z)), 'C' = s(s(s(s(s(s(z)))))), 'D' = s(s(s(s(z))))],
+     [pmul(s(s(z)), s(s(z)), s(s(s(s(z))))), padd(s(s(z)), s(s(s(s(z)))), s(s(s(s(s(s(z)))))))]).
+step(pmul(s(s(z)), s(s(z)), s(s(s(s(z))))),
+     rule(4),
+     ['A' = s(s(z)), 'B' = s(z), 'C' = s(s(s(s(z)))), 'D' = s(s(z))],
+     [pmul(s(s(z)), s(z), s(s(z))), padd(s(s(z)), s(s(z)), s(s(s(s(z)))))]).
+step(pmul(s(s(z)), s(z), s(s(z))),
+     rule(4),
+     ['A' = s(s(z)), 'B' = z, 'C' = s(s(z)), 'D' = z],
+     [pmul(s(s(z)), z, z), padd(s(s(z)), z, s(s(z)))]).
+step(pmul(s(s(z)), z, z), fact(3), [], []).
+step(padd(s(s(z)), s(s(s(s(z)))), s(s(s(s(s(s(z))))))),
+     rule(2),
+     ['A' = s(s(z)), 'B' = s(s(s(z))), 'C' = s(s(s(s(s(z)))))],
+     [padd(s(s(z)), s(s(s(z))), s(s(s(s(s(z))))))]).
+step(peano_answer(factorial_four, s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(10),
+     ['N' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))],
+     [pfact(s(s(s(s(z)))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pfact(s(s(s(s(z)))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(5),
+     ['N' = s(s(s(s(z)))),
+      'Value' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))],
+     [pfac(s(s(s(s(z)))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pfac(s(s(s(s(z)))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(7),
+     ['N' = s(s(s(z))),
+      'Acc' = s(z),
+      'Value' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'Next' = s(s(s(s(z))))],
+     [pmul(s(z), s(s(s(s(z)))), s(s(s(s(z))))),
+      pfac(s(s(s(z))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pmul(s(z), s(s(s(s(z)))), s(s(s(s(z))))),
+     rule(4),
+     ['A' = s(z), 'B' = s(s(s(z))), 'C' = s(s(s(s(z)))), 'D' = s(s(s(z)))],
+     [pmul(s(z), s(s(s(z))), s(s(s(z)))), padd(s(z), s(s(s(z))), s(s(s(s(z)))))]).
+step(pmul(s(z), s(s(s(z))), s(s(s(z)))),
+     rule(4),
+     ['A' = s(z), 'B' = s(s(z)), 'C' = s(s(s(z))), 'D' = s(s(z))],
+     [pmul(s(z), s(s(z)), s(s(z))), padd(s(z), s(s(z)), s(s(s(z))))]).
+step(pmul(s(z), s(s(z)), s(s(z))),
+     rule(4),
+     ['A' = s(z), 'B' = s(z), 'C' = s(s(z)), 'D' = s(z)],
+     [pmul(s(z), s(z), s(z)), padd(s(z), s(z), s(s(z)))]).
+step(pmul(s(z), s(z), s(z)),
+     rule(4),
+     ['A' = s(z), 'B' = z, 'C' = s(z), 'D' = z],
+     [pmul(s(z), z, z), padd(s(z), z, s(z))]).
+step(pmul(s(z), z, z), fact(3), [], []).
+step(padd(s(z), z, s(z)), fact(1), ['A' = s(z)], []).
+step(padd(s(z), s(z), s(s(z))),
+     rule(2),
+     ['A' = s(z), 'B' = z, 'C' = s(z)],
+     [padd(s(z), z, s(z))]).
+step(padd(s(z), s(s(z)), s(s(s(z)))),
+     rule(2),
+     ['A' = s(z), 'B' = s(z), 'C' = s(s(z))],
+     [padd(s(z), s(z), s(s(z)))]).
+step(padd(s(z), s(s(s(z))), s(s(s(s(z))))),
+     rule(2),
+     ['A' = s(z), 'B' = s(s(z)), 'C' = s(s(s(z)))],
+     [padd(s(z), s(s(z)), s(s(s(z))))]).
+step(pfac(s(s(s(z))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(7),
+     ['N' = s(s(z)),
+      'Acc' = s(s(s(s(z)))),
+      'Value' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'Next' = s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))],
+     [pmul(s(s(s(s(z)))), s(s(s(z))), s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))),
+      pfac(s(s(z)), s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pmul(s(s(s(s(z)))), s(s(s(z))), s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))),
+     rule(4),
+     ['A' = s(s(s(s(z)))),
+      'B' = s(s(z)),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'D' = s(s(s(s(s(s(s(s(z))))))))],
+     [pmul(s(s(s(s(z)))), s(s(z)), s(s(s(s(s(s(s(s(z))))))))),
+      padd(s(s(s(s(z)))), s(s(s(s(s(s(s(s(z)))))))), s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))]).
+step(pmul(s(s(s(s(z)))), s(s(z)), s(s(s(s(s(s(s(s(z))))))))),
+     rule(4),
+     ['A' = s(s(s(s(z)))), 'B' = s(z), 'C' = s(s(s(s(s(s(s(s(z)))))))), 'D' = s(s(s(s(z))))],
+     [pmul(s(s(s(s(z)))), s(z), s(s(s(s(z))))),
+      padd(s(s(s(s(z)))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(z)))))))))]).
+step(pmul(s(s(s(s(z)))), s(z), s(s(s(s(z))))),
+     rule(4),
+     ['A' = s(s(s(s(z)))), 'B' = z, 'C' = s(s(s(s(z)))), 'D' = z],
+     [pmul(s(s(s(s(z)))), z, z), padd(s(s(s(s(z)))), z, s(s(s(s(z)))))]).
+step(pmul(s(s(s(s(z)))), z, z), fact(3), [], []).
+step(padd(s(s(s(s(z)))), z, s(s(s(s(z))))), fact(1), ['A' = s(s(s(s(z))))], []).
+step(padd(s(s(s(s(z)))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(z))))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = s(s(s(z))), 'C' = s(s(s(s(s(s(s(z)))))))],
+     [padd(s(s(s(s(z)))), s(s(s(z))), s(s(s(s(s(s(s(z))))))))]).
+step(padd(s(s(s(s(z)))), s(s(s(z))), s(s(s(s(s(s(s(z)))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = s(s(z)), 'C' = s(s(s(s(s(s(z))))))],
+     [padd(s(s(s(s(z)))), s(s(z)), s(s(s(s(s(s(z)))))))]).
+step(padd(s(s(s(s(z)))), s(s(z)), s(s(s(s(s(s(z))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = s(z), 'C' = s(s(s(s(s(z)))))],
+     [padd(s(s(s(s(z)))), s(z), s(s(s(s(s(z))))))]).
+step(padd(s(s(s(s(z)))), s(z), s(s(s(s(s(z)))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = z, 'C' = s(s(s(s(z))))],
+     [padd(s(s(s(s(z)))), z, s(s(s(s(z)))))]).
+step(padd(s(s(s(s(z)))), s(s(s(s(s(s(s(s(z)))))))), s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))),
+      'B' = s(s(s(s(s(s(s(z))))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(z)))))))))))],
+     [padd(s(s(s(s(z)))), s(s(s(s(s(s(s(z))))))), s(s(s(s(s(s(s(s(s(s(s(z))))))))))))]).
+step(padd(s(s(s(s(z)))), s(s(s(s(s(s(s(z))))))), s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = s(s(s(s(s(s(z)))))), 'C' = s(s(s(s(s(s(s(s(s(s(z))))))))))],
+     [padd(s(s(s(s(z)))), s(s(s(s(s(s(z)))))), s(s(s(s(s(s(s(s(s(s(z)))))))))))]).
+step(padd(s(s(s(s(z)))), s(s(s(s(s(s(z)))))), s(s(s(s(s(s(s(s(s(s(z))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = s(s(s(s(s(z))))), 'C' = s(s(s(s(s(s(s(s(s(z)))))))))],
+     [padd(s(s(s(s(z)))), s(s(s(s(s(z))))), s(s(s(s(s(s(s(s(s(z))))))))))]).
+step(padd(s(s(s(s(z)))), s(s(s(s(s(z))))), s(s(s(s(s(s(s(s(s(z)))))))))),
+     rule(2),
+     ['A' = s(s(s(s(z)))), 'B' = s(s(s(s(z)))), 'C' = s(s(s(s(s(s(s(s(z))))))))],
+     [padd(s(s(s(s(z)))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(z)))))))))]).
+step(pfac(s(s(z)), s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(7),
+     ['N' = s(z),
+      'Acc' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'Value' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'Next' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))],
+     [pmul(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(z)), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+      pfac(s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pmul(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(z)), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(4),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(z),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'D' = s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))],
+     [pmul(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))),
+      padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pmul(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))),
+     rule(4),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = z,
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'D' = z],
+     [pmul(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), z, z),
+      padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), z, s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))]).
+step(pmul(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), z, z), fact(3), [], []).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), z, s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))),
+     fact(1),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))],
+     []).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(s(s(s(s(s(s(z))))))))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(z))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(s(z))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(s(s(s(s(s(z)))))))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(z)))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(s(z)))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(s(s(s(s(z))))))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(z))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(s(z))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(s(s(s(z)))))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(z)))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(s(z)))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(s(s(z))))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(z))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(s(z))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(s(z)))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(z)))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(s(z)))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(s(z))))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(z))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(s(z))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(s(z)))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(s(z)))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(s(z))),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(z))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(s(z))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(s(z)),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(z)), s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(s(z)), s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = s(z),
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))]).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))),
+     rule(2),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))),
+      'B' = z,
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))],
+     [padd(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))), z, s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))]).
+step(pfac(s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(7),
+     ['N' = z,
+      'Acc' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'Value' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'Next' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))],
+     [pmul(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+      pfac(z, s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pmul(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), s(z), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     rule(4),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'B' = z,
+      'C' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))),
+      'D' = z],
+     [pmul(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), z, z),
+      padd(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), z, s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))))]).
+step(pmul(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), z, z),
+     fact(3),
+     [],
+     []).
+step(padd(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), z, s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     fact(1),
+     ['A' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))],
+     []).
+step(pfac(z, s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z)))))))))))))))))))))))), s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))),
+     fact(6),
+     ['Acc' = s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(s(z))))))))))))))))))))))))],
+     []).

package/examples/proof/pell-equation.pl

@@ -0,0 +1,295 @@
+pell_answer(solution_5, solution(3363, 2378)).
+pell_answer(solution_8, solution(665857, 470832)).
+pell_answer(check_8, true).
+pell_answer(y_sum_1_to_8, 568344).
+
+clause(2, pell(0, 1, 0), true).
+clause(3,
+       pell(var('N'), var('X'), var('Y')),
+       (var('N') > 0,
+        var('N1') is var('N') - 1,
+        pell(var('N1'), var('X0'), var('Y0')),
+        var('Ax') is 3 * var('X0'),
+        var('By') is 4 * var('Y0'),
+        var('X') is var('Ax') + var('By'),
+        var('Cx') is 2 * var('X0'),
+        var('Dy') is 3 * var('Y0'),
+        var('Y') is var('Cx') + var('Dy'))).
+clause(4,
+       pell_holds(var('N'), true),
+       (pell(var('N'), var('X'), var('Y')),
+        var('X2') is var('X') * var('X'),
+        var('Y2') is var('Y') * var('Y'),
+        var('Twicey2') is 2 * var('Y2'),
+        1 is var('X2') - var('Twicey2'))).
+clause(5, pell_answer(solution_5, solution(var('X'), var('Y'))), pell(5, var('X'), var('Y'))).
+clause(6, pell_answer(solution_8, solution(var('X'), var('Y'))), pell(8, var('X'), var('Y'))).
+clause(7, pell_answer(check_8, true), pell_holds(8, true)).
+clause(8,
+       pell_answer(y_sum_1_to_8, var('Sum')),
+       sumall(var('Y'), (between(1, 8, var('N')), pell(var('N'), anonymous(1), var('Y'))), var('Sum'))).
+
+step(pell_answer(solution_5, solution(3363, 2378)),
+     rule(5),
+     ['X' = 3363, 'Y' = 2378],
+     [pell(5, 3363, 2378)]).
+step(pell(5, 3363, 2378),
+     rule(3),
+     ['N' = 5,
+      'X' = 3363,
+      'Y' = 2378,
+      'N1' = 4,
+      'X0' = 577,
+      'Y0' = 408,
+      'Ax' = 1731,
+      'By' = 1632,
+      'Cx' = 1154,
+      'Dy' = 1224],
+     [5 > 0,
+      4 is 5 - 1,
+      pell(4, 577, 408),
+      1731 is 3 * 577,
+      1632 is 4 * 408,
+      3363 is 1731 + 1632,
+      1154 is 2 * 577,
+      1224 is 3 * 408,
+      2378 is 1154 + 1224]).
+step(5 > 0, builtin, [], []).
+step(4 is 5 - 1, builtin, [], []).
+step(pell(4, 577, 408),
+     rule(3),
+     ['N' = 4,
+      'X' = 577,
+      'Y' = 408,
+      'N1' = 3,
+      'X0' = 99,
+      'Y0' = 70,
+      'Ax' = 297,
+      'By' = 280,
+      'Cx' = 198,
+      'Dy' = 210],
+     [4 > 0,
+      3 is 4 - 1,
+      pell(3, 99, 70),
+      297 is 3 * 99,
+      280 is 4 * 70,
+      577 is 297 + 280,
+      198 is 2 * 99,
+      210 is 3 * 70,
+      408 is 198 + 210]).
+step(4 > 0, builtin, [], []).
+step(3 is 4 - 1, builtin, [], []).
+step(pell(3, 99, 70),
+     rule(3),
+     ['N' = 3,
+      'X' = 99,
+      'Y' = 70,
+      'N1' = 2,
+      'X0' = 17,
+      'Y0' = 12,
+      'Ax' = 51,
+      'By' = 48,
+      'Cx' = 34,
+      'Dy' = 36],
+     [3 > 0,
+      2 is 3 - 1,
+      pell(2, 17, 12),
+      51 is 3 * 17,
+      48 is 4 * 12,
+      99 is 51 + 48,
+      34 is 2 * 17,
+      36 is 3 * 12,
+      70 is 34 + 36]).
+step(3 > 0, builtin, [], []).
+step(2 is 3 - 1, builtin, [], []).
+step(pell(2, 17, 12),
+     rule(3),
+     ['N' = 2,
+      'X' = 17,
+      'Y' = 12,
+      'N1' = 1,
+      'X0' = 3,
+      'Y0' = 2,
+      'Ax' = 9,
+      'By' = 8,
+      'Cx' = 6,
+      'Dy' = 6],
+     [2 > 0,
+      1 is 2 - 1,
+      pell(1, 3, 2),
+      9 is 3 * 3,
+      8 is 4 * 2,
+      17 is 9 + 8,
+      6 is 2 * 3,
+      6 is 3 * 2,
+      12 is 6 + 6]).
+step(2 > 0, builtin, [], []).
+step(1 is 2 - 1, builtin, [], []).
+step(pell(1, 3, 2),
+     rule(3),
+     ['N' = 1,
+      'X' = 3,
+      'Y' = 2,
+      'N1' = 0,
+      'X0' = 1,
+      'Y0' = 0,
+      'Ax' = 3,
+      'By' = 0,
+      'Cx' = 2,
+      'Dy' = 0],
+     [1 > 0,
+      0 is 1 - 1,
+      pell(0, 1, 0),
+      3 is 3 * 1,
+      0 is 4 * 0,
+      3 is 3 + 0,
+      2 is 2 * 1,
+      0 is 3 * 0,
+      2 is 2 + 0]).
+step(1 > 0, builtin, [], []).
+step(0 is 1 - 1, builtin, [], []).
+step(pell(0, 1, 0), fact(2), [], []).
+step(3 is 3 * 1, builtin, [], []).
+step(0 is 4 * 0, builtin, [], []).
+step(3 is 3 + 0, builtin, [], []).
+step(2 is 2 * 1, builtin, [], []).
+step(0 is 3 * 0, builtin, [], []).
+step(2 is 2 + 0, builtin, [], []).
+step(9 is 3 * 3, builtin, [], []).
+step(8 is 4 * 2, builtin, [], []).
+step(17 is 9 + 8, builtin, [], []).
+step(6 is 2 * 3, builtin, [], []).
+step(6 is 3 * 2, builtin, [], []).
+step(12 is 6 + 6, builtin, [], []).
+step(51 is 3 * 17, builtin, [], []).
+step(48 is 4 * 12, builtin, [], []).
+step(99 is 51 + 48, builtin, [], []).
+step(34 is 2 * 17, builtin, [], []).
+step(36 is 3 * 12, builtin, [], []).
+step(70 is 34 + 36, builtin, [], []).
+step(297 is 3 * 99, builtin, [], []).
+step(280 is 4 * 70, builtin, [], []).
+step(577 is 297 + 280, builtin, [], []).
+step(198 is 2 * 99, builtin, [], []).
+step(210 is 3 * 70, builtin, [], []).
+step(408 is 198 + 210, builtin, [], []).
+step(1731 is 3 * 577, builtin, [], []).
+step(1632 is 4 * 408, builtin, [], []).
+step(3363 is 1731 + 1632, builtin, [], []).
+step(1154 is 2 * 577, builtin, [], []).
+step(1224 is 3 * 408, builtin, [], []).
+step(2378 is 1154 + 1224, builtin, [], []).
+step(pell_answer(solution_8, solution(665857, 470832)),
+     rule(6),
+     ['X' = 665857, 'Y' = 470832],
+     [pell(8, 665857, 470832)]).
+step(pell(8, 665857, 470832),
+     rule(3),
+     ['N' = 8,
+      'X' = 665857,
+      'Y' = 470832,
+      'N1' = 7,
+      'X0' = 114243,
+      'Y0' = 80782,
+      'Ax' = 342729,
+      'By' = 323128,
+      'Cx' = 228486,
+      'Dy' = 242346],
+     [8 > 0,
+      7 is 8 - 1,
+      pell(7, 114243, 80782),
+      342729 is 3 * 114243,
+      323128 is 4 * 80782,
+      665857 is 342729 + 323128,
+      228486 is 2 * 114243,
+      242346 is 3 * 80782,
+      470832 is 228486 + 242346]).
+step(8 > 0, builtin, [], []).
+step(7 is 8 - 1, builtin, [], []).
+step(pell(7, 114243, 80782),
+     rule(3),
+     ['N' = 7,
+      'X' = 114243,
+      'Y' = 80782,
+      'N1' = 6,
+      'X0' = 19601,
+      'Y0' = 13860,
+      'Ax' = 58803,
+      'By' = 55440,
+      'Cx' = 39202,
+      'Dy' = 41580],
+     [7 > 0,
+      6 is 7 - 1,
+      pell(6, 19601, 13860),
+      58803 is 3 * 19601,
+      55440 is 4 * 13860,
+      114243 is 58803 + 55440,
+      39202 is 2 * 19601,
+      41580 is 3 * 13860,
+      80782 is 39202 + 41580]).
+step(7 > 0, builtin, [], []).
+step(6 is 7 - 1, builtin, [], []).
+step(pell(6, 19601, 13860),
+     rule(3),
+     ['N' = 6,
+      'X' = 19601,
+      'Y' = 13860,
+      'N1' = 5,
+      'X0' = 3363,
+      'Y0' = 2378,
+      'Ax' = 10089,
+      'By' = 9512,
+      'Cx' = 6726,
+      'Dy' = 7134],
+     [6 > 0,
+      5 is 6 - 1,
+      pell(5, 3363, 2378),
+      10089 is 3 * 3363,
+      9512 is 4 * 2378,
+      19601 is 10089 + 9512,
+      6726 is 2 * 3363,
+      7134 is 3 * 2378,
+      13860 is 6726 + 7134]).
+step(6 > 0, builtin, [], []).
+step(5 is 6 - 1, builtin, [], []).
+step(10089 is 3 * 3363, builtin, [], []).
+step(9512 is 4 * 2378, builtin, [], []).
+step(19601 is 10089 + 9512, builtin, [], []).
+step(6726 is 2 * 3363, builtin, [], []).
+step(7134 is 3 * 2378, builtin, [], []).
+step(13860 is 6726 + 7134, builtin, [], []).
+step(58803 is 3 * 19601, builtin, [], []).
+step(55440 is 4 * 13860, builtin, [], []).
+step(114243 is 58803 + 55440, builtin, [], []).
+step(39202 is 2 * 19601, builtin, [], []).
+step(41580 is 3 * 13860, builtin, [], []).
+step(80782 is 39202 + 41580, builtin, [], []).
+step(342729 is 3 * 114243, builtin, [], []).
+step(323128 is 4 * 80782, builtin, [], []).
+step(665857 is 342729 + 323128, builtin, [], []).
+step(228486 is 2 * 114243, builtin, [], []).
+step(242346 is 3 * 80782, builtin, [], []).
+step(470832 is 228486 + 242346, builtin, [], []).
+step(pell_answer(check_8, true), rule(7), [], [pell_holds(8, true)]).
+step(pell_holds(8, true),
+     rule(4),
+     ['N' = 8,
+      'X' = 665857,
+      'Y' = 470832,
+      'X2' = 443365544449,
+      'Y2' = 221682772224,
+      'Twicey2' = 443365544448],
+     [pell(8, 665857, 470832),
+      443365544449 is 665857 * 665857,
+      221682772224 is 470832 * 470832,
+      443365544448 is 2 * 221682772224,
+      1 is 443365544449 - 443365544448]).
+step(443365544449 is 665857 * 665857, builtin, [], []).
+step(221682772224 is 470832 * 470832, builtin, [], []).
+step(443365544448 is 2 * 221682772224, builtin, [], []).
+step(1 is 443365544449 - 443365544448, builtin, [], []).
+step(pell_answer(y_sum_1_to_8, 568344),
+     rule(8),
+     ['Sum' = 568344],
+     [sumall(Expression, (between(1, 8, N), pell(N, _x, Expression)), 568344)]).
+step(sumall(Expression, (between(1, 8, N), pell(N, _x, Expression)), 568344), builtin, [], []).

package/examples/proof/pendulum-period.pl

@@ -0,0 +1,97 @@
+period_s(pendulum1, 2.0).
+periodError_s(pendulum1, 0.0).
+status(pendulum1, within_period_tolerance).
+reason(pendulum1, "small-angle period matches the two-second target").
+
+clause(1, constant(pi, 3.141592653589793), true).
+clause(2, experiment(pendulum1, length_m, 1.0), true).
+clause(3, experiment(pendulum1, gravity_m_s2, 9.869604401089358), true).
+clause(4, limit(pendulum1, target_period_s, 2.0), true).
+clause(5, limit(pendulum1, tolerance_s, 0.01), true).
+clause(6,
+       period(var('Experiment'), var('Period')),
+       (experiment(var('Experiment'), length_m, var('Length')),
+        experiment(var('Experiment'), gravity_m_s2, var('Gravity')),
+        var('Ratio') is var('Length') / var('Gravity'),
+        var('Root') is var('Ratio') ** 0.5,
+        constant(pi, var('Pi')),
+        var('Twopi') is 2.0 * var('Pi'),
+        var('Period') is var('Twopi') * var('Root'))).
+clause(7,
+       period_error(var('Experiment'), var('Error')),
+       (period(var('Experiment'), var('Period')),
+        limit(var('Experiment'), target_period_s, var('Target')),
+        var('Rawerror') is var('Period') - var('Target'),
+        var('Error') is abs(var('Rawerror')))).
+clause(8,
+       within_period_tolerance(var('Experiment')),
+       (period_error(var('Experiment'), var('Error')),
+        limit(var('Experiment'), tolerance_s, var('Tolerance')),
+        var('Error') < var('Tolerance'))).
+clause(9, period_s(var('Experiment'), var('Period')), period(var('Experiment'), var('Period'))).
+clause(10,
+       periodError_s(var('Experiment'), var('Error')),
+       period_error(var('Experiment'), var('Error'))).
+clause(11,
+       status(var('Experiment'), within_period_tolerance),
+       within_period_tolerance(var('Experiment'))).
+clause(12,
+       reason(var('Experiment'), "small-angle period matches the two-second target"),
+       within_period_tolerance(var('Experiment'))).
+
+step(period_s(pendulum1, 2.0),
+     rule(9),
+     ['Experiment' = pendulum1, 'Period' = 2.0],
+     [period(pendulum1, 2.0)]).
+step(period(pendulum1, 2.0),
+     rule(6),
+     ['Experiment' = pendulum1,
+      'Period' = 2.0,
+      'Length' = 1.0,
+      'Gravity' = 9.869604401089358,
+      'Ratio' = 0.10132118364233778,
+      'Root' = 0.3183098861837907,
+      'Pi' = 3.141592653589793,
+      'Twopi' = 6.283185307179586],
+     [experiment(pendulum1, length_m, 1.0),
+      experiment(pendulum1, gravity_m_s2, 9.869604401089358),
+      0.10132118364233778 is 1.0 / 9.869604401089358,
+      0.3183098861837907 is 0.10132118364233778 ** 0.5,
+      constant(pi, 3.141592653589793),
+      6.283185307179586 is 2.0 * 3.141592653589793,
+      2.0 is 6.283185307179586 * 0.3183098861837907]).
+step(experiment(pendulum1, length_m, 1.0), fact(2), [], []).
+step(experiment(pendulum1, gravity_m_s2, 9.869604401089358), fact(3), [], []).
+step(0.10132118364233778 is 1.0 / 9.869604401089358, builtin, [], []).
+step(0.3183098861837907 is 0.10132118364233778 ** 0.5, builtin, [], []).
+step(constant(pi, 3.141592653589793), fact(1), [], []).
+step(6.283185307179586 is 2.0 * 3.141592653589793, builtin, [], []).
+step(2.0 is 6.283185307179586 * 0.3183098861837907, builtin, [], []).
+step(periodError_s(pendulum1, 0.0),
+     rule(10),
+     ['Experiment' = pendulum1, 'Error' = 0.0],
+     [period_error(pendulum1, 0.0)]).
+step(period_error(pendulum1, 0.0),
+     rule(7),
+     ['Experiment' = pendulum1, 'Error' = 0.0, 'Period' = 2.0, 'Target' = 2.0, 'Rawerror' = 0.0],
+     [period(pendulum1, 2.0),
+      limit(pendulum1, target_period_s, 2.0),
+      0.0 is 2.0 - 2.0,
+      0.0 is abs(0.0)]).
+step(limit(pendulum1, target_period_s, 2.0), fact(4), [], []).
+step(0.0 is 2.0 - 2.0, builtin, [], []).
+step(0.0 is abs(0.0), builtin, [], []).
+step(status(pendulum1, within_period_tolerance),
+     rule(11),
+     ['Experiment' = pendulum1],
+     [within_period_tolerance(pendulum1)]).
+step(within_period_tolerance(pendulum1),
+     rule(8),
+     ['Experiment' = pendulum1, 'Error' = 0.0, 'Tolerance' = 0.01],
+     [period_error(pendulum1, 0.0), limit(pendulum1, tolerance_s, 0.01), 0.0 < 0.01]).
+step(limit(pendulum1, tolerance_s, 0.01), fact(5), [], []).
+step(0.0 < 0.01, builtin, [], []).
+step(reason(pendulum1, "small-angle period matches the two-second target"),
+     rule(12),
+     ['Experiment' = pendulum1],
+     [within_period_tolerance(pendulum1)]).