Замена в матрицах (Pascal -> Matlab)

program z3614; 
uses crt;
const 
n=5; m=7; 
type matr=array[1..n,1..m] of real; 
var a,b:matr; 
 
procedure replace(n,m,k,l:integer; c:real;var u:matr); 
var i,j:integer; 
begin 
for i:=1 to n do 
for j:=1 to m do 
if ((u[i,j]=c) and (j=k)) then 
begin 
u[i,j]:=l; 
end; 
readln;
end;
 
procedure inputmatr(n,m,p:integer; x:char; var u:matr); 
var i,j:integer; 
begin 
writeln('input ',x,' ',n,' * ',m); 
for i:=1 to n do 
for j:=1 to m do 
u[i,j]:=p; 
readln;
end;
 
 
procedure outputmatr(n,m:integer; var f:matr); 
var i,j:integer; 
begin 
writeln('Matriza: '); 
for i:=1 to n do 
begin 
for j:=1 to m do 
write(f[i,j],' '); 
writeln(' '); 
end; 
end; 
 
begin 
vvodmatr(5,7,0,'a',a); 
vivodmatr(5,7,a); 
vvodmatr(4,5,1,'b',b); 
vivodmatr(4,5,b); 
zamena(5,7,3,1,0,a); 
zamena(4,5,5,0,1,b); 
vivodmatr(5,7,a); 
vivodmatr(4,5,b); 
readln; 
end.