>> x=[179 172 183 165 193 186 169 181 176 177]

x =

   179   172   183   165   193   186   169   181   176   177

>> y=[180 175 183 170 189 187 174 183 180 179]

y =

   180   175   183   170   189   187   174   183   180   179

>> c=cov(x,y)

c =

   68.3222   47.7778
   47.7778   34.4444

>> varx=c(1,1)

varx =

   68.3222

>> vary=c(2,2)

vary =

   34.4444

>> maxx=max(x)

maxx =

   193

>> minx=min(x)

minx =

   165

>> maxy=max(y)

maxy =

   189

>> miny=min(y)

miny =

   170

>> prx=mean(x)

prx =

  178.1000

>> pry=mean(y)

pry =

   180

>> k1=c(1,2)/c(1,1)

k1 =

    0.6993

>> q1=pry-k1*prx

q1 =

   55.4545

>> plot(x,y,'.g'),hold on
>> plot([minx,maxx],[k1*minx+q1,k1*maxx+q1]),hold on, shg
>> plot(x,k1*x+q1,'+r'),shg

>> %dalsia priamka

>> k2=c(2,1)/c(2,2)

k2 =

    1.3871

>> q2=prx-k2*pry

q2 =

  -71.5774

>> plot([k2*miny+q2,k2*maxy+q2],[miny,maxy]),shg
>> plot(k2*y+q2,y,'k<'),shg