function vk = fun1(tk)
% vk = fun1(tk)
% 
% Kriterijska funkcija: Simulacija vertikalnega padanja padalca v odvisnosti od zacetne hitrosti
% in zracnega upora.
%
% Vhod:
% tk je koncni cas, do katerega gledamo padanje (pred odprtjem padala!!!)
%
%
% Izhod:
% vk koncna hitrost
%
% (Definirano kriterijsko funkcijo fun1 uporabimo na fzero, da poiscemo niclo, npr.: fzero(@(x)fun1(x),30))
%

caseNo = 3;
m = 105;
c = 1;
S = 1.2;
n = 10000;
zac = [40000; 0];
parametri = [m c S];

t = linspace(0,tk,n);
dY = @(t,Y) diferencialniSistem(t,Y,parametri);
[~,Yres] = ode45(dY,t,zac);

% visine in hitrosti padalca

vk = Yres(end,2)+300;

function dY = diferencialniSistem(t,Y,parametri)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end

end