%Zadanie 1
A = [3.15, 1.32, 4.58, 4.57;
     4.06, -4.02, 4.65, -0.15;
     -3.75, -2.22, -3.42, 3;
     4.13, 0.47, 4.71, -3.58]; % macierz wartości przy niewiadomych
B = [38.72 ; -6.94 ; -0.41 ; -6.02]; % macierz rozwiązań
X = (A^(-1) * B);  % macierz wynikowa
disp(X);
%Zadanie 2
x1 = [-(pi./4) : 0.1 : (pi./4)];
x2 = [(pi/4) : 0.1 : 3 .* exp(1)];
y1 = tan(x1).^2;
y2 = log(x2) - x2.^2 + x2 - 1;
plot(x1,y1);
title('Wykres f(x)');
xlabel('x');
ylabel('y');
hold on;
plot(x2,y2);
legend('tg^2(x)', 'ln(x)-x^2+x-1');
hold on;
%Zadanie 3
x1 = linspace(0, (pi./5), 100)
y1 = linspace(0, exp(1)./5, 100)
% lub:
%x1 = 0 : 0.1 : (pi./5);
%y1 = 0 : 0.1 : exp(1)./5;
[ x, y ] = meshgrid(x1,y1);
f = ((sin(x.*y).^2 + x.^y - 6 .* x - y)) ./ ((exp(1)).^y - cot(x .* y));
surf(x,y,f);
xlabel('x');
ylabel('y');
zlabel('z');
title('Wykres f(xy)');