complex<double> get_r(complex<double> *a, complex<double> *b, complex<double> z, int l){
    b[l] = 0;
    int k;
    for (k = l - 1; k >= 0; k--){
        b[k] = a[k + 1] + z * b[k + 1];
    }

    return a[0] + z * b[0];
}

int main(){

	complex<double> *a = new complex<double>[5];
	complex<double> *b = new complex<double>[5];
	complex<double> *c = new complex<double>[4];

    complex<double> *z = new complex<double>[IT_MAX + 1];
    complex<double> *R = new complex<double>[IT_MAX + 1];
    complex<double> *Rprim = new complex<double>[IT_MAX + 1];

	int l, j, k, n=4;
	double a_min = pow(10,-10);

	a[0] = 16. + 8.I;
	a[1] = -20. + 14I;
	a[2] = 4. - 8I;
	a[3] = -4. + 1I;
	a[4] = 1 + 0I;

    complex<double> z1 = 0 + 1I;
    complex<double> z2 = 1 + 1I;
    complex<double> z3 = -1 -3I;
    complex<double> z4 = 4 + 0I;

    for(int i = 0; i < 5; i++)
        cout << real(a[i]) << " " << imag(a[i]) << endl;

    complex<double> z0 = 0 + 0I;
    R[0] = get_r(a,b,z0, n);
    cout << "\nRj = " << R[0] << endl;
    cout << "\ntablica b\n";
    for(int i = 0; i < 5; i++)
        cout << real(b[i]) << " " << imag(b[i]) << endl;


    Rprim[0] = get_r(a,c, z0, n );
    cout << "\nRj' = " << Rprim[0];
    cout << "\ntablica c\n";
    for(int i = 0; i < 4; i++)
        cout << real(c[i]) << " " << imag(c[i]) << endl;

	return 0;
}