I have found the response to my question, we can find the right solution by mapping and reverse mapping
entree=xlsread('O:\Stage\database.xlsx',1);
sortie=xlsread('O:\Stage\database.xlsx',2);
entree=entree';
sortie=sortie';
N=length(sortie);
net=fitnet(n);
net.numLayers;
[net,TR]=train(net,entree,sortie,'UseParallel','yes');
%%%%%%%% do and extract Mapping settings %%%%
[entree_map,entree_setMap] = mapminmax (entree);
[sortie_map,sortie_setMap] = mapminmax (sortie);
%%% apply mapping setting + regression Modele test + mapping reverse %%%%%%%%%%%%%%%%
e=[1155 10 0.7 1.00017];
emap= mapminmax ('apply',e',entree_setMap);
b1 = net.b {1};
b2 = net.b {2};
IW = net.IW {1,1};
LW = net.LW {2,1};
B1=repmat(b1,1,N);
B2=repmat(b2,1,N);
a1=IW*emap+b1
h1=tansig(a1);
a2=LW*h1+b2;
h3_map=purelin(a2);
h= mapminmax('reverse',h3_map,sortie_setMap)
y=net(e')
, but it is a nother way to do it, and I think is the best: