clear
M = 1;
K = 1;
VarB = [0:0.1:5];
k = 1
N = length(VarB);
for I = 1:N,
    clc
    k = k+1
    B = VarB(I);
    sim('BASIC_MBK_System.mdl');
    n = length(performance);
    J(I) = performance(n);
    U(I) = RMSEf(n);
    E(I) = RMSE(n);
end
clf
plot(VarB,J,VarB,E,VarB,U),grid on, xlabel('B'), ylabel('J, E and U')