| nMax = 2; beta = SetAccuracy1, 2000; |
| gamma = SetAccuracy0, 2000; |
| Forj = 1, j <= NumberOfGuesses, j++, |
| {x = SetAccuracyinitialValuesj, 2000, |
| Do |
| fx = SetAccuracyfx, 2000; |
| k = SetAccuracyx + beta fx, 2000; |
| fk = SetAccuracyfk, 2000; |
| FD = SetAccuracy(fk - fx)/(beta fx), 2000; |
| y = SetAccuracyx - fx/FD, 2000; |
| fy = SetAccuracyfy, 2000; |
| z = SetAccuracyy - (fy/FD)*(1 + ((2 + beta FD)/(1 |
| + beta FD)) (fy/fx)), 2000; |
| fz = SetAccuracyfz, 2000; |
| FDxz = SetAccuracy(fz - fx)/(z - x), 2000; |
| FDxy = SetAccuracy(fx - fy)/(x - y), 2000; |
| FDxk = SetAccuracy(fk - fx)/(k - x), 2000; |
| x = SetAccuracyz - fz/(FDxz + ((FDxk - FDxy)/(k - y) |
| - (FDxk - FDxz)/(k - z) |
| - (FDxy - FDxz)/(y - z))*(x - z) |
| + gamma (z - x)*(z - k)*(z - y)), 2000; |
| , {n, nMax; |
| PrintColumn |
| Nx, 128, |
| Nfx |
| , Background -> LightGreen, LightGray, |
| Frame -> True]]; |
| } |
| ]; |