% TUTORIAL 8

% Question 1
% Define system variables
syms x y
% Define the function
f1 = (x^3) * (y^2) * (12-x-y)

f1 =

% First-order partial derivatives

disp('First-order partiral derivatives: ')

First-order partiral derivatives:

p = diff(f1, x)

p =

q = diff(f1, y)

q =

% Second-order partial derivatives

disp('Second-order partial derivatives: ')

Second-order partial derivatives:

r = diff(p, x)

r =

s = diff(p, y)

s =

t = diff(q, y)

t =

% Solve for critical points

[critical_point_x, critical_point_y] = solve([p == 0, q == 0], [x, y]);
critical_points = [critical_point_x, critical_point_y]

critical_points =

1
for i = 1:length(critical_point_x)
% Current critical point
current_points = [critical_point_x(i), critical_point_y(i)]

% Evaluate second-order derivatives at the critical point

r1 = subs(r, {x, y}, current_points)
s1 = subs(s, {x, y}, current_points)
t1 = subs(t, {x, y}, current_points)

% Hessian determinant
disp('h = r*t - (s^2)')
h1 = r1 * t1 - (s1^2)

% Classification of critical point

if (h1 > 0 && r1 > 0)
disp('The function has a local minimum at this point.');
elseif (h1 > 0 && r1 < 0)
disp('The function has a local maximum at this point.');
elseif (h1 < 0)
disp('The function has a saddle point at this point.');
else
disp('The second derivative test is inconclusive. (Test Fails)');
end
end

current_points =
r1 =
s1 =
t1 =
h = r*t - (s^2)
h1 =
The function has a local maximum at this point.
current_points =
r1 =
s1 =
t1 =
h = r*t - (s^2)
h1 =
The second derivative test is inconclusive. (Test Fails)

% Question 2
% Define system variables
syms x y
% Define the function
f2 = (x^2)+(y^2)-6*x+12

f2 =

2
% First-order partial derivatives
disp('First-order partiral derivatives: ')

First-order partiral derivatives:

p = diff(f2, x)

p =

q = diff(f2, y)

q =

% Second-order partial derivatives

disp('Second-order partial derivatives: ')

Second-order partial derivatives:

r = diff(p, x)

r =

s = diff(p, y)

s =

t = diff(q, y)

t =

% Solve for critical points

[critical_point_x, critical_point_y] = solve([p == 0, q == 0], [x, y]);
critical_points = [critical_point_x, critical_point_y]

critical_points =

for i = 1:length(critical_point_x)
% Current critical point
current_points = [critical_point_x(i), critical_point_y(i)]

% Evaluate second-order derivatives at the critical point

r2 = subs(r, {x, y}, current_points)
s2 = subs(s, {x, y}, current_points)
t2 = subs(t, {x, y}, current_points)

% Hessian determinant
disp('h = r*t - s^2')
h2 = r2 * t2 - (s2^2)

3
 % Classification of critical point
if (h2 > 0 && r2 > 0)
disp('The function has a local minimum at this point.');
elseif (h2 > 0 && r2 < 0)
disp('The function has a local maximum at this point.');
elseif (h2 < 0)
disp('The function has a saddle point at this point.');
else
disp('The second derivative test is inconclusive. (Test Fails)');
end
end

current_points =
r2 =
s2 =
t2 =
h = r*t - s^2
h2 =
The function has a local minimum at this point.

% Question 3
% Define system variables
syms x y
% Define the function
f3 = (x^2)+(y^2)-x*y-2*x-y

f3 =

% First-order partial derivatives

disp('First-order partiral derivatives: ')

First-order partiral derivatives:

p = diff(f3, x)

p =

q = diff(f3, y)

q =

% Second-order partial derivatives

disp('Second-order partial derivatives: ')

Second-order partial derivatives:

r = diff(p, x)

4
r =

s = diff(p, y)

s =

t = diff(q, y)

t =

% Solve for critical points

[critical_point_x, critical_point_y] = solve([p == 0, q == 0], [x, y]);
critical_points = [critical_point_x, critical_point_y]

critical_points =

for i = 1:length(critical_point_x)
% Current critical point
current_points = [critical_point_x(i), critical_point_y(i)]

% Evaluate second-order derivatives at the critical point

r3 = subs(r, {x, y}, current_points)
s3 = subs(s, {x, y}, current_points)
t3 = subs(t, {x, y}, current_points)

% Hessian determinant
disp('h = r*t - s^2')
h3 = r3 * t3 - (s3^2)

% Classification of critical point

if (h3 > 0 && r3 > 0)
disp('The function has a local minimum at this point.');
elseif (h3 > 0 && r3 < 0)
disp('The function has a local maximum at this point.');
elseif (h3 < 0)
disp('The function has a saddle point at this point.');
else
disp('The second derivative test is inconclusive. (Test Fails)');
end
end

current_points =

r3 =
s3 =

5
t3 =
h = r*t - s^2
h3 =
The function has a local minimum at this point.

% Question 4
% Define system variables
syms x y
% Define the function
f4 = (x^2) * (y^3) * (1-x-y)

f4 =

% First-order partial derivatives

disp('First-order partiral derivatives: ')

First-order partiral derivatives:

p = diff(f4, x)

p =

q = diff(f4, y)

q =

% Second-order partial derivatives

disp('Second-order partial derivatives: ')

Second-order partial derivatives:

r = diff(p, x)

r =

s = diff(p, y)

s =

t = diff(q, y)

t =

% Solve for critical points

[critical_point_x, critical_point_y] = solve([p == 0, q == 0], [x, y]);
critical_points = [critical_point_x, critical_point_y]

6
critical_points =

for i = 1:length(critical_point_x)
% Current critical point
current_points = [critical_point_x(i), critical_point_y(i)]

% Evaluate second-order derivatives at the critical point

r4 = subs(r, {x, y}, current_points)
s4 = subs(s, {x, y}, current_points)
t4 = subs(t, {x, y}, current_points)

% Hessian determinant
disp('h = r*t - s^2')
h4 = r4 * t4 - (s4^2)

% Classification of critical point

if (h4 > 0 && r4 > 0)
disp('The function has a local minimum at this point.');
elseif (h4 > 0 && r4 < 0)
disp('The function has a local maximum at this point.');
elseif (h4 < 0)
disp('The function has a saddle point at this point.');
else
disp('The second derivative test is inconclusive. (Test Fails)');
end
end

current_points =

r4 =

s4 =

t4 =

h = r*t - s^2
h4 =

The function has a local maximum at this point.

7
current_points =
r4 =
s4 =
t4 =
h = r*t - s^2
h4 =
The second derivative test is inconclusive. (Test Fails)

% Question 5
% Define system variables
syms x y
% Define the function
f5 = (x^3)+x*(y^2)-12*(x^2)+21*x-2*(y^2)

f5 =

% First-order partial derivatives

disp('First-order partiral derivatives: ')

First-order partiral derivatives:

p = diff(f5, x)

p =

q = diff(f5, y)

q =

% Second-order partial derivatives

disp('Second-order partial derivatives: ')

Second-order partial derivatives:

r = diff(p, x)

r =

s = diff(p, y)

s =

t = diff(q, y)

t =

% Solve for critical points

8
[critical_point_x, critical_point_y] = solve([p == 0, q == 0], [x, y]);
critical_points = [critical_point_x, critical_point_y]

critical_points =

for i = 1:length(critical_point_x)
% Current critical point
current_points = [critical_point_x(i), critical_point_y(i)]

% Evaluate second-order derivatives at the critical point

r5 = subs(r, {x, y}, current_points)
s5 = subs(s, {x, y}, current_points)
t5 = subs(t, {x, y}, current_points)

% Hessian determinant
disp('h = r*t - s^2')
h5 = r5 * t5 - (s5^2)

% Classification of critical point

if (h5 > 0 && r5 > 0)
disp('The function has a local minimum at this point.');
elseif (h5 > 0 && r5 < 0)
disp('The function has a local maximum at this point.');
elseif (h5 < 0)
disp('The function has a saddle point at this point.');
else
disp('The second derivative test is inconclusive. (Test Fails)');
end
end

current_points =
r5 =
s5 =
t5 =
h = r*t - s^2
h5 =
The function has a local maximum at this point.
current_points =
r5 =
s5 =
t5 =
h = r*t - s^2
h5 =
The function has a local minimum at this point.

9
current_points =
r5 =

s5 =
t5 =
h = r*t - s^2
h5 =
The function has a saddle point at this point.
current_points =
r5 =

s5 =
t5 =
h = r*t - s^2
h5 =
The function has a saddle point at this point.

10