欢迎使用CSDN-markdown编辑器

function [tout,yout]= f_cal_kernel(eqs,ini_cdt,tstart,tfinal)    refine = 4;    options = odeset('RelTol',1e-5,'AbsTol',1e-6,...        'Events', @f_cal_kernel_impact_events,...        'InitialStep',1e-6, 'MaxStep',1e-1,...        'Refine',refine);    y0 =ini_cdt;    tout = tstart;    yout = y0;    teout = [];    yeout = [];    ieout = [];    t_mini=1e-2;    while t_mini<(tfinal-tstart)        [t,y,te,ye,ie] = ode23tb(eqs,[tstart tfinal],y0.',options);        nt = length(t);        tout = [tout; t(2:nt)];        yout = [yout; y(2:nt,:)];        teout = [teout; te];        yeout = [yeout; ye];        ieout = [ieout; ie];        y0(1) = 1e-6;        y0(2) = -0.8*y(nt,2);        y0(3) = y(nt,3);        tstart = t(nt);    endend