8 lines
296 B
Text
8 lines
296 B
Text
fcn rump(a,b){ b=BN(b);
|
|
b2,b4,b6,b8:=b.pow(2),b.pow(4),b.pow(6),b.pow(8);
|
|
a2:=BN(a).pow(2);
|
|
r:=( b6*33375 + a2*(a2*b2*11 - b6 - b4*121 - 2)*100 + b8*550 )*ten2n;
|
|
r+=BN(a)*ten2n*100/(2*b);
|
|
tostr(r,66,32)
|
|
}
|
|
println("\n3rd: Rump's example: f(77617.0, 33096.0) = ",rump(77617,33096));
|