function [ig] = information_gain(feat_class1, feat_class2,x)

% compute the information gain for features in feat_class1 and feat_class2
% given the decision boundary x in hessian normal form
% hint: augment the features with the constant one to compute the distance
% to the decision boundary


end