инфа / проги инфа вар 1,2(родионова) / задача 1 2 сем

common c, dd
external f,f1
dimension aa(6), dd(5), v(4)
REAL xc,aa,eps,dd,c,cor,v
integer i
data aa/0.2,1.5,-5,3,16,7/
write(*,*) ' mas a'
16 format(1x,11(f5.2,1x))
write(*,16) (aa(i), i=1,6)
data eps/0.1/
!write(*,*) 'vvedite eps'
!read(*,*) eps
call difpol(aa,dd,6)
p=pol(dd,5,2)
write(*,*) 'proizvod', p
do i=1,4
c=i*p
xc=COR(F1,0.0,1.0,eps)
v(i)=xc
end do
write(*,*) 'koren yravn', v
end

subroutine difpol(a,b,m)
dimension A(m+1),b(m)
real a,b
integer i,m
do i=1,m-1
b(i)=A(i)*(m+1-i)
end do
end

real function pol(a,m,x)
Dimension A(m)
real a
integer i,m, x
p=a(1)
do i=2,m,1
p=p*x+a(i)
end do
pol=p
end

Real function COR(F,A,B,E)
Real F,A,B,E,X
Logical P
XL=A
XP=B
p=.false.
Do while (abs(XL-XP)>E .and..not.p)
XC=(XL+XP)/2
FXL=F(Xl)
FXC=F(Xc)
if (FXL*FXC>0)then
XL=Xc
else
if(FXL*FXC<0)then
XP=Xc
else
p=.true.
end if
end if
end do
COR=Xc
end

real function f(x)
common aa
dimension aa(4)
real x
f=aa(1)*x**5-aa(2)*x**4+aa(3)*x**3-aa(4)*x**2+aa(5)*x+aa(6)
end

real function f1(x)
common c
real x
f1=x**3+c*x-c
end