Update Inverse eigenvalue problem.jl

MikaelSlevinsky · MikaelSlevinsky · commit 455e6a6078db · 2019-10-23T14:07:01.000-05:00
diff --git a/Extras/Inverse eigenvalue problem.jl b/Extras/Inverse eigenvalue problem.jl
@@ -12,14 +12,14 @@ import ApproxFun: bandwidth
 # To answer this question, we use a dual number to observe the sensitivity in λ₀
 # with respect to perturbations in b. The `realpart` of the objective function
 #
-#    f(λ) = λ - 1,
+#    f(λ(b)) = λ(b) - 1,
 #
-# extracts the difference, whereas the `dualpart` evaluates to f'(λ).
+# extracts the difference, whereas the `dualpart` evaluates to ∂f/∂b (λ(b)).
 #
 
 n = 100
-λ = dual(0.0,0.0)
-b = dual(eps(), 1.0)
+λ = dual(0.0, 0.0)
+b = dual(0.0, 1.0)
 v = zeros(Dual{Float64}, n)
 v[1] = 1
 v .+= 0.003.*randn(Float64, n)

Original file line number	Diff line number	Diff line change
`@@ -12,14 +12,14 @@ import ApproxFun: bandwidth`
`12`	`12`	`# To answer this question, we use a dual number to observe the sensitivity in λ₀`
`13`	`13`	# with respect to perturbations in b. The `realpart` of the objective function
`14`	`14`	`#`
`15`		`-# f(λ) = λ - 1,`
	`15`	`+# f(λ(b)) = λ(b) - 1,`
`16`	`16`	`#`
`17`		-# extracts the difference, whereas the `dualpart` evaluates to f'(λ).
	`17`	+# extracts the difference, whereas the `dualpart` evaluates to ∂f/∂b (λ(b)).
`18`	`18`	`#`
`19`	`19`
`20`	`20`	`n = 100`
`21`		`-λ = dual(0.0,0.0)`
`22`		`-b = dual(eps(), 1.0)`
	`21`	`+λ = dual(0.0, 0.0)`
	`22`	`+b = dual(0.0, 1.0)`
`23`	`23`	`v = zeros(Dual{Float64}, n)`
`24`	`24`	`v[1] = 1`
`25`	`25`	`v .+= 0.003.*randn(Float64, n)`