Calcular el punto de interseccion de la recta perpendicular y que corta a otra que pasa por un punto con SAGE
var('u p t pr')
pr=vector(QQ,[2,0,1])# punto de la recta
u=vector(QQ,[0,1,0]) #vector direccion de la recta
p=vector(QQ,[-1,5,-3]) #punto donde pasa la perpendicular de la recta
def rectaperpendicularpunto(u1,p1,pr1,t1):
o1=p1-pr1
o2=o1.cross_product(u1)
o3=o2.cross_product(u1)
o4=[p1[0]+o3[0]*t1,p1[1]+o3[1]*t1,p1[2]+o3[2]*t1]
return(o4)
p1=rectaperpendicularpunto(u,p,pr,t)
p2=[pr[0]+u[0]*t2,pr[1]+u[1]*t2,pr[2]+u[2]*t2]
eq1=p1[0]==p2[0]
eq2=p1[1]==p2[1]
eq3=p1[2]==p2[2]
s1=solve([eq1,eq2,eq3],t,t2,solution_dict=True)
s2=s1[0][t]
s3=s1[0][t2]
p3(t)=p1
p4(t2)=p2
show(p3(s2))
show(p4(s3))
Añadir nuevo comentario