%已有函數 g1=fouryie(x) g=fouradd(x, y) g=fourplus(x, y)；使用函數先定義符號變量。
syms g0 g1 g2 g3 q0 q1 q2 q3 x1 x2;
%%%%%%%%%%%%%%%%%%%%已知數據
A1=[0;0;0;0];A2=[0;-1.8301;3.1699;-0.5];A3=[0;3.1699;11.8301;0.5];
A4=[0;6.8301;11.8301;-1];A5=[0;11.8301;3.1699;0];A6=[0;10;0;-1];
B1=[0;0;0;0];B2=[0;-2;3.6049;0.069];B3=[0;-0.8349; 5.4247;-0.9359];
B4=[0;3.1651;5.6835;0.03];B5=[0;4.3302;3.3462;-0.8966];B6=[0;2.3302;0.2589;0.9659];
V1=[0;0.965923;-0.25883;0];V2=[0;0.6;0.8;0];
L1=9.8080;L2=10.8465;L3=11.0578;L4=12.8826;L5=9.6551;L6=12.3934;
%%%%%%%%%%%%%%%%%%%%%%%%%%
%A=A1;B=B1;L=L1;V0=x1*V1;VJ=V1;%改變數據求相應表達式；
A=A2;B=B2;L=L2;V0=x2*V2;VJ=V2;
%A=A3;B=B3;L=L3;
%A=A4;B=B4;L=L4;
%A=A5;B=B5;L=L5;
%A=A6;B=B6;L=L6;
%%%%%%%%%%%%%%%%
g=[g0; g1; g2; g3];q=[q0; q1;q2;q3];
Temp0=fourTime(q, fouryie(q));
Temp1=fourTime(g, fouryie(q));
Temp2=fourTime(q, fourTime(B, fouryie(q)));
Temp3=fouradd(Temp1, Temp2);
Temp4=fourplus(Temp3,A);
Temp5=fourTime(Temp4, fouryie(Temp4));
Temp60=fourplus(Temp4,V0);
Temp6=fourTime(Temp60,fouryie(Temp60));
Temp7=fourTime(Temp60,fouryie(VJ));
if A==A1
  f1=Temp0(1)-1
  f2=Temp1(1)
else
end
if A==A1|A==A2;
    f3=Temp6(1)-L^2
    f4=Temp7(1)
else
    f3=Temp5(1)-L^2
end
%%%%%%%%%%%%%%%%%%%%