%% Řešení okrajového problému
% Rovnice -u'' +u=\cos x
% Okrajové podmínky u(0)=0,\ u(\frac{\pi}{2})=1

help solveBP1

%

p=@(x)-ones(size(x));
q=@(x) zeros(size(x));
r=@(x) ones(size(x));
f=@(x) cos(x);
I=[0,pi/2];BV=[0 1];
n=10;
[x,y] = solveBP1(p,q,r,f,I,BV,n);
plot(x,y)

%
%% Najděte přesné řešení a porovnejte ho s numerickým, pak zkuste totéž pro krok 0.01 (n=100).

