C++實現二分法求連續一元函式根
阿新 • • 發佈:2020-06-16
本文例項為大家分享了C++實現二分法求連續一元函式根的具體程式碼,供大家參考,具體內容如下
設計一個用二分法求連續一元函式根的通用函式solve
此函式有三個引數:
- 第一個是函式指標,指向所要求根的連續函式
- 第二、三個引數指出根的區間,且確保函式在區間的兩個端點異號
函式的返回值為求得的解
要求編寫main函式如下:
double fun(double x) { double y; y=4*pow(x,3)-6*pow(x,2)+3*x-2; return y; } int main() { cout<<"4*x^3-6*x^2+3*x-2=0在區間(1,2)的根為 x="<<solve(fun,1,2); return 0; }
C++實現:
#include <iostream> #include <cmath> using namespace std; double solve(double (*fun)(double x),double a,double b); double fun(double x); int main() { cout << "4*x^3-6*x^2+3*x-2=0在區間(1,2)的根為 x=" << solve(fun,2); return 0; } double solve(double (*fun)(double x),double b) { double i = b - a; double c = (a + b) / 2; while (i > 0.0000001) { i = b - a; if (fun(c) == 0)return c; if (fun(c) * fun(a) < 0) { b = c; c = (a + b) / 2; } else { a = c; c = (a + b) / 2; } } return c; } double fun(double x) { double y; y = 4 * pow(x,3) - 6 * pow(x,2) + 3 * x - 2; return y; }
總結:
- 函式與指標的結合
- 注意返回的型別與要求
以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支援我們。