Help

Hi guys, I got errors at this point of running the codes. Can anyone help?
M=3
Y=[126;126;126;127;128;133;126;126;126;126;126;126;126;126;126;125;124;124;130;130]


function [h, mu, sigmas, P] = em(M, Y);

[T,D] = size(Y); % T=20, D=1
Tmu=1.0e-03; % Tmu=0.0010
h = ones(T,M)/M;
P = ones(M,1)/M;

for j=1:M,
k = round(rand(1)*T);
mu(j,1) = Y(k,1);
sigmas(j,1) = ones(1,D);
end

Which language is that. Doesn't seem to be like C or BASIC.