%% Script file to test the convertNet program

X=rand(10,2);  %10 points in R^2
Y=rand(10,3);  %10 points in R^3

net=newrb(X',Y');
[C,W,b,spr]=convertNet(net);

Z1=sim(net,X');  %These are the outputs from the network.

Z1'

%Manually construct the RBF:
A=edm(X,C);
Phi=rbf1(A,1,spr);
Z2=Phi*W+repmat(b,10,1);

Z2