нужна помощь с программой, не хочет работать
+нужно перевести программу на паскаль
могу скинуть условие если нужно
Код:
External F
integer c
real s,F,A,B,fint,E
common /z/c
E=0.001
B=1.3
do c=1,5,1
A=cor(F,0,1,E)
!a=c*0.2
S=fint(F,A,B,E)
write(*,10)S
10 format(x,'площадь равна=',F6.2)
pause
enddo
END
Real function F(x)
real x
common /z/c
F=x**3+c*x-c
end
real function fint(f,a,b,eps)
external f
real x,a,b,eps,sp,ss,h,d,f
integer n
n=100
sp=0.
11 ss=0.
h=(b-a)/n
x=a
do x=a,b,h
ss=ss+f(x)*h
end do
d=abs(sp-ss)
sp=ss
n=2*n
if (d>=eps) goto 11
fint=ss
end
real function COR(f,a,b,e)
external F
real x,e,f,xl,xp,u,v,d
integer a,b
logical P
xl=a
xp=b
U=f(xl)
v=f(xp)
p=.false.
do while(abs(xl-xp)>e.and..not.p)
x=xp-v*(xl-xp)/(u-v)
d=f(x)
if((u*d)>0)then
xl=x
u=d
else
if((u*d)<0)then
xp=x
v=d
else
p=.true.
endif
endif
enddo
cor=x
return
end