Numerical integration/Adaptive Simpson's method: Difference between revisions
Numerical integration/Adaptive Simpson's method (view source)
Revision as of 15:16, 17 May 2023
, 1 year ago→{{header|Pascal}}: Do not redundantly pass around f.
(Added Pascal before →{{header|Perl}}) |
(→{{header|Pascal}}: Do not redundantly pass around f.) |
||
Line 1,076:
depth : integer) : real;
procedure quad_asr_simpsons_ (
var
begin
m := (a + b) / 2;
Line 1,085 ⟶ 1,084:
end;
function quad_asr_ (
tol, whole, m, fm : real;
depth : integer) : real;
Line 1,094 ⟶ 1,092:
delta, tol_ : real;
begin
quad_asr_simpsons_ (
quad_asr_simpsons_ (
delta := left + right - whole;
tol_ := tol / 2;
Line 1,101 ⟶ 1,099:
quad_asr_ := left + right + (delta / 15)
else
quad_asr_ := (quad_asr_ (
left , lm, flm, depth - 1)
+ quad_asr_ (
right, rm, frm, depth - 1))
end;
Line 1,110 ⟶ 1,108:
fa := f(a);
fb := f(b);
quad_asr_simpsons_ (
quad_asr := quad_asr_ (
end;
|