valor de contorno

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

N=20
h=.1
x=[0:h:2]'
y=zeros(N+1,1)
A = zeros(N+1,N+1)
b = zeros(N+1,1)
A(1,1) = 1
b(1) = 0
for i=2:N
 A(i,i-1) = -1/h^2
 A(i,i) = 2/(h^2) + 1/10000
 A(i,i+1) = -1/h^2
 
 b(i) = 100
end
A(N+1,N+1) = 1
b(N+1) = 10
y = inv(A)*b
disp([x,y])