%% PVIs -- ejemplos
%% Logistica
clear 
clc
close all

a = 2; % parametro para f
f = @(t,y) a*y*(1-y); % y'(t)= f(t,y(t))
intervalo = [0, 10]; % en que intervalo queremos resolver el PVI
y0 = 0.1; % condicion inicial

[t,y] = ode45(f, intervalo, y0); % ode45 es una caja negra (por ahora)
plot(t,y, '-o','LineWidth',2)


%% Masa-resorte
clear 
clc
close all

m = 1; % masa
k = 1; % cte del resorte
mu = 0.5; % coef de rozamiento con el piso
A = [0 1; -k/m -mu/m];
f = @(t,y) A*y; % y'(t)= A*y(t)
intervalo = [0, 10]; 
y0 = [0.5; 0]; % posicion y velocidad iniciales (en ese orden)

[t,y] = ode45(f, intervalo, y0);
plot(t,y, '-o','LineWidth',2)