“…with c S = {1, 2, 3, 5, 6, 8, 9} and c D = {4, 7, 10}. The benefit of the proposed cost function, J exc , was highlighted by comparing the produced exciting motions with those obtained from the minimization of the condition number of the total BP regressor J cond = cond(W b ) which is commonly used in industrial robotics [11]. The two optimization processes were to find the joint angle values at N K = 4 via points, equally distributed over the whole trajectory of T F = 2s, that minimize the above mentioned criteria under the constraint that the B-spline interpolated joint trajectories respect the joint limitations q − j ≤ BS j (kT s, q j ) ≤ q + j , k = 1, ..., (n − 1).…”