function [curvature] = Uk(k, n, coords)
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
    curvature = 0;
    for i = 1:n
       if i == k
       else
           if i == 1
               curvature = curvature + 15*Ujk(i, k, coords);
           else
               curvature = curvature + Ujk(i, k, coords);
           end
       end
    end
end