texto: metodos numericos para ecuaciones diferenciales ordinarias

La ecuación de movimiento de la masa es
que se puede escribir como
F2 = nx ′
N = mg
mx ′′ = F − F1 − F2
mx ′′ = F − kx − nx ′
mx ′′ + nx ′ + kx = µN[1 − b(v0 − x ′ ) + c(v0 − x ′ ) 3 ]
Esta ecuación diferencial se puede llevar a un sistema de ecuaciones diferenciales haciendo
x1 = x
x2 = x ′
x ′ 1 = x2 (3.13)
x ′ 2 = − k
m x1 − n
m x2 + µN[1 − b(v0 − x2) + c(v0 − x2) 3 ]
m
(3.14)
Los datos del problema son b = 0.3, c = 0.1, m = 1kg,g = 9.8m/s 2 , u = 0.6,
k = 1600N/m, n = 0.1kg/s
Programa banda.m
function y=banda(r)
%Programa que determina el desplazamiento de una banda transportadora
clc;global v0;v0=r;
format long
a=0;
b=8;
m=2; %numero de ecuaciones
N=1000; %numero de subintervalos
c(1)=0;%condiciones iniciales
c(2)=0;
h=(b-a)/N;
t=a;
for j=1:m,
30