199 Real
zero(0), one(1), two(2), p1(0.1);
207 std::cout << std::right << std::setw(20) <<
"CHECK REGRET: v'(eps) is correct? \n";
208 std::cout << std::right << std::setw(20) <<
"t"
209 << std::setw(20) <<
"v'(x)"
210 << std::setw(20) <<
"(v(x+t)-v(x-t))/2t"
211 << std::setw(20) <<
"Error"
213 for (
int i = 0; i < 13; i++) {
216 diff = (vy-vx)/(two*t);
217 err = std::abs(diff-dv);
218 std::cout << std::scientific << std::setprecision(11) << std::right
219 << std::setw(20) << t
220 << std::setw(20) << dv
221 << std::setw(20) << diff
222 << std::setw(20) << err
234 std::cout << std::right << std::setw(20) <<
"CHECK REGRET: v''(eps) is correct? \n";
235 std::cout << std::right << std::setw(20) <<
"t"
236 << std::setw(20) <<
"v''(x)"
237 << std::setw(20) <<
"(v'(x+t)-v'(x-t))/2t"
238 << std::setw(20) <<
"Error"
240 for (
int i = 0; i < 13; i++) {
243 diff = (vy-vx)/(two*t);
244 err = std::abs(diff-dv);
245 std::cout << std::scientific << std::setprecision(11) << std::right
246 << std::setw(20) << t
247 << std::setw(20) << dv
248 << std::setw(20) << diff
249 << std::setw(20) << err
262 std::cout << std::right << std::setw(20) <<
"CHECK REGRET: v'(0) is correct? \n";
263 std::cout << std::right << std::setw(20) <<
"t"
264 << std::setw(20) <<
"v'(x)"
265 << std::setw(20) <<
"(v(x+t)-v(x-t))/2t"
266 << std::setw(20) <<
"Error"
268 for (
int i = 0; i < 13; i++) {
271 diff = (vy-vx)/(two*t);
272 err = std::abs(diff-dv);
273 std::cout << std::scientific << std::setprecision(11) << std::right
274 << std::setw(20) << t
275 << std::setw(20) << dv
276 << std::setw(20) << diff
277 << std::setw(20) << err
289 std::cout << std::right << std::setw(20) <<
"CHECK REGRET: v''(0) is correct? \n";
290 std::cout << std::right << std::setw(20) <<
"t"
291 << std::setw(20) <<
"v''(x)"
292 << std::setw(20) <<
"(v'(x+t)-v'(x-t))/2t"
293 << std::setw(20) <<
"Error"
295 for (
int i = 0; i < 13; i++) {
298 diff = (vy-vx)/(two*t);
299 err = std::abs(diff-dv);
300 std::cout << std::scientific << std::setprecision(11) << std::right
301 << std::setw(20) << t
302 << std::setw(20) << dv
303 << std::setw(20) << diff
304 << std::setw(20) << err
317 std::cout << std::right << std::setw(20) <<
"CHECK REGRET: v'(-eps) is correct? \n";
318 std::cout << std::right << std::setw(20) <<
"t"
319 << std::setw(20) <<
"v'(x)"
320 << std::setw(20) <<
"(v(x+t)-v(x-t))/2t"
321 << std::setw(20) <<
"Error"
323 for (
int i = 0; i < 13; i++) {
326 diff = (vy-vx)/(two*t);
327 err = std::abs(diff-dv);
328 std::cout << std::scientific << std::setprecision(11) << std::right
329 << std::setw(20) << t
330 << std::setw(20) << dv
331 << std::setw(20) << diff
332 << std::setw(20) << err
344 std::cout << std::right << std::setw(20) <<
"CHECK REGRET: v''(-eps) is correct? \n";
345 std::cout << std::right << std::setw(20) <<
"t"
346 << std::setw(20) <<
"v''(x)"
347 << std::setw(20) <<
"(v'(x+t)-v'(x-t))/2t"
348 << std::setw(20) <<
"Error"
350 for (
int i = 0; i < 13; i++) {
353 diff = (vy-vx)/(two*t);
354 err = std::abs(diff-dv);
355 std::cout << std::scientific << std::setprecision(11) << std::right
356 << std::setw(20) << t
357 << std::setw(20) << dv
358 << std::setw(20) << diff
359 << std::setw(20) << err