Update fractional_to_ordinary.jl

Author: ChrisRackauckas

test/fractional_to_ordinary.jl

@@ -21,8 +21,7 @@ sys = fractional_to_ordinary(eqs, x, α, 10^-7, 1)
 prob = ODEProblem(sys, [], tspan)
 sol = solve(prob, radau5(), saveat=timepoint, abstol = 1e-10, reltol = 1e-10)
 
-time = 0
-while(time <= 1)
+for time in 0:0.1:1
     @test isapprox(expect(time, α), sol(time, idxs=x), atol=1e-7)
     time += 0.1
 end
@@ -35,10 +34,8 @@ sys = fractional_to_ordinary(eqs, x, α, 10^-7, 1; matrix=true)
 prob = ODEProblem(sys, [], tspan)
 sol = solve(prob, radau5(), saveat=timepoint, abstol = 1e-10, reltol = 1e-10)
 
-time = 0
-while(time <= 1)
+for time in 0:0.1:1
     @test isapprox(expect(time, α), sol(time, idxs=x), atol=1e-7)
-    time += 0.1
 end
 
 α = 0.9
@@ -49,10 +46,8 @@ sys = fractional_to_ordinary(eqs, x, α, 10^-7, 1)
 prob = ODEProblem(sys, [], tspan)
 sol = solve(prob, radau5(), saveat=timepoint, abstol = 1e-10, reltol = 1e-10)
 
-time = 0
-while(time <= 1)
+for time in 0:0.1:1
     @test isapprox(expect(time, α), sol(time, idxs=x), atol=1e-7)
-    time += 0.1
 end
 
 # Testing for example 2 of Section 7
@@ -88,4 +83,4 @@ msys = linear_fractional_to_ordinary([3, 2.5, 2, 1, .5, 0], [1, 1, 1, 4, 1, 4],
 mprob = ODEProblem(sys, [], tspan)
 msol = solve(prob, radau5(), abstol = 1e-5, reltol = 1e-5)
 
-@test isapprox(expect(5000), msol(5000, idxs=x_0), atol=1e-5)
+@test isapprox(expect(5000), msol(5000, idxs=x_0), atol=1e-5)