>>63
ごめん、式間違ってた
math ライブラリの fabs も使わないようにしたバージョン
#include <stdio.h>
/* ニュートン法による解法 */
int main()
{
double num = 1729.3, cbrt, eps = 1e-10;
cbrt = num;
while(cbrt*cbrt*cbrt - num > eps || cbrt*cbrt*cbrt - num < -eps){
cbrt = (num + 2*cbrt*cbrt*cbrt) / 3 / cbrt / cbrt;
printf("%25.15f\n", cbrt);
}
printf("%f^(1/3) = %.15f\n", num, cbrt);
return 0;
}
eps が 1e-12 あたりまでなら収束する