tipo=input('ingrese el tipo de aproximacion (continua (C) o discreta (D))=','s' ); 
switch tipo
case 'D'
    x=input('ingrese las coordenadas x ');
    y=input('ingrese las coordenadas y ');
    m=length(x);
    n=input('ingrese el grado del polinomio por el cual quiere aproximar ');
    A=zeros(n+1,m);
    for j=0:n;
        A(j+1,:)=(x.^(j));
    end
    c=(A*A')\(A*y');
    disp('los coeficientes son [an ... a0]=');P=fliplr(c')
    cond=condest(A*A')%estima el numero de condicion de la matriz
    error=norm(y'-A'*c)%error cuadratico medio
    v=x(1):0.0001:x(m);
    w=polyval(P,v);
    plot(x,y,'o',v,w,'r');
case 'C'



end