Changes between Version 5 and Version 6 of doc/maths/remez/tutorialfixingparameters
 Timestamp:
 Jan 9, 2012, 11:38:19 PM (11 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

doc/maths/remez/tutorialfixingparameters
v5 v6 13 13 {{{ 14 14 #!cpp 15 const double a1 = 9.99999976 6327004083636926544075452569800e1;15 const double a1 = 9.999999765898820673279342160490060830302e1; 16 16 }}} 17 17 … … 68 68 { 69 69 RemezSolver<3, real> solver; 70 solver.Run("1e1000", real::R_PI_2 * real::R_PI_2, f, g, 30);70 solver.Run("1e1000", real::R_PI_2 * real::R_PI_2, f, g, 40); 71 71 return 0; 72 72 } … … 87 87 88 88 {{{ 89 Final error: 4.613764035583927492630791542066631050746e989 Step 8 error: 4.618689007546850899022101933442449327546e9 90 90 Polynomial estimate: 91 x**0*1.66666570 8792613295297026208726412532085e192 +x**1*8.333017 134192478312472752663154642556843e393 +x**2*1.98066 0618990173078599247774582918018345e494 +x**3*2.6000 38299369332084157603784882057661080e691 x**0*1.666665709650470145824129400050267289858e1 92 +x**1*8.333017291562218127986291618761571373087e3 93 +x**2*1.980661520135080504411629636078917643846e4 94 +x**3*2.600054767890361277123254766503271638682e6 95 95 }}} 96 96 … … 101 101 double fastsin2(double x) 102 102 { 103 const double a3 = 1.66666570 8792613295297026208726412532085e1;104 const double a5 = 8.333017 134192478312472752663154642556843e3;105 const double a7 = 1.98066 0618990173078599247774582918018345e4;106 const double a9 = 2.6000 38299369332084157603784882057661080e6;103 const double a3 = 1.666665709650470145824129400050267289858e1; 104 const double a5 = 8.333017291562218127986291618761571373087e3; 105 const double a7 = 1.980661520135080504411629636078917643846e4; 106 const double a9 = 2.600054767890361277123254766503271638682e6; 107 107 108 108 return x + x*x*x * (a3 + x*x * (a5 + x*x * (a7 + x*x * a9))));