MODELO NEURONAL DE FITZHUGH-NAGUMO

V' (t) = -V*(V - (Subscript[V, 1]))*(V - (Subscript[V, 2])) - W + 
        E;    V (0) = 0
\(W' \((t)\) = ϵ*\((V - 
            C*W)\); \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ W \((0)\) = 0\)

f[V_, W_] := -V*(V - (Subscript[V, 1]))*(V - (Subscript[V, 2])) - W + 
        e
g[V_, W_] := ϵ*(V - c*W)

(Subscript[V, 1]) := 2/10
(Subscript[V, 2]) := 1
e := 23/100
ϵ := 2/100
c := 0.5

f[V, W]
g[V, W]

Solve[{f[V, W], g[V, W]} == {0, 0}, {V, W}]

Sólo existe un punto de equilibrio en V≥0, W≥0, de valor  aproximado (0.110603 , 0.221206)

Simplify[f1[V_, W_] = D[f[V, W], V]]
f2[V_, W_] = D[f[V, W], W]
g1[V_, W_] = D[g[V, W], V]
g2[V_, W_] = D[g[V, W], W]

jaco=(({{f1[0.110603, 0.221206], f2[0.110603, 0.221206]}, {g1[0.110603, 0.221206], g2[0.110603, 0.221206]}}))

N[Eigenvalues[jaco]]

<<Graphics`PlotField`

PlotVectorField[{f[V,W], g[V,W]},{V, 0, 0.5}, {W,0, 0.7}]