% Řešení okrajové úlohy -(p(x)u')' + q(x)u = f(x)
% Příklad 1: -(cos(x)u')'+sin(x)u=x^2, u(0)=0, u(pi/2)=0
p=@(x)cos(x);
q=@(x)sin(x);
f=@(x)x.^2;
I=[0,pi/2];
BV=[0,0];
[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,':')