function [bi,bieq] = fun_3_c_bivillk(x)
bieq=pi.*x(1).^2.*x(2)-1;
bi=[];