x:[1,2,3,4,5];

h:[1,2,10,h4,1];

s:read("Gib die Standardabweichung ein");

v:s**2;

n:length(x);

N:sum(h[i],i,1,n);

p:h/N;

m:sum(p[i]*x[i],i,1,n);

g:v=sum(p[i]*(x[i]-m)**2,i,1,n);

l:solve(g,h4);