BVProblem with solver RadauIIa5: Two incompatible array sizes to boundary condition function

General Usage
2

In case anyone is interested, here is the corrected version:

#  This example using BVProblem with solver RadauIIA5 is based on
# an attmpt to reproduce the analytic solution... 
#   y(t) = exp(-t^2) + 10.

using BoundaryValueDiffEqFIRK  

function f!(dy_dt, y, par, t)

#  It seems that t is provided as a vector.
#  (A vector might have been expected only for an "orthogonal collocation method".)

    dy_dt[1] = .-2.0.*t.*(y[1] .- 10.0)  

end

function initialGuess(t)
#  And here, t seems to be a vector
# The arguments of this function appear to be undocumented.

# We use the correct solution but without the constant as the initial guess 

    [exp.(-t.*t)] 

end

function bc!(resid, y, par, t)   #  Boundary conditions 
# Here again t seems to be a vector, which seems inevitable, 
#  as the documentation says we can constrain the entire tspan.
    
    resid[1] = abs.(y[1,1] .- 11.0)     # Gave quick result 

end

tspan = (0.0, 10.0)

bvp = BVProblem(f!, bc!, initialGuess , tspan, 0.0)

sol = solve(bvp, RadauIIa5(); dt = 0.1, abstol=1e-2, reltol=1e-2)