Berichtgenerator - Berichte generieren mit Hilfe von Gforth — ein Kinderspiel; 3 Listings: in.fs m.dep bio.txt
GTK-Listings - Widgets zum Anfassen - GUI-Skripting mit Forth und GTK+; Listing 1+2.
Quellcode Bibliothek - Ein einfaches Quellcode–Bibiliothekssystem; 3 Listings: SAMPLE.LIB APPL.F LIBRARY.ANS
\ Euler 9 .solution ( mka)
decimal
vocabulary euler
euler definitions
\ Satz des Pythagoras: a^2 + b^2 = c^2
: pyt ( a b -- c^2 )
dup * swap dup * + ;
\ Quadratwurzel ziehen: Heron-Verfahren (Heron von Alexandria)
\ Die Iterationsvorschrift lautet: xn+1 = ( xn + a/xn ) / 2
\ In RPN: a xn / xn + 2 / --> xn+1
\ Weil nur die natürlichen Werte der Wurzel aus c^2 gebraucht werden,
\ wird ganzzahlig kalkuliert.
\ Quelle: forth-ev wiki
\ http://www.forth-ev.de/wiki/doku.php
: sqrt ( a -- s )
DUP 0= ?EXIT
DUP >R
BEGIN
R@ OVER / OVER + 2/
DUP ROT - ABS 1 <=
UNTIL R> DROP ;
\ Wir brauchen nur die Werte von c, für die c eine natürliche Zahl ist.
: nsqrt? ( a b -- f )
pyt
dup sqrt ( -- n s )
dup * ( -- n s*s )
- 0= ; ( -- f ) \ if zero, s is natural number.
\ Ist a+b+c=1000?
: 1000? ( a b -- f )
2dup pyt sqrt
+ + 1000 = ;
\ Huebsche Ausgabe der gefundenen Loesungen (!)
: .solution ( a b -- )
cr ." -------------------"
cr 2dup swap ." a=" . ." b=" .
2dup pyt sqrt ." c=" .
cr 2dup 2dup pyt sqrt + + ." a+b+c=" .
cr 2dup pyt sqrt * * ." abc=" . ;
\ Teste a,b bis zum Wert a=b=xx
variable xx 1000 xx !
: p9 ( -- )
xx @ 0 DO
xx @ 0 DO
i j nsqrt? IF
i j 1000? IF i j .solution THEN
THEN
LOOP
LOOP
cr ." -------------------" ;
words cr cr
( finis)