[VA] SRC #012d - Then and Now: Area

[J-F Garnier]

Update of my initial solution to include the "tiny" area at y= -4.05 .. -4.08 :

Program line changes or additions:
241 ! RHS<1 FOR Y=-4.05 !
242 F=FNROOT(-4,-4.05,FNR1(FVAR))
243 G=FNROOT(-4.05,-4.1,FNR1(FVAR))
244 I4=INTEGRAL(G,F,E*10000,FNF(IVAR)) @ DISP I4;IBOUND
250 S=I0+I2+I1+I3
260 T1=TIME-T1 @ DISP "AREA=";S+I4;"in";T1
265 DISP "(MAIN=";S;" TINY=";I4;")"
310 X2=ABS(Y+X)^(1/3) @ IF Y<-X THEN X2=-X2

New result (HP-71B):
RHS=1 for Ymax= 2.82740261413
LHS=RHS for X= 0 and Y= .831971149978
Ix IBOUND
5.42070674749E-4 5.42070677958E-13
1.23815989579 1.23815999491E-12
.023696731798 2.36967317996E-13
.814227669493 8.14227500886E-13
7.1976193049E-5 7.19937671652E-13
AREA= 2.07669834394 in 5594 s (about 1h33min)
(MAIN= 2.07662636775 TINY= 7.1976193049E-5 )

J-F