 function [Profile,ParamForSave]=PR_func(counter,imu,isig)% Add function in standart way
dist=2^(-16);
x=0:dist:1-dist;
load(['Rand/gaus_rand_N' num2str(counter) '.mat'])
len=length(x);

Profile=0;%Profile shape
ParamForSave=[];
for iloop=1:5

Amplitude=abs(Amplitude_rand(iloop));
Sigma=abs(Sigma_rand(iloop)+0.5)*0.000000006*isig;%explain
mu=(mean(x)/len+(Mu_rand(iloop)-0.5)*dist*0.0015*imu); %generate array of random value - our coefficients
 
ParamForSave=[ParamForSave,Amplitude,mu,Sigma];
Profile=Profile+Amplitude.*exp(-((x/len-mu).^2./(2.*((Sigma)^2))));

end
Profile=Profile/sum(Profile);



