I&#39;m helping a colleague with some numerical work; one thing I need to do is to create a small procedure which produces the difference between an integral and its trapezoidal approximation; such as:<br><br><span style="font-family: courier new,monospace;">dt(f,a,b):=float(abs(integrate(f(s),s,a,b)-(f(a)+f(b))*(b-a)/2));</span><br style="font-family: courier new,monospace;">
<br>But this doesn&#39;t work; or at least, when I tried it earlier I got a &quot;Too many contexts&quot; error message.&nbsp; Is there an approved technique for writing a procedure which uses functions as inputs?<br><br>Thanks,<br>
Alasdair<br>