function y=modelnafunkcija(par,t,perioda) %y=modelnafunkcija(par,t,perioda) %izracuna vrednost funkcije, s katero modeliramo kolicino co2 v atmoesferi %glede na parametre 'par' ob casu 't'. 'perioda' je fiksen parameter, %ki ga podamo in naj bi predstavljal periodo soncnih peg % / computes the value of the function used to model the concentration %of co2 a time 't'. 'par' contains the values of the parameters, 'perioda' %is an additional parameter a=par(1);%posebej izpisemo parametre, ki nastopajo v modelu / we write out the parameters b=par(2); c=par(3); d=par(4); e=par(5); f=par(6); g=par(7); h=par(8); k=par(9); %l=par(10); t=t-1979; y=a+b*(t)+c*cos(2*pi*t)+d*sin(2*pi*t)+e*cos(4*pi*t)+f*sin(4*pi*t)+g*cos(2*pi*t/perioda)+h*sin(2*pi*t/perioda)+k*t.^2;