数値計算プログラム

２関数積分も基本的な考えは同じです。
入れ子構造になります。
１関数と同じ台形則が、一番平凡で普通です。

f$:="(p^2*cos(q))" ２変数関数
input a,b
input c,d
input e,l
x1:=(b-a)*2
y1:=(d-c)*2
i:=1
repeat
x1:=x1*0.5
y1:=y1*0.5
y:=c
v:=0
i1:=2^(i-1)
j:=1
repeat
p:=a
q:=y
v2:=val(f$)
p:=b
v2:=v2+val(f$)
if i1:=1 then label1
x:=a
k:=1
repeat
c:=x+x1
p:=x
v2:=v2+2*val(f$)
k:=k+1
until k>=i1-1
label1
if j=1 or j=i1+1 then v2:=v2*0.5
v:=v+2*v2
y:=y+y1
j:=j+1
until j>=i1+1
v:=v*(b-a)*(d-c)/4^i
if i=1 then label2
if abs(v-v1) label2
v1:=v
i:=i+1
until i>=l
i:=i-1
label3
end