function dx = morrisLecar(t,x,params)
if length(params) == 0
    params = [-0.01, 0.15, 0.1, 0.05, -0.5, -0.7, 1, 0.5, 2, 1.2 0.005];
end
V1 = params(1);
V2 = params(2);
V3 = params(3);
V4 = params(4);
El = params(5);
Ek = params(6);
ECa = params(7);
gl = params(8);
gk = params(9);
gCa = params(10);
mu = params (11);

V = x(1);
w = x(2);
u = x(3);

dx = zeros(3,1);
dx(1) = gerjesztes(t, V, u) - gl*(V-El)-gk*w*(V-Ek)-gCa*mInf(V,V1,V2)*(V-ECa);
dx(2) = lambda(V,V3,V4)*(wInf(V,V3,V4)-w);
dx(3) = slow(t,u,V,mu);