function [coef_cijk] = cijk(bij, bik, bjk, qi, qj, qk)
%CIJK Summary of this function goes here
%   Detailed explanation goes here
coef_cijk = simplify(0.5*(diff(bij,qk) + diff(bik,qj)-diff(bjk,qi)))
end