function v=fdms(primer, h)
close
%% ucitavamo podatke
run(primer)
%% formiramo mrezu
x=a:h:b;
n=length(x);
%% matrica sistema
%poddijagonala
a=-1/h^2-p(x)/(2*h);
a(n)=0;
%dijagonala
b=2/h^2+q(x);
b(1)=1;b(n)=1;
%naddiajgonala
c=-1/h^2+p(x)/(2*h);
c(1)=0;
% spdiags, na zalost radi tako da uzima prvi element vektora a, i izbacuje
% poslednji. Kod vektora c uzima poslednji element ali izbacuje prvi.
% Zbog toga moramo pomeriti elemente
a=[a(2:n),a(1)];
c=[c(n),c(1:n-1)];
M=spdiags([a',b',c'],-1:1,n,n);
% full(M) ispisaće čitavu matricu M
disp(M) % ispisuje samo poddijagonalne,dijagonalne i naddijagonalne elemente
%spy(M)
F=f(x');
F(1)=A;
F(n)=B;
v=M\F;
if nargout==0
plot(x,ue(x),'b')
hold on
plot(x,v,'--or')
hold off
end
end
