Вычислить с заданной точностью значение заданной комплексной функции

#include <iostream>
using namespace std;
class Complex
{
public: double a, b;
public: Complex(){}
public: Complex(double a_, double b_)
{
    a = a_;
    b = b_;
}
        Complex operator*(Complex z)
        {
            return Complex(this->a*z.a - this->b*z.b, this->a*z.b + this->b*z.a);
        }
        Complex operator+(Complex z)
        {
            return Complex(this->a+z.a, this->b+z.b);
        }
        Complex operator/(double z)
        {
            return Complex(this->a/z, this->b /z);
        }
};
int main()
{
    const int eps = 0.0001;
    int i = 1;
    Complex z(2.0, 4.0),temp,sh(0.0,0.0);
    temp = z;
    sh = temp;
    do
    {
        temp = temp*z*z;
        temp = temp/((i + 1.0)*(i + 2.0));
        i += 2;
        sh = sh + temp;
    } while (temp.a*temp.a+temp.b*temp.b > eps*eps);
    cout << sh.a << " " << sh.b;
    system("pause");
    return 0;
}