function Szelo_plot(P1, P2, a, b)
    %% A P1, P2 pontokon átmenő egyenes egyenlete:
    %p = (P1(2) -P2(2))/(P1(1) - P2(1))*(x - P1(1)) + P1(2)
    m = 5;   % a szelőről szép grafikon - m = 2 is jó?!
    xx = linspace(a,b,m);
    pxx = (P1(2) -P2(2))/(P1(1) - P2(1))*(xx - P1(1)) + P1(2);
    plot([P1(1),P2(1)],[P1(2),P2(2)],'markersize', 24, 'r.',xx,pxx)  %% saját színnel jobb lenne...
endfunction