Вычислить по схеме Горнера значение многочлена с комплексными коэффициентами

program l3
    implicit none
    integer re_i , n , i
    real x , y
    real,allocatable :: a( : )
    real,allocatable :: b( : )
    write (*,*) 'Vvedite n'
    read(*,*) n
    allocate (a(n))
    allocate (b(n))
    call random_number(a)
    do i=1, n
        a(i)=100*a(i)-50
        a(i)=int(a(i))
    enddo
write (*,*) a
write(*,*) '----------------'
    call random_number(b)
    do i=1, n
        b(i)=100*b(i)-70
        b(i)=int(b(i))
    enddo
write(*,*) b
end
subroutine plus (x1,x2,y1,y2,xp,yp)
    real x1,x2,y1,y2,xp,yp
    xp=x1+x2
    yp=y1+y2
    end
     subroutine mult (x1,x2,y1,y2,xp,yp,i)
xp=x1*x2-y1*y2
yp=x2*y1+x1*y2
end

program l3
    implicit none
    integer n , i
    real a, b, x, y
    complex z, poly
    open (6, file = 'Data', status = 'old')
! В файле data иссходные данные расположены в следующем порядке:
! n
! x
! a(n), b(n)
! a(n-1), b(n-1)
! a(1), b(1)
! в файле как минимум, должны быть введены n, x и y
! по окончании исходных данных управление передается на метку 3
! Ниже использована экономичная схема Горнера с обраным ходом    
    z = 0
    read (2,*) n
    read (2,*) x, y
    z = cmplx(y,y)
    read (2,*, end=3) a, b
    poly = cmplx(a,b)
    do i = n-1, 1
        read (2,*) a, b
        poly = z*poly + cmplx(a,b)
    end do
3   write (*,*) 'Sum of polynom=', poly
end program l3