形式
Plain text
投稿日時
2019-01-15 23:40
公開期間
無期限
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  4. #define EPS 1.0e-8
  5. #define POW 2
  6. double Newton(double);
  7. double Func(double, double);
  8. double dFunc(double);
  9. int main(){
  10. double Ln = Newton(2.0);//L_1 = sqrt(2)
  11. printf("%12.10f\n",Ln);
  12. int n = 1;
  13. int pow_base = 2;
  14. int n_max = 20;
  15. while( n < n_max ){
  16. printf( "%12.10f\n", pow_base * Ln );
  17. Ln = Newton( 2-Newton(4 - pow(Ln, 2)) );
  18. n += 1;
  19. pow_base *= POW;
  20. }
  21. return 0;
  22. }
  23. double Newton(double C){
  24. double x = C;
  25. double dx = -Func(x, C) / dFunc(x);
  26. while( fabs(dx) > EPS ){
  27. x += dx;
  28. dx = -Func(x, C) / dFunc(x);
  29. }
  30. return x;
  31. }
  32. double Func(double x, double C){
  33. return x*x-C;
  34. }
  35. double dFunc(double x){
  36. return 2.0 * x;
  37. }
ダウンロード 印刷用表示

このコピペの URL

JavaScript での埋め込み

iframe での埋め込み

元のテキスト