clc
close all
clear all


% Exercise 2
A = [1 exp(.5) 0; 3 4 .2]
B = [0 1; 3 -2; log(5) 4]
size(A)
size(B)
A(2, 2) = -3
A(:, 2) = []
B(1, :) = []
f = @(x)(x^(1/3) + 1).^2;
figure
fplot(f)
xlabel('x')
ylabel('y')
grid
D = f(B)

% Exercise 3
X = -3:.1:3;
Y = X;
[x, y] = meshgrid(X, Y);
q = @(x1, x2) x1.^2 + x2.^2 - 4*x1.*x2
z = q(x, y);
figure
mesh(x, y, z)
xlabel("x_1")
ylabel("x_2")
zlabel("x_3")