bisection method to find roots

[ssushchy]

I've worked out a method to solve the roots and too print the root. but I want to show the values at each iteration step. every time i seem to do something i mess the rest of my code up. how would i go about doing that? this is what i have currently.

#include<iostream>
#include <stdio.h>
#include <math.h>
using namespace std;

double f(double x);
double bisection(double xa, double xb, double eps);

int main()
{
double x0 = 0.0, xn = 1.0, eps = 0.000001;
int n = 10000;

cout<< "Root value: "<< bisection(x0, xn, eps)<<endl;

return 0;
}

double f(double x)
{
double y;

y = x-(.2*sin(x))-0.5;
return y;
}

double bisection(double xa, double xb, double eps)
{
double xm, fa, fm;

fa = f(xa);
        while (xb - xa > eps)
{
        xm = (xa + xb) / 2;
        fm = f(xm);
        if (fm == 0.0)
break;
        else if (fa*fm > 0)
{
        xa = xm;
        fa = fm;
}
        else {
        xb = xm;
}
}
        return xm;
}