close all; clear all;clc; monakveti_x=[-2 0.5]; monakveti_y=[0.5 1]; rad=1; wrewiris_centri_x=-0.5; wrewiris_centri_y=0.5; fi=0:0.01:2*pi; wrew_x=wrewiris_centri_x+rad*cos(fi); wrew_y=wrewiris_centri_y+rad*sin(fi); plot(monakveti_x,monakveti_y) axis equal grid on hold on plot(wrew_x,wrew_y) x_1=monakveti_x(1)-wrewiris_centri_x; x_2=monakveti_x(2)-wrewiris_centri_x; y_1=monakveti_y(1)-wrewiris_centri_y; y_2=monakveti_y(2)-wrewiris_centri_y; r=rad; [x,y,result]=circlelineintersection(x_1,x_2,y_1,y_2,r) x=x+wrewiris_centri_x; y=y+wrewiris_centri_y; plot(x,y,'o')