function [y,z] = AplicacionPoincare(xyz, x0)


id = find(xyz(2:end,1) > x0 & xyz(1:end-1,1) <= x0);

P1 = xyz(id,:);
P2 = xyz(id+1,:);

t = (x0 - P1(:,1))./(P2(:,1) - P1(:,1));

y = P1(:,2) + t.*(P2(:,2) - P1(:,2));
z = P1(:,3) + t.*(P2(:,3) - P1(:,3));


end

Published with MATLAB® 7.12