HP42s first major program (Double Integral) Best way to approach?

[Werner]

Easy fix

Code:

00 { 197-Byte Prgm }
01▸LBL "BORE"
02 MVAR "a"
03 MVAR "b"
04 MVAR "c"
05▸LBL 10
06 VARMENU "BORE"
07 RTN
08 RCL "a"
09 XEQ 03
10 LSTO "t"
11 RCL "b"
12 XEQ 03
13 RCL- "t"
14 GTO 10
15▸LBL 03
16 RCL "c"
17▸LBL "HV"
18 X=Y?
19 GTO 00
20 X>Y?
21 X<>Y
22 LSTO "d"
23 X<>Y
24 LSTO "D"
25 ÷
26 XEQ "EK"
27 RCL- ST Y
28 X<>Y
29 LASTX
30 +
31 RCL "d"
32 X↑2
33 ×
34 X<>Y
35 RCL "D"
36 X↑2
37 ×
38 +
39 RCL× "D"
40 3
41 ÷
42 RTN
43▸LBL 00
44 3
45 Y↑X
46 LASTX
47 ÷
48 STO+ ST X
49 RTN
50▸LBL "EK"
51 X↑2
52 SIGN
53 LASTX
54 -
55 SQRT
56 SIGN
57 LASTX
58 STO+ ST Y
59 SQRT
60 STO+ ST X
61 XEQ "AGM2"
62 PI
63 X<>Y
64 ÷
65 STO× ST Y
66 X<>Y
67 4
68 ÷
69 RTN
70▸LBL "AGM2"
71 LSTO "B"
72 SIGN
73 +/-
74 LSTO "T"
75 R↓
76 LSTO "A"
77 X↑2
78▸LBL 02
79 RCL "A"
80 RCL× "B"
81 SQRT
82 X<> "B"
83 RCL- "A"
84 2
85 STO× "T"
86 ÷
87 STO+ "A"
88 X↑2
89 RCL× "T"
90 RCL+ ST Y
91 X≠Y?
92 GTO 02
93 RCL "B"
94 END

Werner