function A2P4
% Checks TwoPtBVP

a = 0;
b = 1;
n = 500;

uvals = TwoPtBVP(n,a,b,@p,@q,@r,@s);
x = linspace(a,b,n);
plot(x,uvals)
shg

clc
% Want to see a snapshot of uvals to check for correctness...
uvals(1:40:n)


function y = p(x)
y = 2 + x*cos(20*pi*x);

function y = q(x)
y = 20*sin(1/(.1+x)^2);

function y = s(x)
y = -30+70*x*(x-1)*cos(5*pi*x);

function y = r(x)
y = 1000*(x-.3)*(x-.7);
