In this tutorial you will get program for bisection method in C and C++.
To find a root very accurately Bisection Method is used in Mathematics. Bisection method algorithm is very easy to program and it always converges which means it always finds root.
Bisection Method repeatedly bisects an interval and then selects a subinterval in which root lies. It is a very simple and robust method but slower than other methods.
It is also called Interval halving, binary search method and dichotomy method.
Bisection Method calculates the root by first calculating the mid point of the given interval end points.
Bisection Method Procedure
The input for the method is a continuous function f, an interval [a, b], and the function values f(a) and f(b). The function values are of opposite sign (there is at least one zero crossing within the interval). Each iteration performs these steps:
1. Calculate the midpoint c = (a + b)/2
2. Calculate the function value at the midpoint, function(c).
3. If convergence is satisfactory (that is, a – c is sufficiently small, or f(c) is sufficiently small), return c and stop iterating.
4. Examine the sign of f(c) and replace either (a, f(a)) or (b, f(b)) with (c, f(c)) so that there is a zero crossing within the new interval.
Pros and Cons
Advantage of the bisection method is that it is guaranteed to be converged and very easy to implement.
Disadvantage of bisection method is that it cannot detect multiple roots and is slower compared to other methods of calculating the roots.
Program for Bisection Method in C
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 |
#include<stdio.h> //function used is x^3-2x^2+3 double func(double x) { return x*x*x - 2*x*x + 3; } double e=0.01; double c; void bisection(double a,double b) { if(func(a) * func(b) >= 0) { printf("Incorrect a and b"); return; } c = a; while ((b-a) >= e) { c = (a+b)/2; if (func(c) == 0.0){ printf("Root = %lf\n",c); break; } else if (func(c)*func(a) < 0){ printf("Root = %lf\n",c); b = c; } else{ printf("Root = %lf\n",c); a = c; } } } int main() { double a,b; a=-10; b=20; printf("The function used is x^3-2x^2+3\n"); printf("a = %lf\n",a); printf("b = %lf\n",b); bisection(a,b); printf("\n"); printf("Accurate Root calculated is = %lf\n",c); return 0; } |
Output
a = -10.000000
b = 20.000000
Root = 5.000000
Root = -2.500000
Root = 1.250000
Root = -0.625000
Root = -1.562500
Root = -1.093750
Root = -0.859375
Root = -0.976563
Root = -1.035156
Root = -1.005859
Root = -0.991211
Root = -0.998535
Accurate Root calculated is = -0.998535
Program for Bisection Method in C++
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 |
#include<iostream> using namespace std; //function used is x^3-2x^2+3 double func(double x) { return x*x*x - 2*x*x + 3; } double e=0.01; double c; void bisection(double a,double b) { if(func(a) * func(b) >= 0) { cout<<"Incorrect a and b"; return; } c = a; while ((b-a) >= e) { c = (a+b)/2; if (func(c) == 0.0){ cout << "Root = " << c<<endl; break; } else if (func(c)*func(a) < 0){ cout << "Root = " << c<<endl; b = c; } else{ cout << "Root = " << c<<endl; a = c; } } } int main() { double a,b; a=-10; b=20; cout<<"The function used is x^3-2x^2+3\n"; cout<<"a = "<<a<<endl; cout<<"b = "<<b<<endl; bisection(a,b); cout<<"\n"; cout<<"Accurate Root calculated is = "<<c<<endl; return 0; } |
This article is submitted by Rahul Maheshwari. You can connect with him on facebook.
Comment below if you have any queries regarding above program for bisection method in C and C++.
What is using namespace Std; ?
Using Using namespaces used to compile cout, cin, Endl
cos typing std:: every line is so annoying and hussle.
This program is be compiled in dev promgram so using namespace std; sould be define so say this program is c++
sir how can write a program using bisection method of function x-cos
how i can write a program using bisection method of function x-cosx
Bisection Method in C# :
using System;
namespace Application1
{
class Program
{
public double c;
public double func(double x)
{
return x * x * x – 2 * x * x + 3;
}
public void bisection(double a, double b, double e)
{
Program func = new Program();
if (func.func(a) * func.func(b) >= 0)
{
Console.WriteLine(“Incorrect a and b”);
return;
}
c = a;
while ((b – a) >= e)
{
c = (a + b) / 2;
if (func.func(c) == 0.0)
{
Console.WriteLine(“Root = ” + c);
break;
}
else if (func.func(c) * func.func(a) < 0)
{
Console.WriteLine("Root = " + c);
b = c;
}
else
{
Console.WriteLine("Root = " + c);
a = c;
}
}
}
public static void Main(string[] args)
{
double a, b, e;
Console.WriteLine("Enter the desired accuracy:");
e = Convert.ToDouble(Console.ReadLine());
Console.WriteLine("Enter the lower limit:");
a = Convert.ToDouble(Console.ReadLine());
Console.WriteLine("Enter the upper limit:");
b = Convert.ToDouble(Console.ReadLine());
Program bisec = new Program();
bisec.bisection(a, b, e);
}
}
}
how can i write c++ program for bisection method using class and object..?????
Thanks. I really needed it. Thank you