Matlab - generacion matriz

function [yp]=yprimas(y,t)

% y(1)= x1

% y(2)= x1´

% y(3)= x2

% y(4)= x2´

y=[2 3 4 5]

mu=1/82.45;

nu=1-mu;

yp=y;

h=201

% distancias del satélite a la tierra y a la luna,

d1=((y(1)+mu)^2+y(3)^2)^1.5;

d2=((y(1)-nu)^2+y(3)^2)^1.5;

%ecuaciones

yp(1)=y(2);

yp(2)=2*y(4)+y(1)-nu*(y(1)+mu)/d1-mu*(y(1)-mu)/d2;

yp(3)=y(4);

yp(4)=-2*y(2)+y(3)-nu*(y(3)+mu)/d1-mu*(y(3)-mu)/d2;

matriz=[];

for i=1:201 %puedes reemplzar el 201 por h

for j=1:4

matriz(i,j)=rand; %aquí donde va rand colocas la operación que vas a realizar 

end

end

a={@(x) x+2,@(x) x+1};

a{1}(2)    %aquí uso la primera función y le asigno un valor de 2 a la x

function [yp]=yprimas(y,t)

% y(1)= x1

% y(2)= x1´

% y(3)= x2

% y(4)= x2´

y=[2 3 4 5]

mu=1/82.45;

nu=1-mu;

yp=y;

h=201

% distancias del satélite a la tierra y a la luna,

d1=((y(1)+mu)^2+y(3)^2)^1.5;

d2=((y(1)-nu)^2+y(3)^2)^1.5;

%ecuaciones

yp(1)=y(2);

yp(2)=2*y(4)+y(1)-nu*(y(1)+mu)/d1-mu*(y(1)-mu)/d2;

yp(3)=y(4);

yp(4)=-2*y(2)+y(3)-nu*(y(3)+mu)/d1-mu*(y(3)-mu)/d2;

T=zeros(h,4)

for i=1:201

  for j=1:4

    T(i,j)=yp(i)----aqui es que no se como llamar las 4 "yp( )" que tengo mas arriba

  end

end

end