(muH*Fzl/(2*FRl)))*muH*Fzl/(2*FRl); if FRr<=muH*Fzr/2 f1r=1; else f1r=Jr; end if FRl<=muH*Fzl/2 f1l=1; else f1l=Jl; end flon1=f1r*Cx1*s1; flon2=f1l*Cx2*s2; flat1=f1r*Cy1*a1; flat2=f1l*Cy2*aa2; % cc=10; % aa=