Nojrony2

function ans = neutrony2(x)
%     X = [1 2 1 2; 1 1 2 2];
%     Y = [0 1 1 0 ];
%     net = newff([0 3; 0 3], [2 1], {'logsig' 'logsig'},'trainlm');
%     net.trainParam.goal = 0.00001;
%     net.trainParam.epochs = 200;
%     net.trainParam.show = 1;
%     net = init(net);
%     net = train(net,X,Y);
%     [X1 Y1] = meshgrid(0:0.1:3);
%     Z = Y1;
%     Z(:) = sim(net, [X1(:)';Y1(:)']);
%     surf(X1,Y1,Z);
%     X = -1:0.05:1;
%     Y = sin(3*pi*X).^(2).*sin(pi*X);
%     XW = -1:0.01:1;
%     YW = sin(3*pi*XW).^(2).*sin(pi*XW);
%     net = newff([-1 1],[10 1],{'tansig' 'purelin'},'trainlm');
%     net.trainParam.goal = 0.0000001;
%     net.trainParam.epochs = 50000;
%     %net.trainParam.show = 100;
%     net = init(net);
%     net = train(net,X,Y);
%     YS = sim(net,XW);
%     plot(XW,YW,'r',X,Y,'gO',XW,YS,'b');

% a b c h
% f = V1 - V2
% f = (a^3) - ((b*c)*((1/3)*h))

%Dane łóczonce
X = rand(4,500);
Y = zeros(1,500);
for i=1:500
    Y(i) = (X(1,i)^3) - ((X(2,i)*X(3,i))*((1/3)*X(4,i)));
end

% Źeć nełronowa
net = newff([0 1; 0 1; 0 1; 0 1], [60 1], {'tansig' 'tansig'},'trainlm');
net.trainParam.goal = 1e-5;
net.trainParam.epochs = 800;
net.trainParam.show = 50;
net = init(net);
net = train(net,X,Y);

%Cpiur dezdoffy
XT = rand(4,50);
for i=1:50
    YT(i) = (XT(1,i)^3) - ((XT(2,i)*XT(3,i))*((1/3)*XT(4,i)));
end

%Testoffanie Źeći
YS = sim(net,X);
plot(1:500, abs(YS-Y), 'b', 1:500, genius(1:100,1e-5), 'g');
YSS = zeros(1,500);
for i=1:500
   YSS(i) = abs(Y(i) - YS(i));
end
ans = sum(YSS)/500;
end