echo on
%
% Numericke chyby v Matlabu
%
% Pr. 1: Chyba v reprezentaci cisel
pause
c=0.1;              % chyba je asi na 16. desetinnem miste
pause
u=c*ones(1,100);    % vektor u obsahuje 100 krat jednu desetinu
pause
sum(u)-10           % vysledek by mel byt nula
pause
%
% Vypocet preorganizujeme
%
A=reshape(u,10,10); % preskaldame vektor u do matice
pause
sum(sum(A))-10      % vysledek je asi 10 krat presnejsi nez predtim
pause
% Nasleduje prikaz 'keyboard', pro pokracovani skriptu zadejte prikaz 'dbcont'.
keyboard
clc
%
% Pr.2: Nasobne koreny polynomu
%
format long
P=poly(1:4)         % polynom s koreny 1,2,3,4
pause
roots(P)            % koreny jsou nalezeny pomerne presne
pause
P=poly([1 1 1 1])   % polynom se ctyrnasobnym korenem rovnym 1
pause
roots(P)            % tady je chyba pomerne velka
pause
format short
% Nasleduje prikaz 'keyboard', pro pokracovani skriptu zadejte prikaz 'dbcont'.
keyboard
clc
%
% Pr.3: Problemy s grafickym zobrazenim
f=@(x) x.^3-9*x.^2+27*x-27;       % funkce (x-3)^3
pause
xx=linspace(3-0.00005,3+0.00005,201); % body pro zobrazeni v okoli bodu 3, velikosti okoli je 10^-4
pause
fx=f(xx);       % funkcni hodnoty
pause
plot(xx,fx)     % graf funkce
pause
close
%
% Totez pomoci polynomu
%
P=[1 -9 27 -27];
pause
Px=polyval(P,xx);   % hodnoty polynomu
pause
plot(xx,Px)         % je to trochu lepsi
pause
close
%
% Zkusime upravit funkci
%
f=@(x) (x-3).^3;
pause
fx=f(xx);           % funkcni hodnoty
pause
plot(xx,fx)         % je to vyrazne lepsi
pause
% Nasleduje prikaz 'keyboard', pro pokracovani skriptu zadejte prikaz 'dbcont'.
keyboard
close
clc
echo off
clear all