Skip to main content

fixed point iteration method using C++


F(x)=X2-6X+8
G(X)=X=(X2+8)/6;
We have to initiate a value then put it on G(X) and get value of G(X) and then put this on F(X).it continues until it F(X)=0 or we get tolerable fraction into F(X) to output or finishing maximum iteration step. 
///Bismillahir Rahmanir Rahim
///Author: Tanvir Ahmmad
///CSE,Islamic University,Bangladesh
#include<iostream>
#include<math.h>
using namespace std;
double func(double num)
{
    return ((num*num)-(6*num)+8);
}
double value(double num)
{
    return ((8+(num*num))/6);
}
int main()
{
    double guess,tol,val,fund,ans;
    int it_time,i;
    cout<<"Enter initial guess: ";
    cin>>guess;
    cout<<"Enter tolerable error: ";
    cin>>tol;
    cout<<"Enter maximum iteration: ";
    cin>>it_time;
    for(i=1; i<=it_time; i++)
    {
        val=value(guess);
        if(abs(func(val))>=tol)
        {
            ans=val;
            fund=func(val);
            cout<<"Iteration-"<<i<<":     x =   "<<val<<" and f(x) =   "<<fund<<endl;
            guess=val;
        }
        else break;
    }
    if(i==(it_time+1)) cout<<"Using "<<it_time<<" step root is : "<<ans<<endl;
    else cout<<"Root is : "<<ans<<endl;
    return 0;
}


Comments

Popular posts from this blog

Codeforces round 1676(A. Lucky?)

Just count the  first three  &  last three  number ///Bismillahir Rahmanir Rahim ///Author: Tanvir Ahmmad ///CSE,Islamic University,Bangladesh #include< iostream > #include< cstdio > #include< algorithm > #include< string > #include< cstring > #include< sstream > #include< cmath > #include< cstring > #include< vector > #include< queue > #include< map > #include< set > #include< stack > #include< vector > #include< iterator > #include   < functional >   ///sort(arr,arr+n,greater<int>()) for decrement of array /*every external angle sum=360 degree   angle find using polygon hand(n) ((n-2)*180)/n*/ ///Floor[Log(b)  N] + 1 = the number of digits when any number is repres...

Codeforces round 1661(B. Getting Zero)

using   second operation  the answer is maximum is  15  because  2 15 =32768 so we need to  pre-calculate all the illegible input 0 to 32768   using  second operation Then  just increment 1(using first operation)  from input check is it  minimum is not ! ///Bismillahir Rahmanir Rahim ///Author: Tanvir Ahmmad ///CSE,Islamic University,Bangladesh #include< iostream > #include< cstdio > #include< algorithm > #include< string > #include< cstring > #include< sstream > #include< cmath > #include< cstring > #include< vector > #include< queue > #include< map > #include< set > #include< stack > #include< vector > #include< iterator > #include   < functional >   ///sort(arr,arr+n,greater<int...

Codeforces round 1676(E. Eating Queries)

Just use the lower bound & check is it greater than the sum or not Lower bound:https://www.geeksforgeeks.org/upper_bound-and-lower_bound- for -vector- in -cpp-stl/ ///Bismillahir Rahmanir Rahim ///Author: Tanvir Ahmmad ///CSE,Islamic University,Bangladesh #include< iostream > #include< cstdio > #include< algorithm > #include< string > #include< cstring > #include< sstream > #include< cmath > #include< cstring > #include< vector > #include< queue > #include< map > #include< set > #include< stack > #include< vector > #include< iterator > #include   < functional >   ///sort(arr,arr+n,greater<int>()) for decrement of array /*every external angle sum=360 degree   angle find using polygon hand(n) ((n-2)*180...