I am trying to find a root of an equation using Interval Newton method in Julia. I am using “IntervalArithmetic.jl” package in “Pluto.jl” notebook. The code seems logically correct (checked a few times) but it is giving an error
syntax: "for" at /Users/dycobs/.julia/pluto_notebooks/Mini revelation.jl#==#4e903eeb-0172-4cd3-b5e9-21d235f99d28:13 expected "end", got "NXk"
The code is here
begin using IntervalArithmetic, Plots, ForwardDiff, PlutoUI, FiniteDiff f(x) = a - b * x l = 10 m = 10 a1 = vec(mince(100..200, l)) b1 = vec(mince(1..2, m)) d(f, x) = ForwardDiff.derivative(f, x) local Xnew = 0..300 nn = 10 NXk = (0..0) for i in 1:nn for j in 1:l a = a1[j] for k in 1:m b = b1[k] NXk_new = mid(Xnew) - (f(mid(Xnew)) / d(f, Xnew)) NXk = NXk ∪ NXk_new end end Xkplus NXk ∩ Xnew Xnew = Xkplus @show(interval(i), Xnew) end end
Could you please explain the possible reason behind the error?