Luento 2, ti 25.2.05

clear all; close all
n=5; h=1/(n+1);
   % n sisäpistettä, n+1 jakoväliä, n+2 jakopistettä reunoineen. 
A=-diag(ones(n-1,1),-1)+2*diag(ones(n,1))-diag(ones(n-1,1),1);

% Isommilla n:n arvoilla on syytä tallettaa harvana, spdiags
%  Kts. help sparse
x=h:h:1-h;   % Sisäsolmut (n kpl).
f=inline('(3*x+x.^2).*exp(x)','x')
%u=inline('','') % Hae tämä Maplesta.%
b=h^2*f(x); b=b'  % Oltava pystyvektori.
v=A\b;
x=[0,x,1]; v=[0;v;0];
%fplot(u,[0,1],'r'); hold on
plot(x,v,'o')
axis([-.1 1.1 -.1 .5])
grid on