#
# canola - canon canola 1614p emulator
# Copyright (C) 2011 Peter Miller
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License, version 3, as
# published by the Free Software Foundation.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
# General Public License for more details.
#
# You should have received a copy of the GNU General Public License along
# with this program. If not, see <http://www.gnu.org/licenses/>.
#
; ------------------  Card-Pocket-Begin  -----------------------------------
; Title: ST-4: Summations of Paired Data
; Color: green
; Style: blank
; ------------------  Card-Pocket-End  -------------------------------------
;
; How-to-Operate:
;     [START]
;  x1 [START] y1 [START]
;  x2 [START] y2 [START]
;  x3 [START] y3 [START]
;   :    :     :    :
;  xn [START] yn [START]
;     [START] => n
;                sigma(x)
;                sigma(x**2)
;                sigma(y)
;                sigma(x*y)
;                sigma(y**2)
;
; Example:
;      [START]
;  5.3 [START] 6.7 [START]
;  3.5 [START] 4.9 [START]
;  8.3 [START] 6.1 [START]
;  7.3 [START] 8.2 [START]
;      [START] =>  4      ; n
;                 24.4    ; sigma(x)
;                162.52   ; sigma(x**2)
;                 25.9    ; sigma(y)
;                163.15   ; sigma(x*y)
;                173.35   ; sigma(y**2)
;
    CM14        ; sigma(x)
    CM13        ; sigma(x**2)
    CM12        ; sigma(y)
    CM11        ; sigma(x*y)
    CM10        ; sigma(y**2)
    CM9         ; n
    ENT
FJ $614
    M14
    *
    +=
    M13
    RV
    *
    ENT
    SM8
    M12
    +=
    M11
    RM8
    *
    +=
    M10
    1
    M9
    RM9
    EJ $614
    PRINT       ; n
    RM14
    PRINT       ; sigma(x)
    RM13
    PRINT       ; sigma(x**2)
    RM12
    PRINT       ; sigma(y)
    RM11
    PRINT       ; sigma(x*y)
    RM10
    PRINT       ; sigma(y**2)