function xy = schwerpt(f, a, b)
A = integral(f, a, b);
xs = integral(@(t) t .* f(t), a, b) / A;
ys = integral(@(t) f(t).^2, a, b) / A / 2;
xy = [xs; ys];
end
