clear all

x = linspace(0,1,20);
t = linspace(0,2,5);

m = 0;
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);

    u = sol(:,:,1);

%%
surf(x,t,u)    
    title('Numerical solution computed with 20 mesh points')
    xlabel('Distance x')
    ylabel('Time t')

%%
figure
    plot(x,u(end,:))
    title('Solution at t = 2')
    xlabel('Distance x')
    ylabel('u(x,2)')
 
%%
    xout=x(1:3)+.18;
    j=4;
    [uout,DuoutDx] = pdeval(m,x,u(j,:),xout)