Caida fricción sobre resorte amortiguado

 function total2

clc
 % solucion caida libre
  time1=0:1/10:5;
condini=[100 0];
x=0*time1;
z=0*time1;  %Resuelvo las ecuaciones diff
[t,y]=ode45(@cfr,time1,condini); %solucion masa resorte
  time2=5:1/10:35;
alturafinal=y(length(y(:,1)),1);
velocidadfinal=y(length(y(:,2)),2);
condicionesiniciales=[alturafinal velocidadfinal];
x2=0*time2;
z2=0*time2;  % desarrolla las ecuaciones diff rk 4 orden
[t2,y2]=ode45(@fk,time2,condicionesiniciales); % Unir programa a simular
  tt=[time1 time2];
yt=[y;y2];
xt=[x x2];
zt=[z z2] % graficar resultado
  for i=1:length(tt)
pause(1/100)
subplot(2,2,1)
plot(tt,yt(:,2))subplot(2,2,2)
plot(tt,yt(:,1))subplot(2,2,3)
plot3(xt(i),zt(i),yt(i,1),'o r','MarkerFaceColor','m')
grid on
axis([-2 2 -2 2 min(yt(:,1)) max(yt(:,1))]) end
end  function rk=cfr(t,y)
% Ecuaciones diferenciales caida libre
  m=2;
b=0.01;
g=-9.81;
rk=[y(2);g-(b/m)*y(2)];
 end  function rk=fk(t,y)
m=2;
k=10;
c=1;
g=9.81;
rk=[y(2);g-(k/m)*y(1)-(c/m)*y(2)];
 end