Matlab - Crear una matriz

k=input('ingrese la cantidad de puntos a generar: ')

for i = 1 : k

    x(i) = (i-1)/10;

    D(i,1)=x(i);

endfor

for j = 1 : k

    y(j)=erf(x(j));

    D(j,2)=y(j);

endfor

disp('Los datos generados son:')

disp('      x        y')

disp(D)

%La matriz que contiene los datos se llama D

%Fila uno corresponde a las abscisas

%Fila dos corresponde a las ordenadas al origen

n=input('Ingrese el grado impar del polinomio: ');

k=input('Ingrese la cantidad de puntos a generar: ');

  for i = 1 : k                                                     %Genera los datos que usamos para hacer el ajuste por cuadrado mínimos 

    x(i) = (i-1)/10;

    y(i)=erf(x(i));

      for w = 1 : (n+1)/2                                           %Genera la matriz X para el calculo de las ecuaciones normales

          X(i,w+1)=(x(i))^(2*w-1);

          X(1:i,1)=1;

      endfor

   endfor

X;

y=y';

v=size(X);

  for j = 2 : v(2)                               %Genera submatriz a partiir de la matrz X para calcular los coeficientes de todos los polinomios

    x=X(1:k,1:j)';

    A=x*x';

    b=x*y;

    C(j,1:j) = fliplr((A\b)');                   %Guarda los coeficientes en una matriz C

    C(1,j) = j-1;

  endfor

P1(x)=c0*x + c1

P3(x)=c0*x^3 +c1*x + c2

P5(x)=co*x^5+c1*x^3+c2*x+c3

P7(x)=co*x^7+c1*x^5+c2*x^3+c3*x+c4

P9(x)=co*x^9+c1*x^7+c2*x^5+c3*x^3+c4*x+c5

for i = 1 : k %Genera los datos que usamos para hacer el ajuste por cuadrado mínimos 

x(i) = (i-1)/10;

y(i)=erf(x(i));

endfor