function lampoCoo7Neu
% Teht. 7/8
% Neumannin RE't

% Lämpöyhtälö: c=vakio, m=0, f(x,t,u,dudx)=dudx, s=0
% (AE) : u(x,t0)= u0(x)

close all
global x t u     % Ennen kutsua myös istunnossa, jos haluat jälkeenpäin käyttää.
m = 0;
L=10; x = linspace(0,L,20);
%t = linspace(0,1,5);
t = linspace(0,10,40);
sol = pdepe(m,@lampo,@AE,@RE,x,t);
% Extract the first solution component as u.

u = sol(:,:,1);

% A surface plot is often a good way to study a solution.
surf(x,t,u)    
title('Numerical solution computed with 20 mesh points.')
xlabel('Distance x')
ylabel('Time t')

%[xi,ti,Ui]=harvenna(x,t,u,3,2);
%figure
%surfc(xi,ti,Ui);
% A solution profile can also be illuminating.
figure
plot(x,u(end,:))
title(['Ratkaisu hetkellä t = ',num2str(t(end))])
xlabel('Paikka x')
ylabel('u(x,Tloppu)')
% --------------------------------------------------------------
% Lämpöyhtälö: c=1 (huomaa: c on vasemmalla puolella, eli pi^2 u_t = u_xx
function [c,f,s] = lampo(x,t,u,DuDx)
c = 1;
f = DuDx;
s = 0;
% --------------------------------------------------------------
function u0 = AE(x)
%u0 = exp(-(x-5).^2);    % Teht. 7
u0 = .01*(x.^4/4-17*x.^3/3+35*x.^2)-1;    % Teht. 7
% --------------------------------------------------------------
function [pl,ql,pr,qr] = RE(xl,ul,xr,ur,t)
% Neumannin RE:t  
pl = 0;                    % Vasemmassa reunassa derivaatta = 0
ql = 1;     
%pr=0;
pr = 0;                   % Oikeassa reunassa derivaatta u_x = t  
qr = 1;                   
