Matlab - regresión lineal inicie en cero

x=[1 2 3 4 5 6 7 8]; %ejemplo

y=[4 9 11 15 18 26 34 45]; %ejemplo

plot(x-x(1),y-y(1))

function capas

D1= input('ingrese vector de Distancia de ida D1= ');

T1= input ('ingrese vector de Tiempo de ida T1= ');

subplot(1,2,1)

plot(D1,T1,'.','Color','r')

grid on

hold on

xlabel('Distancia (m)')

ylabel('Tiempo (ms)')

title('Dromocronas datos de ida')

c=input ('ingrese numero de rectas para tiempo de ida c= ');

for k=1:+1:c

    l=input ('ingrese vector con inicio y fin de recta k l= ');

    P=l(1,2)-l(1,1)+1;

    r=l(1,2);

    s=l(1,1);

    S1=0;

    S2=0;

    S3=0;

    S4=0;

    for i=s:+1:r

        n=P;

        B1=D1(1,i);

        S1=S1+B1;

        B2=T1(1,i);

        S2=S2+B2;

        B3=D1(1,i)*T1(1,i);

        S3=S3+B3;

        B4=(D1(1,i))^2;

        S4=S4+B4;

    end

    m(1,k)=((n*S3)-(S1*S2))/((n*S4)-(S1^2))

    b(1,k)=(S2/P)-(m(1,k)*(S1/P))

    V(1,k)=(1/m(1,k))*1000

    X1=0:+1:D1(1,r)+10;

    Y1=b(1,k)+m(1,k)*X1;

    subplot(1,2,1)

    plot (X1,Y1,'r')

    grid ('on')

    hold on

end

for o=1:+1:c-1

    X(1,o)=abs(b(1,o+1)-b(1,o))/abs(m(1,o+1)-m(1,o));

end