1function [v_CE, dA, F_m_out, F_CE, tau_m] = loopMuscles(m, qrel, l_CE, A, S)
2
3
4
5
6nMuscles = length(A)/2;
7v_CE = zeros(2*nMuscles,1);
8dA = zeros(2*nMuscles,1);
9F_m_out = zeros(2*nMuscles,1);
10F_CE = zeros(2*nMuscles,1);
11tau_m = zeros(6,1);
12
13refi = [repmat([1 2 3],7,1);repmat([4 5 6],7,1)];
14qmus = [repmat(qrel(1,:),7,1);repmat(qrel(2,:),7,1)];
15mCount = [1:7 1:7];
16
17for iMuscle = 1:nMuscles*2
18 [v_CE(iMuscle), dA(iMuscle), tau, F_m_out(iMuscle), F_CE(iMuscle)] = ...
19 muscle(...
20 m(mCount(iMuscle)), ...
21 qmus(iMuscle,:)', ...
22 l_CE(iMuscle), A(iMuscle), S(iMuscle));
23 tau_m(refi(iMuscle,:)) = tau_m(refi(iMuscle,:))+tau;
24end
25end
26