[u1g u2g] = meshgrid(linspace(0,1,10), linspace(0,1,10));
u1 = u1g(:);
u2 = u2g(:);
input = [u1 u2];
u1f = u1*15-5;
u2f = u2*15;
output = (u2f-(5.1/(4*pi^2))*u1f.^2 + 5*u1f/pi-6).^2 + 10*(1-1/(8*pi))*cos(u1f)+10;
output = output + 0.1*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: lolimot