0414卵の名無しさん
2020/08/08(土) 20:21:06.61ID:Fonqs1WQhttps://i.imgur.com/iBRWjep.jpg
プログラムして20度以外でも答がでるようにしてみた。
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)