function dunc_demo(n, p)
% demo the distance to uncontrollability codes
A = randn(n, n);
B = randn(n, p);
[upper,lower,z] = dunc_trisect(A, B, 1e-8)
[upper,lower,z,localmins] = dunc_fast(A, B, 1e-8)
