Octave

Rumpの例題

Rumpの例題 a = 77617, b = 33096; に対して f = (333.75 - a^2) * b^6 + a^2 * (11 * a^2 * b^2 -121 * b^2 -121 * b^4 - 2) + 5.5 * b^8 + a / (2 * b) を計算する. fの正しい結果は f = -0.827396... 一方で,これらをOctaveやFortranで計算してみると..…