port(:,1)=[5*ones(50,1);15*ones(10,1);30*ones(20,1)];
port(:,2)=[10*ones(40,1);0*ones(10,1);20*ones(30,1)];
port(:,3)=[0*ones(70,1);80*ones(5,1);120*ones(5,1)];
port(:,4)=[10*ones(30,1);20*ones(35,1);0*ones(15,1)];
port(:,5)=[20*ones(50,1);0*ones(30,1)];
base=calculabase2(port);
target=[0:20:1580]';
aprox=base*base'*target;
plot(target,'o',"markersize", 20,aprox,'xr',"markersize", 20)