using Plots
k=(1:60)
## Start with linear
seq = 0.85.^k
plot(k, seq, yscale=:log10, linecolor=colorant"black", linewidth=2, label="Linear",
ylims=(10 .^ (-8.), 10 .^ (-0.)))
## yaxis
ylabel!("Error")
## Algebraic
plot!(k, 1 ./ k, linecolor=RGB(0,0.8,0), linewidth=2, label="Algebraic")
## Superlinear
seq = (2*k).^(-0.05*k)
plot!(k, seq, linecolor=colorant"blue", linewidth=2, label="Superlinear")
## Quadratic
k=1.:20.
seq = (0.99).^(2 .^ k)
plot!(k, seq, linecolor=colorant"magenta", linewidth=2, label="Quadratic")
## Cubic
seq = (0.99).^(3 .^ k)
plot!(k, seq, linecolor=colorant"yellow", linewidth=2, label="Cubic")
k=(1:60)
## Start with linear
seq = 0.85.^k
plot(k, seq, yscale=:log10, linecolor=colorant"black", linewidth=2, label="Linear",
ylims=(10 .^ (-8.), 10 .^ (-0.)))
## yaxis
ylabel!("Error")
## Algebraic
plot!(k, 1 ./ k.^2, linecolor=RGB(0,0.8,0), linewidth=2, label="Algebraic")
## Superlinear
seq = (2*k).^(-0.05*k)
plot!(k, seq, linecolor=colorant"blue", linewidth=2, label="Superlinear")
## Quadratic
k=1.:20.
seq = (0.99).^(2 .^ k)
plot!(k, seq, linecolor=colorant"magenta", linewidth=2, label="Quadratic")
## Cubic
seq = (0.99).^(3 .^ k)
plot!(k, seq, linecolor=colorant"yellow", linewidth=2, label="Cubic")
xlabel!("Iteration number")