plot ib-vbe curves as function of six values of vce for a common emitter bjt in matlab

main.m
% Define parameters
Ib_values = [0.1, 0.2, 0.3, 0.4, 0.5, 0.6]; % Set of Ib values
Vce_values = 0:2:10; % Set of Vce values

% Constants
Vt = 26e-3; % Thermal voltage

figure; % Create a new figure for the plot

for i = 1:length(Ib_values)
    Ib = Ib_values(i);
    Ic = 100*Ib; % Assume beta = 100
    Ies = Ic/1000;
    
    % Calculate IB-VBE curve
    VBE = linspace(0, 0.8, 100);
    IB = Ies*(exp(VBE/Vt) - 1);
    
    % Calculate IC-IE curves for different VCE values
    for j = 1:length(Vce_values)
        VCE = Vce_values(j);
        IC = beta*(IB - Ies*(exp(VBE/Vt) - 1)) - (1+beta)*VCE;
        plot(IB, IC, 'DisplayName', sprintf('Vce = %dV', VCE));
        hold on;
    end
end

xlabel('IB (Base current)');
ylabel('IC (Collector current)');
title('Common Emitter BJT IB-VBE Curves at different VCE values');
legend;
grid on;
843 chars
33 lines