1function [v_CE, dA, F_m_out, F_CE, tau_m] = loopMuscles(m, qrel, l_CE, A, S)
 2%LOOPMUSCLES Summary of this function goes here
 3%   Detailed explanation goes here
 4
 5%#codegen
 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