@@ -150,13 +150,13 @@ The method return an array upar with the best fit parameters with their errors.
@. model(x,p) = p[1] + p[2] * exp(-(p[3]-p[1])*x)
fit_routine(model, ydata, param=3)
"""
function fit_routine(model::Function,ydata::Array{uwreal},param::Int64=3)
function fit_routine(model::Function,ydata::Array{uwreal},param::Int64=3;wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}},Nothing}=nothing)