hfst函数模块: function out = hfst(u1,u2,r,h) d=r*h; d0=h*d; y=u1+h*u2; a0=sqrt(d*d+8*r*abs(y)); a=0; out1 delta=10; if isempty(y1k_1) y1k_1=0; end if isempty(y2k_1) y2k_1=0; end y1k=y1k_1+h*y2k_1; % hfst
persistent v1_1 v2_1 h=0.1; delta=10; if isempty(v1_1) v1_1=0; end if isempty(v2_1) v2_1=0; end %hfst