| > | restart; |
| > | A:=Matrix([[25,18,10,5.],[20,16,15,6],[15,20,22,13],[8,10,7,15]]); |
| > |
| > | Sm:=add(add(A[i,j],i=1..4),j=1..4); |
| > |
| > | for j to 4 do Ay[j]:=add(A[i,j],i=1..4); end; |
| > | for j to 4 do Ax[j]:=add(A[j,i],i=1..4); end; |
![Ay[1] := 51](images/tv143.gif){kind=small}
![Ay[2] := 52](images/tv144.gif){kind=small}
![Ay[3] := 75](images/tv145.gif){kind=small}
![Ay[4] := 86.](images/tv146.gif){kind=small}
![Ax[1] := 15.](images/tv147.gif){kind=small}
![Ax[2] := 56](images/tv148.gif){kind=small}
![Ax[3] := 69](images/tv149.gif){kind=small}
![Ax[4] := 124](images/tv1410.gif){kind=small}
| > | for i to 4 do |
| > | for j to 4 do B[i,j]:=Ax[i]*Ay[j]/Sm; end:end: |
| > | Sm:=add(add((A[i,j]-B[i,j])^2/B[i,j],i=1..4),j=1..4); |
| > |
| > |
| > |