- 234 名前：デフォルトの名無しさん mailto:sage [2009/02/03(火) 14:58:41 ]
- >>233
オイラー法ってこんなか？
#include <stdio.h>
#include <math.h>
double newton(){
double x=1.0;
int i;
for(i=0;i<50;i++){
x=x-(x*x-2)/(2*x);
}
return x;
}
double euler(){
int n=50;
double x=1.0,h=(double)(2-1)/n;
int i;
for(i=0;i<n;i++){
x=x+h/(2*x);
}
return x;
}
int main(){
double temp=newton();
printf("sqrt(2)：%f,newton:%f,差:%f\n",sqrt(2),temp,sqrt(2)-temp);
temp=euler();
printf("sqrt(2)：%f,euler:%f,差:%f\n",sqrt(2),temp,sqrt(2)-temp);
return 0;
}
|
 |
