function [xbar, s] = mitstandfor(x) n=length(x); xbar=0; s=0; for i=1:n xbar=xbar+x(i); end xbar=xbar/n; if n==1 disp('Standardabweichung nicht definiert') else for i=1:n s=s+(x(i)-xbar)^2; end s=s/(n-1); s=sqrt(s); end