質問編集履歴

2

2020/05/24 07:49

投稿

hacch
hacch

スコア15

test CHANGED
File without changes
test CHANGED
@@ -105,3 +105,97 @@
105
105
 
106
106
 
107
107
  ```
108
+
109
+
110
+
111
+ ```C
112
+
113
+
114
+
115
+ #include<stdio.h>
116
+
117
+ #include<math.h>
118
+
119
+ double f(double x);
120
+
121
+ double fd(double x);
122
+
123
+ double power(double x, int n);
124
+
125
+ int main(void){
126
+
127
+ double eps=1.0e-10;
128
+
129
+ double a=2.0;
130
+
131
+ double x;
132
+
133
+ int k;
134
+
135
+ x=a;
136
+
137
+ for(k=0;fabs(f(x)/fd(x))>eps;k++){
138
+
139
+ x=x-f(x)/fd(x);
140
+
141
+ printf("k=%d x=%.7f\n",k,x);
142
+
143
+ }
144
+
145
+ printf("ans=% .7f \n",x);
146
+
147
+ return 0;
148
+
149
+ }
150
+
151
+ double f(double x){
152
+
153
+ double x6,x4,x3,y;
154
+
155
+ x3=power(x,3);
156
+
157
+ x4=x3*x;
158
+
159
+ x6=x3*x3;
160
+
161
+ y=x6-7*x4+11*x3-10;
162
+
163
+ return y;
164
+
165
+ }
166
+
167
+ double fd(double x){
168
+
169
+ double x5,x3,x2,y;
170
+
171
+ x2=power(x,2);
172
+
173
+ x3=x2*x;
174
+
175
+ x5=x2*x3;
176
+
177
+
178
+
179
+ y=6*x+6*x+6*x+6*x+6*x
180
+
181
+ -28*x+(-28*x)+(-28*x)
182
+
183
+ +33*x+33*x;
184
+
185
+ return y;
186
+
187
+ }
188
+
189
+
190
+
191
+ double power(double x,int n){
192
+
193
+ if(n==0) return 1;
194
+
195
+ else return x*power(x,n-1);
196
+
197
+ }
198
+
199
+
200
+
201
+ ```

1

2020/05/24 07:49

投稿

hacch
hacch

スコア15

test CHANGED
File without changes
test CHANGED
@@ -11,6 +11,12 @@
11
11
  2箇所分からないところがコメントしてあります。
12
12
 
13
13
  分かる方いらっしゃいましたら教えていただけると嬉しいです。。
14
+
15
+  
16
+
17
+ ①は式はあっているはずですが、不安なので質問しました。
18
+
19
+ ②はニュートン法のために返却値として返したいですが、どのような式が入るのかよく分からなかったので、ここが出来ないと①が合っているかも断定できなかったので質問しました。
14
20
 
15
21
 
16
22