Z Skrypty dla studentów Ekonofizyki UPGOW

Oscylator bistabliny

Rozważmy oscylator nieliniowy:

\(m \ddot x = -\frac{dU(x)}{dx}-\gamma x ,\)

który możemy zaimplementować jako funkcję w matlabie:

function dx = ODEbistable(X,T)
    global gama;
    dx = zeros(2,1);
    dx(1) = X(2);
    dx(2) = -gama*X(2)-4*X(1).^3+8*X(1);
    return
end

Baseny przyciągania w nieliniowym oscylatorze.