n=100;
p=primes(n); %row vector
X=zeros(n-1,length(p));
for i=1:n-1
   v=factor(i+1);
   for j=1:length(v)
      X(i,:)=X(i,:)+(p==v(j));
   end
end
lcm=1;
for i=1:length(p)
    lcm=lcm*p(i).^(max(X(:,i)));
end
lcm