update bilinear test case operating point

def get_x0_bl_simple():

    # nonlinear operating point: 

    # for failure case make recycle really large, other streams very small

    x0_bl = {

        'm1': 0.2,

        'm2': 0.2,

        'm3': 0.5,

        'm4': 0.1,

        'm5': 0.4,

        'm6': 0.1,

        'm7': 0.3,

        'm8': 0.3,

        'x1': 0.01, 

        'x2': 0.01,

        'x3': 0.01,

        'x4': 0.01,

        'x5': 0.01,

        'x6': 0.01,

        'x7': 0.01,

        'x8': 0.01,

    }

    return x0_bl

    x0 = {}

    x0['m1'] = 60

    x0['m8'] = 100

    x0['m2'] = x0['m1']

    x0['m3'] = x0['m2'] + x0['m8']

    x0['m4'] = 20

    x0['m5'] = x0['m3'] - x0['m4'] 

    x0['m6'] = 40

    x0['m7'] = x0['m8']

    x0['x1'] = 80

    x0['x2'] = x0['x1']

    x0['x5'] = 40

    x0['x6'] = x0['x5']

    x0['x7'] = x0['x5']

    x0['x8'] = x0['x7']

    x0['x3'] = (x0['m2']*x0['x2']+x0['m8']*x0['x8'])/x0['m3']

    x0['x4'] = (x0['m3']*x0['x3']-x0['m5']*x0['x5'])/x0['m4']

    # print("f3=", x0['m3']*x0['x3'],"f5=", x0['m5']*x0['x5'])

    # for i, iconstr in constr_bl.items():

    #     print(i, "  ", iconstr, "  ",iconstr.subs(x0))

    # for key, val in x0_bl.items():

    #     x0_bl[key] = int(val)

    return x0