質問編集履歴
2
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -51,4 +51,51 @@
|
|
|
51
51
|
else return x*power(x,n-1);
|
|
52
52
|
}
|
|
53
53
|
|
|
54
|
+
```
|
|
55
|
+
|
|
56
|
+
```C
|
|
57
|
+
|
|
58
|
+
#include<stdio.h>
|
|
59
|
+
#include<math.h>
|
|
60
|
+
double f(double x);
|
|
61
|
+
double fd(double x);
|
|
62
|
+
double power(double x, int n);
|
|
63
|
+
int main(void){
|
|
64
|
+
double eps=1.0e-10;
|
|
65
|
+
double a=2.0;
|
|
66
|
+
double x;
|
|
67
|
+
int k;
|
|
68
|
+
x=a;
|
|
69
|
+
for(k=0;fabs(f(x)/fd(x))>eps;k++){
|
|
70
|
+
x=x-f(x)/fd(x);
|
|
71
|
+
printf("k=%d x=%.7f\n",k,x);
|
|
72
|
+
}
|
|
73
|
+
printf("ans=% .7f \n",x);
|
|
74
|
+
return 0;
|
|
75
|
+
}
|
|
76
|
+
double f(double x){
|
|
77
|
+
double x6,x4,x3,y;
|
|
78
|
+
x3=power(x,3);
|
|
79
|
+
x4=x3*x;
|
|
80
|
+
x6=x3*x3;
|
|
81
|
+
y=x6-7*x4+11*x3-10;
|
|
82
|
+
return y;
|
|
83
|
+
}
|
|
84
|
+
double fd(double x){
|
|
85
|
+
double x5,x3,x2,y;
|
|
86
|
+
x2=power(x,2);
|
|
87
|
+
x3=x2*x;
|
|
88
|
+
x5=x2*x3;
|
|
89
|
+
|
|
90
|
+
y=6*x+6*x+6*x+6*x+6*x
|
|
91
|
+
-28*x+(-28*x)+(-28*x)
|
|
92
|
+
+33*x+33*x;
|
|
93
|
+
return y;
|
|
94
|
+
}
|
|
95
|
+
|
|
96
|
+
double power(double x,int n){
|
|
97
|
+
if(n==0) return 1;
|
|
98
|
+
else return x*power(x,n-1);
|
|
99
|
+
}
|
|
100
|
+
|
|
54
101
|
```
|
1
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -5,6 +5,9 @@
|
|
|
5
5
|
プログラムはこちらです。
|
|
6
6
|
2箇所分からないところがコメントしてあります。
|
|
7
7
|
分かる方いらっしゃいましたら教えていただけると嬉しいです。。
|
|
8
|
+
|
|
9
|
+
①は式はあっているはずですが、不安なので質問しました。
|
|
10
|
+
②はニュートン法のために返却値として返したいですが、どのような式が入るのかよく分からなかったので、ここが出来ないと①が合っているかも断定できなかったので質問しました。
|
|
8
11
|
|
|
9
12
|
```C
|
|
10
13
|
|