findzeros.m

来自「Resolution of the Cronig-Penney s predat」· M 代码 · 共 13 行

function [yz, indices] = findzeros(y)
    
    k = 1;
    yz = [];
    indices = [];
    for j = 1 : length(y) - 1
        if y(j) * y(j + 1) <= 0
            yz(k) = y(j + 1);
            indices(k) = j + 1;
            k = k + 1;
        end
    end