The documentation of xlabel explains with an example, how a multi-line output works. Did you read it? Then the only problem is the font size, which is explained there also:
plot((1:10).^2)
xlabel({'\fontsize{15}Population','\fontsize{20}(in thousands)'})
