function?[sJL_in]?=?carin（smm_inpvmax）
[ab]=size（s{1}）;
s_in=（rand（ab）s_new=（s{1}==0）.*s_in;
JL_in（1）=sum（sum（s_new））;%實(shí)際進(jìn)入的
JL_in（2）=sum（sum（s_in））;%本該進(jìn)入的
s{1}=s{1}+s_new;
s{2}=s{2}+（m==1）.*s_new.*randi（vmaxab）;
s{2}=s{2}-（m==2）.*s_new.*randi（vmaxab）;
s{3}=s{3}+（m==3）.*s_new.*randi（vmaxab）;
s{3}=s{3}-（m==4）.*s_new.*randi（vmaxab）;
end