>>593
(1)をやってみたけど問題文の値と結果が全然異なるから違ってたらスマン
あとぼるじょあをよろしく

#include <iostream>
#include <cmath>
#define sgn(a)((a) == 0.0 ? 0.0 : (a) < 0.0 ? -1.0 : 1.0)
//f(x) = (x^2-x+1)^2 - x^2(2x+log x)
double func(double x)
{
return pow((x*x - x + 1.0), 2) - x*x * (2*x + log(x));
}
//f'(x) ≒ f(x+ε)-f(x-ε) / 2ε
double func2(double x)
{
double eps = pow(10.0, -5) * sgn(x)*x;
return func(x + eps) - func(x - eps) / 2*eps;
}
int main()
{
using namespace std;

cout << " x f'(x) true f'(x)" << endl;
for (double x = 0.1; x < 4.1; x += 0.1) {
printf("%.1f %10.5f %10.5f¥n", x, func(x), func2(x));
}
return 0;
}