Kurs:Computational Physics/Pseudocode:Differentialquotient

// programm numdiff tabelliert den sin'(x) am Punkt x=1 als Funktion der Schrittweite
programm numdiff
for h = 1 to 16 step 1
h= $10^{h}$
DQ= cos(x); #exakter Differentialquotient
DQ1= (sin(x+h)-sin(x))/h ; // asymmetrische Formel
DQ2= (sin(x+h)-(sin(x-h))/2h; // symmetrische Formel
y= [h; (DQ1-DQ)/DQ; (DQ2-DQ)/DQ]; //Vergleichende Tabelle
next
output "Schrittweite h", "Asymmetrische Formel", "Symmetrische Formel"
//Genauigkeit der Variablen codieren in der jeweiligen Programmiersprache !
output y
end numdiff