samer@16: % quadf - scalar quadratic form
samer@16: %
samer@16: % quadf :: [Size], real -> [Size].
samer@16: % quadf :: [Size] -> [Size].
samer@16: %
samer@16: % Scalar function: quadf(x,s) = -x^2/(2*s), s defaults to 1.
samer@16: function y=quadf(x,s), if nargin<2, y=-(x.^2)/2; else y= -(x.^2)/(2*s); end