Quadgk integration

[tbeason]

you need

function unknotpotential(a,b,c)
ans = quadgk(t->unknotintegralfunction(a,b,c,t),0,2*pi,order=quadorder)
end

(untested)

Also it looks like you are missing a definition for quadorder