derivação

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

//malha
N = 11
h = 1/(N-1)
x = 0:h:1
//sistema discreto
A = zeros(N,N)
b = zeros(N,1)
A(1,1) = 1
b(1) = 1
for i = 2:N-1
 A(i,i-1) = -1/h^2
 A(i,i) = 2/h^2 + 1
 A(i,i+1) = -1/h^2
 
 b(i) = exp(-x(i))
end
A(N,N) = 1
b(N) = 2
u = inv(A)*b