function v=fdm1(primer, h)
%% 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;
M=diag(b)+diag(a(2:n),-1)+diag(c(1:n-1),1);
%disp(M);
F=f(x');
F(1)=A;
F(n)=B;
v=M\F; % standardno rešavamo sistem
if nargout==0
plot(x,ue(x),'b')
hold on
plot(x,v,'--or')
hold off
end
end
