home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / l / l217 / 2.ddi / EXAMPLES / CH11EX01.PRO < prev    next >
Encoding:
Prolog Source  |  1990-03-26  |  808 b   |  29 lines
  1. /*
  2.    Copyright (c) 1986, 90 by Prolog Development Center
  3. */
  4.  
  5. predicates
  6.      solve(real, real, real)
  7.      reply(real, real, real)
  8.      mysqrt(real, real, real)
  9.      equal(real, real)
  10.  
  11. clauses
  12.      solve(A, B, C) :- D = B*B-4*A*C,  reply(A, B, D), nl.
  13.      reply(_, _, D) :- D < 0, write("No solution"), !.
  14.      reply(A, B, D) :- D = 0, X=-B/(2*A), write("x=", X), !.
  15.      reply(A, B, D) :- 
  16.          mysqrt(D, D, SqrtD),
  17.          X1 = (-B + SqrtD)/(2*A),
  18.          X2 = (-B - SqrtD)/(2*A),
  19.          write("x1 = ", X1," and x2 = ", X2).
  20.  
  21.      mysqrt(X, Guess, Root) :-
  22.          NewGuess = Guess-(Guess*Guess-X)/2/Guess,
  23.          not(equal(NewGuess, Guess)), !,
  24.          mysqrt(X, NewGuess, Root).
  25.  
  26.      mysqrt(_, Guess, Guess).
  27.  
  28.      equal(X, Y) :- X/Y > 0.99999 , X/Y < 1.00001.
  29.