def modified_euler(f, x0, y0, h, steps):

    for _ in range(steps):

        y_pred = y0 + h * f(x0, y0)

        y0 += 0.5 * h * (f(x0, y0) + f(x0 + h/2, y+y_pred/2))

        x0 += h

    return y0

def modified_euler(f, x0, y0, h, steps):

    for _ in range(steps):

        y_pred = y0 + h * f(x0, y0)

        y0 += 0.5 * h * (f(x0, y0) + f(x0 + h, y_pred))

        x0 += h

    return y0

def modified_euler(f, x0, y0, h, steps):

    for _ in range(steps):

        y_pred = y0 + h * f(x0, y0)

        y0 += h * f(x0 + h, y_pred)

        x0 += h

    return y0

def modified_euler(f, x0, y0, h, steps):

    for _ in range(steps):

        y_pred = y0 + h * f(x0, y0)

        y0 += 0.5 * h * (f(x0, y0) + f(x0 + h/2, y_pred))

        x0 += h

    return y0

def modified_euler(f, x0, y0, h, steps):

    for _ in range(steps):

        y_pred = y0 + h * f(x0, y0)

        y0 += 0.5 * h * (f(x0, y0) + f(x0 + h/2, y+y_pred))

        x0 += h

    return y0