2. Cvičení

Materiály ze cvičení

Bohužel může zase nastat problém s přílišným počtem stažení.

Chcete-li vyzkoušet svou funkci na řešení soustavy rovnic Ax=b, kód je:

n  = 5;
CC = 0;DD = 0;
for kk = 1:1e4
 A = round(5*rand(n));
 b = round(5*rand(n,1));
 A = A - diag(diag(A));

  %rank check
 if rank(A)~=n;DD=DD+1;end;

 %%%%%%%%%%%%%%%%%%%%%%%%%%
 %%YOUR CODE HERE%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%

  % solution check
 xM = A\b;
 if round(1e3*xM)~= round(1e3*xYours);CC=CC+1;end;
end
% number of fails
disp(CC-DD); 

K dobrovolnému samostudiu

Existence a jednoznačnost řešení

Pro ty z vás, které ještě pořád zajímá proč má soustava obyčejných diferenciálních rovnic za předpokladu, že vektorové pole (pravá strana) je spojité vždycky řešení a proč je toto řešení při zavedení počáteční podmínky jednoznačné doporučuji,