%% kmin/k_min/kMin per tolleranza desiderata SEZIONE AUREA
phi = (1+sqrt(5))/2;
% errStimato = 2/phi^k
k = ceil( log(2/tol) / log(phi))
%% Calcolo iterazioni bisezione dato errore prefissato
clc; clear; close all;
% numero iterazioni per errore prefissato
syms k; assume(k,'real');
sol = solve(err_k == (b-a)/(2^(k+1)), k);
disp(simplify(sol));
%% Matrice tridiagonale
elementi = [1, 2, 3];
n = 6;
A = diag(elementi(1)*ones(n-1,1), -1) + ...
diag(elementi(2)*ones(n, 1), 0) + ...
diag(elementi(3)*ones(n-1,1), 1)
%% Matrice pentadiagonale
elementi = [1, 2, 3, 4, 5];
n = 6;
A = diag(elementi(1)*ones(n-2,1), -2) + ...
diag(elementi(2)*ones(n-1,1), -1) + ...
diag(elementi(3)*ones(n, 1), 0) + ...
diag(elementi(4)*ones(n-1,1), 1) + ...
diag(elementi(5)*ones(n-2,1), 2)
%% Simmetrica e Definita Positiva
if(isequal(A, A'))
disp('simmetrica')
end
if all(eig(A) > 0)
disp('definita positiva')
end
