fn rationalApproxExtended80(z: f80) f80 z: f80214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231fn rationalApproxExtended80(z: f80) f80 {
const pS0: f80 = 1.66666666666666666631e-01;
const pS1: f80 = -4.16313987993683104320e-01;
const pS2: f80 = 3.69068046323246813704e-01;
const pS3: f80 = -1.36213932016738603108e-01;
const pS4: f80 = 1.78324189708471965733e-02;
const pS5: f80 = -2.19216428382605211588e-04;
const pS6: f80 = -7.10526623669075243183e-06;
const qS1: f80 = -2.94788392796209867269e+00;
const qS2: f80 = 3.27309890266528636716e+00;
const qS3: f80 = -1.68285799854822427013e+00;
const qS4: f80 = 3.90699412641738801874e-01;
const qS5: f80 = -3.14365703596053263322e-02;
const p = z * (pS0 + z * (pS1 + z * (pS2 + z * (pS3 + z * (pS4 + z * (pS5 + z * pS6))))));
const q = 1.0 + z * (qS1 + z * (qS2 + z * (qS3 + z * (qS4 + z * qS5))));
return p / q;
}