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