RosettaCodeData/Task/Cumulative-standard-deviation/PARI-GP/cumulative-standard-deviation.parigp
2023-07-01 13:44:08 -04:00

19 lines
328 B
Text

newpoint(x)={
myT=x;
myS=0;
myN=1;
[myT,myS]/myN
};
addpoint(x)={
myT+=x;
myN++;
myS+=(myN*x-myT)^2/myN/(myN-1);
[myT,myS]/myN
};
addpoints(v)={
print(newpoint(v[1]));
for(i=2,#v,print(addpoint(v[i])));
print("Mean: ",myT/myN);
print("Standard deviation: ",sqrt(myS/myN))
};
addpoints([2,4,4,4,5,5,7,9])