x=0:0.01:1;
y=4*x./(1+10*(x.^2));
x=x';
y=y';
plot(x,y,'blue');
hold on;
n=4; %grado de los polinomios
V=zeros(length(x),n+1);
for k=1:n+1
    V(:,k)=(x.^(k-1));
end
thetals= V\y; %coeficientes del polinomio de cuadrados minimos
%matlab calcula cuadrados mínimos cuando no hay solucion, es decir
%inv(V'*V)*V'*y pero de formas mas eficientes.
yls=V*thetals; %aproximacion por least squares de la función.
plot(x,yls,'red');
hold off;