clear all
clc
close all
m = 0.2533; %masa
k = 10; %rigidez
c = 0.1592; %amortiguamiento
u(1) = 0; %desplazamiento inicial
u1(1) = 0; %velocidad inicial
u2(1) = 0;
p_0 = 0; %fuerza inicial
dt = 0.1; %intervalo de tiempo
t = 0:dt:1;
p = [0 5 8.6603 10 8.6603 5 0 0 0 0 0];
a1 = (6*m/(dt^2)) + (3*c/dt)
a2 = (6*m/dt) + (2*c)
a3 = (2*m) + (dt*c/2)
k1 = k + a1
for i = 1:1:10
p1(i+1) = (p(i+1) + (a1*u(i)) + (a2*u1(i)) +(a3*u2(i))); %2.1
u(i+1) = (p1(i+1))/k1; %2.2
u1(i+1) = (3*(u(i+1)-u(i))/dt) - (2*u1(i)) - (dt*u2(i)/2); %2.3
u2(i+1) = (6*(u(i+1)-u(i))/(dt^2)) - (2*u2(i)) - (6*u1(i)/dt); %2.4
end
