牛顿法解方程

题目要求：

应用牛顿法解方程：

#include "stdio.h"#include "math.h"float func(float x){    return pow(x,4)+pow(x,3)+1;}float x(float a,float b,int k){    return a+k*(b-a)/4  ;}float ING(float a,float b){    return ((b-a)/90)*(7*func(x(a,b,0))+32*func(x(a,b,1))+    12*func(x(a,b,2))+32*func(x(a,b,3))+7*func(x(a,b,4)));}main(){    float a,b;    printf("Please input the low & high limitation and the accuracy\n");    printf("Low limitation:");    scanf("%f",&a);    printf("High limitation:");    scanf("%f",&b);    printf("The result of integration is %f",ING(a,b));    getche();}

运行结果：

运行结果

欧拉方法求解微分方程

题目要求：

已知微分方程的初值问题如下：

求y(1)的值。

#include "stdio.h"#include "math.h"double SQRT(double a){				/*迭代法开方运算*/    double xx = a,x = 0.0; 			/*迭代初值*/    while(fabs(xx - x)>0.00001){        x = xx;        xx = 0.5*(x + a / x) ;    }    return xx;						/*返回迭代结果*/}main(){    double a,r;    printf("Please input a digit for squart\n");    scanf("%lf",&a);    r = SQRT(a);    printf("Sqrt(%f) = %f\n",a,r);    getche();}

运行结果：

运行结果