利用solve函数求解直线上到某点的距离为L的点的座标

%% 利用solve函数求解直线上到某点的距离为L的点的座标
P0=[4 5 6];%定义点
v1=[1 2 3];%定义直线向量
L=-10;%定义到固定点P的距离，负号表示方向
syms t
t=solve('(v1(1)^2+v1(2)^2+v1(3)^2)*t^2-100','t')%利用solve函数求解非线性方程,t为参数式直线方程的参数
t=eval(t)%将符号解转化为数值解
if t(1)*L>0
   t=t(1);
else
   t=t(2);
end
x=P0(1)+v1(1)*t;%直线的参数方程，将t代入，可得点的座标，下同
y=P0(2)+v1(2)*t;
z=P0(3)+v1(3)*t;
P1=[x y z]

计算结果为：

t =

   -2.6726

P1 =

    1.3274   -0.3452   -2.0178