neural_network_function.matlab

x1=-1.0; x2=+1.0; y1=-1.5; y2=+4.5; z1=-1.0; z2=+1.0;
nx=7; ny=9; nz=5;
step_x=(x2-x1)/(nx-1); step_y=(y2-y1)/(ny-1); step_z=(z2-z1)/(nz-1);
step_min = min([step_x step_y step_z]); 
[x,y,z]=meshgrid([x1:step_x:x2], [y1:step_y:y2], [z1:step_z:z2]); 
f= 4.*x+x.^2.*y-z.*y;

isosurface(x,y,z,f), title('PS. Press<enter>');
pause;
close;

xx=reshape(x,1,nx*ny*nz); 
yy=reshape(y,1,nx*ny*nz); 
zz=reshape(z,1,nx*ny*nz);

ff= 4.*xx+xx.^2.*yy-zz.*yy;
p=[xx; yy; zz];

goal = 0.001;
spread = 1.0*step_min;

net = newrb(p,ff,goal,spread);
pause;
close;

net.layers{1}.size
smlt=sim(net,p); 
[ff' smlt']

smltr=reshape(smlt,ny,nx,nz);
isosurface(x,y,z,smltr), title('AS. Press<enter>');
pause;
close;