Function File: flow = corrector (flow, tol, opts, init)

applies a newton/levenberg-marquardt corrector to a flow state, in order to reach local lift/circulation balance. tol specifies tolerance. Reads the options minit, maxit, use_fsolve and silent from opts. Returns empty matrix if not successful.