function y=c9mbvp1(x)
% c9mbvp1 - MATLAB function to boundary value problem

% Copyright (c) Dingyu Xue, Northeastern University, China
% Last Modified 18 May, 2022
   W=get_param(gcs,'ModelWorkspace'); assignin(W,'a',x)
   try
      txy=sim('bp6_model.slx'); y0=txy.yout; 
      y=abs(y0(end)-exp(-2));
   catch, y=10; 
   end
end