def METHOD_NAME(linear_solver, nfe_x=100, nfe_t=200):
    m = build_burgers_model(nfe_x=nfe_x, nfe_t=nfe_t)
    nlp = PyomoNLP(m)
    sqp(nlp, linear_solver)
    load_solution(m, nlp)
    return pe.value(m.obj)