%函數(shù)f1存放于f1.m
function?x=f1（ijpr）
global?x;
x=char（x）;
if（j<=i）
???return;
else
????q=0;
for?t=i:j?%對于區(qū)間[ij]自上而下求累加概率值
????q=p（t）+q;y（t）=q;?
end
for?t=i:j%把所有自上而下的累加概率值與該區(qū)間總概率值減該累加概率值之差取絕對值存在一數(shù)組
v（t）=abs（y（t）-（q-y（t）））;?
end
for?t=i:j
????if（v（t）==min（v））?%求該