[u1g u2g] = meshgrid(linspace(0,1,10), linspace(0,1,10));
u1 = u1g(:);
u2 = u2g(:);
input = [u1 u2];
[N n] = size(input);
c = zeros(1,n);
sigma = 0.25;
hheight = 1;
gheight = 0.75;
den = 0.1;
for i = 1:n
den = den+(1-input(:,i))./n;
end
hyperbel = 0.1./den;
x = zeros(N,1);
for i = 1:2
x = x + ((input(:,i)-c(i))/sigma).^2;
end
gauss = exp(-0.5.*x);
output = hheight.*hyperbel + (gheight-hyperbel).*gauss;
output = output + 0.01*randn(size(output))*(max(output)-min(output));
[LMNBest AllLMN] = LMNTrain([input output]);
Current training method: lolimotQuad
Current training method: lolimotSparseQuad
Current training method: lolimot
Current training method: hilomotSparseQuad
Current training method: hilomotQuad
Current training method: hilomot
Suggested model: lolimotQuad