O restart;
2
3
O f1 := x1 Cx2 C7; f2 := x1 Cx2 C1;
2
3
f1 := x1 Cx2 C7
f2 := x1 Cx2 C1
2
2
O G := f1 Cf2 ;
2
3
G := x1 Cx2 C7
2
(1)
C x1 Cx2 C1
2
(2)
O plot3d G, x1 = K3 ..3, x2 = K3 ..3 ;
-3
1,500
-1
x1
1,000
1
3
500
0
-3
O e1 := f1 = 0; e2 := f2 = 0;
-2
-1
0
x2
1
2
e1 := x12 Cx23 C7 = 0
e2 := x1 Cx2 C1 = 0
O fsolve e1, e2 , x1, x2 ;
x2 = K2.000000000, x1 = 1.000000000
O f11 := diff f1, x1 ; f12 := diff f1, x2 ;
f11 := 2 x1
2
f12 := 3 x2
O f21 := diff f2, x1 ; f22 := diff f2, x2 ;
3
(3)
(4)
(5)
f21 := 1
f22 := 1
(6)
O x10 := 0; x20 := K2;
O
O
O
O
O
O
O
O
O
x10 := 0
x20 := K2
e1 := subs x1 = x10, x2 = x20 , f11 $ s1 Csubs x1 = x10, x2 = x20 , f12 $ s2 = Ksubs
= x10, x2 = x20 , f1 ;
e1 := 12 s2 = 1
e2 := subs x1 = x10, x2 = x20 , f21 $ s1 Csubs x1 = x10, x2 = x20 , f22 $ s2 = Ksubs
= x10, x2 = x20 , f2 ;
e2 := s1 Cs2 = 1
sol1 := solve e1, e2 , s1, s2 ;
1
11
sol1 := s2 =
, s1 =
12
12
11
1
x10 := evalf 0 C
; x20 := evalf K2 C
;
12
12
x10 := 0.9166666667
x20 := K1.916666667
e1 := subs x1 = x10, x2 = x20 , f11 $ s1 Csubs x1 = x10, x2 = x20 , f12 $ s2 = Ksubs
= x10, x2 = x20 , f1 ;
e1 := 1.833333333 s1 C11.02083334 s2 = K.799189811
e2 := subs x1 = x10, x2 = x20 , f21 $ s1 Csubs x1 = x10, x2 = x20 , f22 $ s2 = Ksubs
= x10, x2 = x20 , f2 ;
e2 := s1 Cs2 = K0.
sol1 := solve e1, e2 , s1, s2 ;
sol1 := s2 = K0.08698664603, s1 = 0.08698664603
11
1
x10 := evalf 0 C
C0.0869866 ; x20 := evalf K2 C
K0.0869866 ;
12
12
x10 := 1.003653267
x20 := K2.003653267
(7)
x1
(8)
x1
(9)
(10)
(11)
x1
(12)
x1
(13)
(14)
(15)