摘要 讨论 acosPi(x) = acos(x)/π 的错误计算。
由计算机的错误计算(七十八)知,IEEE 754-2019中也含有函数 acosPi(x) = acos(x)/π . 因此,本节探讨该函数的计算精度问题。
例1. 计算 acosPi(0.9999999999) = acos(0.9999999999)/π .
类似于(七十八)中的 asinPi(x)函数,似乎没有哪种语言内置了函数 acosPi(x) . 那么,不妨用Java直接编程计算吧:
import java.lang.Math;
public class acosPi{
public static void main(String[] args) {
double x = 0.9999999999;
double result = Math.acos(x)/Math.PI;
System.out.println(result);
}
}
则运行后输出 4.501581767054305E-6 .
另外,若在Go语言下直接编程计算:
package main
import (
"fmt"
"math"
)
func main() {
x := 0.9999999999
value := math.Acos(x)/math.Pi
fmt.Printf("%v", value)
}
则输出 4.501581767147819e-06 .
从上可知,二者的后6位数字不相同。这时,是否以为二者的前10位数字是正确的?
事实上,正确值是 0.4501581580823044e-5(ISRealsoft 提供)。因此,二者均只有7位正确数字,错误率均为 9/16 = 56.25% .