diferecas-centrais

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

N = 20
h = 2/N
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) + 
 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])