equations = :[
        iA = 0.0

        Lff = yLff(phi) # function table with linear interpolation
        Laa = yLaa(phi)
        Laf = yLaf(phi)
        Lfa = yLfa(phi)

        ΨF = Lff*iF + Laf*iA
        ΨA = Laa*iA + Lfa*iF

        dΨAdt = der(ΨA)
        dΨFdt = der(ΨF)

        uF =  dΨFdt + iF*Rf

        # mechanisch
        der(phi) = ωmech
        tau = (ydLaa(phi)*iA^2 + ydLff(phi)*iF^2)*0.5 + ydLfa(phi)*iF*iA

        Inertia*der(ωmech) = 0.0 - Load(time)
    ]