Lagrangen funktion gradientti on laskettu kirjassa:
> clos28ex3:=[3*x^2+2*y*z-2*x*lambda,2*x*z-2*y*lambda,2*x*y-2*z-2*z*lambda,x^2+y^2+z^2-1]; 2 clos28ex3 := [3 x + 2 y z - 2 x lambda, 2 x z - 2 y lambda, 2 2 2 2 x y - 2 z - 2 z lambda, x + y + z - 1]
> joo:=allvalues(generalsolve(clos28ex3),`d`); joo := [{lambda = 0, y = -1, z = 0, x = 0}, {lambda = 0, y = 1, z = 0, x = 0}, {lambda = -3/2, y = 0, z = 0, x = -1}, {lambda = 3/2, y = 0, z = 0, x = 1}, {z = 1, y = 0, x = 0, lambda = -1}, {z = -1, y = 0, x = 0, lambda = -1}, {x = -2/3, lambda = -4/3, y = 1/3, z = 2/3}, {x = -2/3, lambda = -4/3, z = -2/3, y = -1/3}, 1/2 1/2 {x = -3/8, lambda = 1/8, y = - 3/16 22 , z = 1/16 22 }], [{lambda = 0, y = -1, z = 0, x = 0}, {lambda = 0, y = 1, z = 0, x = 0}, {lambda = -3/2, y = 0, z = 0, x = -1}, {lambda = 3/2, y = 0, z = 0, x = 1}, {z = 1, y = 0, x = 0, lambda = -1}, {z = -1, y = 0, x = 0, lambda = -1}, {x = -2/3, lambda = -4/3, y = 1/3, z = 2/3}, {x = -2/3, lambda = -4/3, z = -2/3, y = -1/3}, 1/2 1/2 {x = -3/8, lambda = 1/8, y = 3/16 22 , z = - 1/16 22 }] > joo2:=convert(joo[1],set) union convert(joo[2],set); joo2 := {{lambda = 0, y = 1, z = 0, x = 0}, {lambda = 0, y = -1, z = 0, x = 0}, {lambda = 3/2, y = 0, z = 0, x = 1}, {lambda = -3/2, y = 0, z = 0, x = -1}, {z = 1, y = 0, x = 0, lambda = -1}, {z = -1, y = 0, x = 0, lambda = -1}, {x = -2/3, lambda = -4/3, y = 1/3, z = 2/3}, {x = -2/3, lambda = -4/3, z = -2/3, y = -1/3}, 1/2 1/2 {x = -3/8, lambda = 1/8, y = 3/16 22 , z = - 1/16 22 }, 1/2 1/2 {x = -3/8, lambda = 1/8, y = - 3/16 22 , z = 1/16 22 }}
> solve(convert(clos28ex3,set)); {{lambda = 0, y = 1, z = 0, x = 0}, {lambda = 0, y = -1, z = 0, x = 0}, {lambda = 3/2, y = 0, z = 0, x = 1}, {lambda = -3/2, y = 0, z = 0, x = -1}, {z = 1, y = 0, x = 0, lambda = -1}, {z = -1, y = 0, x = 0, lambda = -1}, {x = -2/3, lambda = -4/3, y = 1/3, z = 2/3}, {x = -2/3, lambda = -4/3, z = -2/3, y = -1/3}, 1/2 1/2 {x = -3/8, lambda = 1/8, y = 3/16 22 , z = - 1/16 22 }, 1/2 1/2 {x = -3/8, lambda = 1/8, y = - 3/16 22 , z = 1/16 22 }}