Matlab - ¿Ayuda con programa en Matlab?

Hola gente, me gustaria que me ayudaran con un programa el metodo de Jacobi Y el de Gauss-seidel en Matlab, tengo esto pero no se que le falta :
--------------------------------------…
clc, clear all, close all
%% ingreso de datos
A= input('ingrese la matriz de coheficientes')
B= input('ingrese el vector de soluciones')
%% ingreso de calculos
N=size(A);
n=N(1,1);
M=eye(n);
Mu=eye(n);

for j=1:n-1
M=eye(n)
for i=j+1:n
M(i,j)=-A(i,j)/A(j,j);
end
Mu=M*Mu;
U=M*A;
A=U;
end
U
L=inv(Mu)
Y=Mu*B
X=inv(U)*Y