function [z]=brok(k,om,R)

z=4*(om-k).^2-R^2*((2*k-1).^2-exp(-4*k));


