labels = 'ABCD';
for i = 1:4
    subplot(220+i);
    xlabel('Bias voltage (V)','fontsize',14);
    ylabel('C* (nF/cm^2)','fontsize',14);
    title(labels(i),'fontsize',14);
end