Use the analytic solution Y(x) to compute y1 in problem6_method.

Fix a bug (bad fraction) in problem6_method1.
Remove the partition_delta function; it's stupid in light of [a : delta : b].

diff --git a/partition.m b/partition.m
index 6568be38b39ba0580abef16b63d189ff4121e2fe..a08423671716d741302f67944344e3fbf76236db 100644 (file)
--- a/partition.m
+++ b/partition.m
@@ -24,5 +24,5 @@ function [p,delta] = partition(integerN, a, b)
   ## towards `b`.
   delta = (b - a)/integerN;
 
-  p = partition_delta(a,b,delta);
+  p = [a : delta : b];
 end

diff --git a/partition_delta.m b/partition_delta.m
deleted file mode 100644 (file)
index febb4dc..0000000
--- a/partition_delta.m
+++ /dev/null
@@ -1,19 +0,0 @@
-function [p] = partition_delta(a, b, delta)
-  ## Partition the interval [a,b] into subintervals of length delta. We
-  ## do not require that a<b.
-  ##
-  ## INPUTS:
-  ##
-  ##   * ``a`` - The "left" endpoint of the interval to partition.
-  ##
-  ##   * ``b`` - The "right" endpoint of the interval to partition.
-  ##
-  ##   * ``delta`` - The length of the subintervals.
-  ##
-  ## OUTPUTS:
-  ##
-  ##   * ``p`` - The resulting partition, as a vector of length (b-a)/delta.
-  ##
-
-  p = [a : delta : b];
-end