[X_ucz, Y_ucz] = meshgrid(-20:0.2:20,-20:0.2:20);
punkt_uczenie = [X_ucz(:)'; Y_ucz(:)'];

[X_test, Y_test] = meshgrid(-20:0.2:20,-20:0.2:20);
punkt_test = [X_test(:)'; Y_test(:)'];

funkcja_uczenie = 0.5+((sin(sqrt(X_ucz.^2+Y_ucz.^2)-0.5)).^2)./((1+0.001*(X_ucz.^2+Y_ucz.^2)).^2);
funkcja_test = 0.5+((sin(sqrt(X_test.^2+Y_test.^2)-0.5)).^2)./((1+0.001*(X_test.^2+Y_test.^2)).^2);

figure(1);
mesh(X_test,Y_test,funkcja_test);
title('Wykres funkcji aproksymowanej');
axis([-20 20 -20 20 0 2]);

net = newff([-20 20; -20 20],[10 10 1],{'tansig' 'tansig' 'purelin'},'trainlm');
net.trainParam.epochs = 100;
net.trainParam.goal = 1e-10;

net = init(net);
[net tr] = train(net,punkt_uczenie,funkcja_uczenie(:)');

Z = sim(net,punkt_test);
Z = reshape(Z,size(X_test));

figure(3);
mesh(X_test,Y_test,Z);
title('Wykres funkcji aproksymującej');
axis([-20 20 -20 20 0 2]);