rectangle_rule was erron incl the final endpoint. Bump to 0.2.1

Author: fusion809

Project.toml

@@ -1,7 +1,7 @@
 name = "FunctionIntegrator"
 uuid = "7685536e-2581-4f83-bef1-2ba363c9cb91"
 authors = ["Brenton Horne <[email protected]>"]
-version = "0.2.0"
+version = "0.2.1"
 
 [deps]
 FastGaussQuadrature = "442a2c76-b920-505d-bb47-c5924d526838"

src/Rectangle.jl

@@ -5,18 +5,16 @@ numerically approximates:
 
 ``\\displaystyle\\int_a^b f(x) dx``
 
-using the [rectangle rule](https://en.wikipedia.org/wiki/Riemann_sum#Left_Riemann_sum) with ``N+1`` steps (including the endpoints). Whilst the type mentioned in the function definition for ``N`` is just 'Number' (as opposed to 'Integer'), this is just so that scientific notation can be used to define it (as scientific notation gives the type 'Float64'), an error message will be printed if it is not a positive integer. 
+using the [rectangle rule](https://en.wikipedia.org/wiki/Riemann_sum#Left_Riemann_sum) with ``N`` gridpoint values. Whilst the type mentioned in the function definition for ``N`` is just 'Number' (as opposed to 'Integer'), this is just so that scientific notation can be used to define it (as scientific notation gives the type 'Float64'); an error message will be printed if it is not a positive integer. 
 """
 function rectangle_rule(f::Function, N::Number, a::Number, b::Number)
     N = convert(Int64, N);
     h = (b-a)/N;
     y = 0;
     x = a;
-    for i=1:N+1
+    for i=1:N
         y += h*f(x)
-        if i < N+1
-            x = x + h;
-        end
+        x = x + h;
     end
     return y
 end

test/logoverx.jl

@@ -17,7 +17,7 @@ printstyled("Integrating log(x)/x from 1 to e and comparing the result to the an
     printstyled("Running: radau_quadrature\n"; color = :magenta)
     @time @test radau_quadrature(x -> log.(x).*x.^(-1), 8, 1, exp(1)) ≈ 0.5
     printstyled("Running: rectangle_rule\n"; color = :magenta)
-    @time @test rectangle_rule(x -> log.(x).*x.^(-1), 4.38e7, 1, exp(1)) ≈ 0.5
+    @time @test rectangle_rule(x -> log.(x).*x.^(-1), 4.3800001e7, 1, exp(1)) ≈ 0.5
     printstyled("Running: simpsons_rule\n"; color = :magenta)
     @time @test simpsons_rule(x -> log.(x).*x.^(-1), 100, 1, exp(1)) ≈ 0.5
     printstyled("Running: trapezoidal_rule\n"; color = :magenta)