ド底辺シリツ医には解けないであろう問題。

https://i.imgur.com/iBRWjep.jpg

プログラムして２０度以外でも答がでるようにしてみた。

intsect <- function(a,b,c,d){
a1=Re(a) ; a2=Im(a)
b1=Re(b) ; b2=Im(b)
p=(a2-b2)/(a1-b1)
c1=Re(c) ; c2=Im(c)
d1=Re(d) ; d2=Im(d)
q=(c2-d2)/(c1-d1)
if(p==q) return(NA)
else{
x= ((p*a1 - a2) - (q*c1 - c2))/ (p-q)
y= p*x - ( p*a1 - a2)
x + 1i*y
}
}
foo <- function(N20=20){

plot(NULL,type='n',xlim=c(0,2),ylim=c(0,2),asp=1,ann=F)