Matlab - SIMULACIÓN PÉNDULO FOUCAULT

function dx = focaultPendulum(t,x)

%define your constants

lat = 0*pi/180;

g = 9.81;

L = 5;

Omega = 7.27*10^-5;

dx = zeros(4,1);

dx(1) = x(3);

dx(2) = x(4);

dx(3) = 2*(Omega)*x(4)*sin(lat) - (g/L)*x(1);

dx(4) = -2*(Omega)*x(3)*sin(lat) - (g/L)*x(2);

end

[T,Y] = ode45(@focaultPendulum,[0 20],[0.7 0.5 0.6 0.3]);

plot(T,Y(:,1));

hold all

plot(T, Y(:,2))

plot(T, Y(:,3))

plot(T, Y(:,4))