24/04/25 07:33:36.44 PiWgohuV.net
>>484
複素点 a, b, c, dでa,b や c,dが
２直線を形成しない座標であったり、平行なときを場合分けして
a,bを結ぶ直線とc,dを結ぶ直線の交点を返す関数を修正。
intsect[a_,b_,c_,d_] :=(
If[(a-b)(c-d)==0,Return["Not two lines."]];
;
a1=Re[a] ; a2=Im[a];
b1=Re[b] ; b2=Im[b];
c1=Re[c] ; c2=Im[c];
d1=Re[d] ; d2=Im[d];
;
mxn11=Det[{{a1,a2},{b1,b2}}];
mxn12=a1-b1;
mxn21=Det[{{c1,c2},{d1,d2}}];
mxn22=c1-d1;
mxn=Det[{{mxn11,mxn12},{mxn21,mxn22}}];
mxd=Det[{{a1-b1,a2-b2},{c1-d1,c2-d2}}];
;
If[mxd==0,Return["Two lines are pararell."]];
x=mxn/mxd;
;
myn11=mxn11;
myn12=a2-b2;
myn21=nxn21;
myn22=c2-c2;
myn=Det[{{myn11,myn12},{myn12,myn22}}];
myd=mxd;
y=myn/myd;
;
x+y*I
)
intsect[0,0,1,1I]
intsect[0I,1+0I,0+1I,1+1I]
intsect[0I,2+0I,1+1I,1+2I]
intsect[0I,2I,-1+1I,1+1I]
intsect[0,1+2I,3+4I,5+6I]