function?y=gneiguzhangzi（tx）
global?mi?mb?mo?c1?c2?c3?c4?dita1?dita2?dita3?dita4;
epson=0.096;
n=16;
l=10;
rb=0.018;
ri=0.060;
Fd=4215;
Fmax=8430;
Fr=10;
wi=2.4.*pi;
wm=wi.*ri./（ri+rb）;
wb=（wi.*ri）./（2.*rb）;
dd=wi.*t+2*pi/5;
for?k=1:1:n;
????cc=wm.*t+2.*（k-1）.*pi./16;
????cc=mod（cc2*pi）;
????if?0????????Acita=1;
????elseif??2*pi/5????????Acita=0;?????
????else
????????Acita=1;
????end
end
F=Fmax.*（1-（1-cos（cc））./2/epson）.^1.5;

for?k=1:1:n;
if?2.*k.*pi./wi?????R=1;
else
????（2.*k.*pi+2.*2*pi/5）./wi????R=0;
end
end
T=2.*pi.*（ri+rb）./（n.*wi.*rb）;
for?k=0:1:n-1;
if?k.*T????S=1;
else
????k.*T+n.*T.*asin（l./2.*ri）./pi