% 
clear all
% Heatequation semidiscterized with Robin boundary conditions
%
global m eta1 eta0
eta1=1000; eta1=0;
eta0=100; eta0=0;
m=100;
mm=m+2;
xx=[0:1/(m+1):1];
T=1/2;
y0=2000*(1+sin(2*pi*xx))'; u0=y0;

plot(xx,u0)
pause
options=odeset('AbsTol',1e-10,'RelTol',1e-9);
[tt, yy]=ode15s('funcHeNeu',[0 T],y0,options);
for i=1:size(yy,1)
    plot(xx,yy(i,:))
    drawnow
end    
mesh(xx,tt,yy)
xlabel('rom')
ylabel('tid')