%7DOF pulley system
%input parameters
mass = .2663; %kg, mass of pulley on tensioner spring
J = [0, 41.5, 13.1, 2.63, 42.1, 17.6, 2.08]/10000;  %kg-m^2, rotational inertias
r = [81.25, 64.5, 70.60, -41.15, 30, 67.5, -38.1]/1000; %m, radii of pulleys
x = [0, 261.5, 252, 90.3, 86, 0, 128]/1000; %m, x location of pulleys
y = [0, 60, 234, 251.1, 354, 167.5, 126]/1000; %m, y location of pulleys

EA = 65000; %N 
kT = 9340; %N/m, stiffness of tensioner spring, assumed to be linear

%frequency parameters for graphing
Ldiv = 10; %increment for omega in rpm
Lom = 3000; %max omega in rpm

%inital response plot
plotOn = 1; %shows plot
[wnrpm,V,Amat,list] = MDOFpulley(J,r,x,y,mass,EA,kT,Ldiv,Lom,plotOn);