利用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