wks.addCol();
wcol(wks.ncols)[L]$ = Peak1;  // wks.ncols gives the index of last column
wcol(wks.ncols)[U]$ = Area;
wcol(wks.ncols)[O]$= "ave(B,5)";

wks.addCol();
wcol(wks.ncols)[L]$ = Peak1;
wcol(wks.ncols)[U]$ = Error;
wcol(wks.ncols)[O]$= "sqrt((ave(B^2,5) - ave(B,5)^2))"; // std = sqrt((E(X^2) - E(X)^2))

import numpy as np
import originpro as op
 
wks = op.find_sheet()
nn = wks.cols

wks1 = op.new_sheet()
wks1.cols = 0

for i, idx in enumerate(range(1, nn, 9)):  # loop over second col of every 9 columns
    data = np.asarray(wks.to_list(idx))
    nc = len(data) // 5
    data = np.split(data[:5*nc],nc)
    wks1.from_list(wks1.cols, np.mean(data,1), f"Peak{i+1}", 'Area')
    wks1.from_list(wks1.cols, np.std(data,1), f"Peak{i+1}", 'Error')