function [b]= computeI() %%%%%%%%%%%%%%%%%%%%%%%%%%%% d01 = load('./test-0100'); d02 = load('./test-0090'); d03 = load('./test-0080'); d04 = load('./test-0070'); d05 = load('./test-0060'); d06 = load('./test-0050'); d07 = load('./test-0040'); d08 = load('./test-0030'); d09 = load('./test-0020'); d10 = load('./test-0010'); d11 = load('./test-0000'); x = [d01(:,2);d02(:,2);d03(:,2);d04(:,2);d05(:,2);d06(:,2);d07(:,2);d08(:,2);d09(:,2);d10(:,2);d11(:,2)]; y = [ones(length(d01),1)]; y2 = [y*100;y*90;y*80;y*70; y*60;y*50;y*40;y*30;y*20;y*10;y*0]; A = [x ones(length(x),1)]; b = A\y2;