不能使用numpy库，给我实现用四阶隐形Adams法求解区间[a,b]的常微分初值问题（柯西问题），并且给我代码，我需要输入f epsilon a y(a) b

好的，以下是使用四阶隐式Adams法求解常微分初值问题的代码：

def implicit_adams(f, eps, a, ya, b):
    """
    :param f: 函数f
    :param eps: 精度
    :param a: 区间左端点
    :param ya: 初值y(a)
    :param b: 区间右端点
    """
    h = eps**(1/4)  # 步长
    t = [a, a+h]  # 时间节点
    w = [ya, ya+h*f(a, ya)]  # 解节点
    while t[-1] < b:
        tn = t[-1] + h
        wn = w[-1] + h*(1/720)*(251*f(t[-1]+3*h/2, w[-1]+3/2*k3) + \
                                  646*f(t[-1]+h/2, w[-1]+1/2*k2) - \
                                  264*f(t[-1]+h, w[-1]+k1) + \
                                  106*f(tn, w[-1]+h/720*(251*k4-646*k3+264*k2)))
        # 计算误差
        err = eps*abs(wn/(b-a))
        # 判断是否需要调整步长
        if err <= eps:
            t.append(tn)
            w.append(wn)
        h = 0.9*h*(eps/err)**(1/4)
    return t, w

其中，k1、k2、k3、k4是隐式Adams法中的求解系数，计算方式如下：

# 计算k1、k2、k3、k4
k1 = f(t[-1], w[-1])
k2 = f(t[-1]+h/2, w[-1]+h/2*k1)
k3 = f(t[-1]+h/2, w[-1]+h/2*k2)
k4 = f(t[-1]+h, w[-1]+h*k3)

你可以按照以下方式调用函数：

t, w = implicit_adams(f, eps, a, ya, b)

其中，f为函数名，eps为精度，a为区间左端点，ya为初值y(a)，b为区间右端点。函数返回时间节点t和解节点w，可以通过它们来得到解函数。

需要注意的是，这里使用的是四阶隐式Adams法，因此初始点需要通过其他方法（如龙格-库塔法）得到。另外，该方法的实现需要依赖于函数f的实现，因此需要确保函数f在区间[a,b]内是连续可导的。