clc;
clear all;
close all;
A=zeros(5,5)
for j=1:5
    for i=1:5
        if i>3
            A(i,j)=sin(pi/i)^2-0.94;
        else
            A(i,j)=((i-j)^2*0.983)^0.5;
        end
    end
end
A
K7=length(A(1,:))
SrZn=sum(A')'/K7
MIN=min(A)
MAX=max(A)
A1=A;
A1(:,[1 4])=A1(:,[4 1])
X=diag(A)
P = sort(X)
%% SLAU
A6 = [3 -1 0;-2 1 1;2 -1 4]
B=[5;0;15]
X=A6^-1*B
Z=A6*X


%%zadanie2 interpolyacia-approximacia
X=[0.324 0.645 0.966 1.287 1.609 1.930 2.251 2.572 2.893]
Y=[-2.052 -1.756 -1.076 -0.284 0.982 2.209 4.013 5.796 8.011]
plot(X,Y, '.r', 'MarkerSize', 20)
hold on
grid on
xx=X(1):0.1:X(end);
%%xx=X(1):0.1:2.924;
yy=spline(X,Y,xx);
figure(1)
plot(xx,yy)


%%zadanie3 MNK
X=[0.324 0.645 0.966 1.287 1.609 1.930 2.251 2.572 2.893]
Y=[-2.052 -1.756 -1.076 -0.284 0.982 2.209 4.013 5.796 8.011]

figure(2)
plot(X,Y, '*r', 'MarkerSize', 10)
hold on
grid on
xxx=X(1):0.1:2.924;
p=polyfit(X,Y,2)
yyy=p(1)*xxx.*xxx+p(2)*xxx+p(3);
figure(2)
plot(xxx,yyy)

%%zadanie4 integral
X=[0.324 0.645 0.966 1.287 1.609 1.930 2.251 2.572 2.893]
Y=[-2.052 -1.756 -1.076 -0.284 0.982 2.209 4.013 5.796 8.011]
Q=trapz(X,Y)
%%zadanie4 integral
Q1=quad('exp(-(X.^2)/2)',-2,2,1e-4)