program opnieuw; { This software has been designed and is CopyLefted by Han de Bruijn: (===) @-O^O-@ #/_\# ### Least Squares best fit to a Hyperbola y = B/(x-A) of the Year/Light-speed data from relevant tables in * The Atomic Constants, Light, and Time * http://www.setterfield.org/report/report.html } uses Math, tabellen, rekenen; procedure Gewichten(var w : lijst); var k,L : integer; begin L := Length(w); for k := 0 to L-1 do begin w[k] := 1; end; end; procedure Format; const delta : double = 1; { 1 2 3 4 5 6 7 8 9 0 1 2 3 } reeks : array[1..13] of integer = (1,2,4,5,6,7,8,10,12,13,14,27,29); var licht,jaar,w : lijst; keer,maal : integer; A,B,R,fout : double; uit : TextFile; kop,rood,door : string; eps : punt; begin AssignFile(uit,'tabel.htm'); Rewrite(uit); Writeln(uit,''); Writeln(uit,''); Writeln(uit,'
'); kop := ' '; Writeln(uit,kop); Writeln(uit,' '); Writeln(uit,' '); Writeln(uit,' '); Writeln(uit, ' '); Writeln(uit, ' '); Writeln(uit, ' '); Writeln(uit, ' '); Writeln(uit, ' '); Writeln(uit, ' '); Writeln(uit,' '); Writeln(uit,' '); for keer := 1 to 13 do begin SetLength(jaar,0); SetLength(licht,0); Gegevens(keer,jaar,licht); SetLength(w,Length(jaar)); Gewichten(w); Writeln(uit,' '); Writeln(uit,' '); uitkijken(jaar,licht,A,B); rood := ''; door := ''; if A < 0 then rood := ''; if A < 0 then door := ''; Writeln(uit,' '); kijkuit(jaar,licht,A,B,R); rood := ''; door := ''; if A < 0 then rood := ''; if A < 0 then door := ''; Writeln(uit,' '); Jacquelin(jaar,licht,A,B); rood := ''; door := ''; if A < 0 then rood := ''; if A < 0 then door := ''; Writeln(uit,' '); Raaklijn(jaar,licht,w,A,B,R); rood := ''; door := ''; if A < 0 then rood := ''; if A < 0 then door := ''; Writeln(uit,' '); uitkijken(jaar,licht,A,B); if A < 0 then begin rood := ''; door := ''; for maal := 0 to 20 do begin NewtonHdB(jaar,licht,A,B,eps); { Writeln(eps.x,eps.y); } if (eps.x < delta) then begin rood := ''; door := ''; Break; end; end; Writeln(uit,' '); end else begin Writeln(uit,' '); end; uitkijken(jaar,licht,A,B); if A < 0 then begin rood := ''; door := ''; for maal := 0 to 20 do begin Leibovici(jaar,licht,A,B,fout); { Writeln(fout); } if (fout < delta) then begin rood := ''; door := ''; Break; end; end; Writeln(uit,' '); end else begin Writeln(uit,' '); end; Writeln(uit,' '); Writeln(uit,' '); end; Writeln(uit,' '); Writeln(uit,'
tabelmethode 1methode 2methode 3methode 4methode 5methode 6grafiek
',keer,'.txt'+rood,A:10:0,door+''+rood, A:10:0,' ± ',R:5:0,door+''+rood,A:10:0,door+''+rood, A:10:0,' ± ',R:5:0,door+''+rood,A:10:0,door+''+rood,A:10:0,door+'',keer,'.jpg
'); Writeln(uit,''); Writeln(uit,''); CloseFile(uit); end; begin Format; end.