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