% Řešení okrajové úlohy -(p(x)u')' + q(x)u = f(x)
% Příklad 1: -(x*u')'+x^2*u=sin(x)cos(x), u(0)=1, u(1)=2
p=@(x)x;
q=@(x)x.^2;
f=@(x)sin(x).*cos(x);
I=[0,1];
BV=[1,2];
[x1,u1]=solveBP_FEM(p,q,f,I,BV,100);    % přibližné řešení 1
[x2,u2]=solveBP_FEM(p,q,f,I,BV,1000);   % přibližné řešení 2
[x3,u3]=solveBP_FEM(p,q,f,I,BV,10000);  % přibližné řešení 3
plot(x1,u1,x2,u2,'--',x3,u3,':')