Anonymous user
Numerical integration: Difference between revisions
→{{header|Ada}}
MaiconSoft (talk | contribs) No edit summary |
|||
Line 165:
Sum : Scalar := F(A) + F(B);
X : Scalar := 1.0;
begin
while X <= Scalar (N) - 1.0 loop
Sum := Sum + 2.0 * F (A + X * (B - A) / Scalar (N));
Line 175:
function Simpsons (A, B : Scalar; N : Positive) return Scalar is
H : constant Scalar := (B - A) / Scalar (N);
begin
for I in
else
Sum_E := Sum_E + F (A + H * Scalar (I));
end if;
end loop;
return (H /
end Simpsons;
end Integrate;</lang>
| |||