评价网 > 科教论文 > 理学毕业论文 > 数学论文 > 正文

二分法求解单变量非线性方程及其应用与实现(3)

2013-07-15 01:12
end
max1 = 1 + round((log(b - a) - log(delta))/log(2));
for k = 1:max1
    c = (a + b)/2;
    yc = fevel(’f’, c);
    if yc == 0 a = c; b = c; break,
    elseif yb * yc > 0
        b = c; yb = yc;
    else
        a = c; ya = c;
    end
    if (b - a) < delta, break
    end
end
k, c = (a + b)/2, err = abs(b - a), yc = feval(‘f’, c)

6.2 %基本二分法的C语言实现
　　方程式为：f(x) = 0，示例中f(x) = 1+x-x^3 
　　使用示例： 
　　input a b e: 1 2 1e-5 
　　solution: 1.32472 
　　源码如下： 
　　#include <stdio.h> 
　　#include <stdlib.h> 
　　#include <math.h> 
　　#include <assert.h> 
　　double f(double x) 
　　{ 
　　return 1+x-x*x*x; 
　　} 
　　int main() 
　　{ 
　　double a = 0, b = 0, e = 1e-5; 
　　printf("input a b e: "); 
　　scanf("%lf%lf%lf", &a, &b, &e); 
　　e = fabs(e); 
　　if (fabs(f(a)) <= e) 
　　{ 
　　printf("solution: %lg\n", a); 
　　} 
　　else if (fabs(f(b)) <= e) 
　　{ 
　　printf("solution: %lg\n", b); 
　　} 
　　else if (f(a)*f(b) > 0) 
　　{ 
　　printf("f(%lg)*f(%lg) > 0 ! need <= 0 !\n", a, b); 
　　} 
　　else 
　　{ 
　　while (fabs(b-a) > e)  （转载自http://zw.nseac.coM科教作文网）
　　{ 
　　double c = (a+b)/2.0; 
　　if (f(a)* f ( c ) < 0) 
　　b = c; 
　　else 
　　a = c; 
　　} 
　　printf("solution: %lg\n", (a+b)/2.0);