f$:="p^2*cos(q)"
input a,b,c.d;
input e,l;
x1:=b-a;
y1:=d-c;
p:=a;
q:=c;
v1:=val(f$):
p:=b;
v1:=v1+val(f$);
q:=d;
v1:=v1+val(f$);
p:=a;
v1:=v1+val(f$);
v:=v1;
for i:=1 to l-1;
y2:=y1;
x2:=x1;
y1:=y1*0.5;
x1:=x1*0.5;
y:=c;
for j:=1 to 2^(i-1)+1;
q:=y;
x:=a+x1;
v2:=0;
for k:=1 to 2^(i-1);
p:=x;
v2:=v2+val(f$);
x:=x2;
next k;
if j=1 or j=2^(i-1)+1 then v2:=v2*0.5;
v:=v+v2*4;
y:=y+y2;
next j;
y:=c+y1;
for j:=1 to 2^(i-1);
q:=y;
x:=a;
v2:=0;
for k:=1 to 2^i+1;
p:=x;
if k=1 or k=2^I+1 then v2:=v2+2*val(f$);
goto l1;
v2:=v2+4*val(f$);
label l1;
x:=x+x1;
next k;
v:=v+v2;
y:=y+y2;
next j;
v3:=v*(b-a)*(d-c)/4^(i+1);
if abs(v3-v1)
next i;
i:=i-1;
l2;
print v3;
end;