function XC=montecarloerrpLBUFRR(n) %This function takes n random guesses with in a logical interval %This was the starting point for fitting, but linfit was used in the final %version, which includes a linear optimisation setp from each montecarlo %guess k=1; XC=[1,1,1,1,1,1]; M=[(1/25)*0.001.*rand(n,1),(1/25)*0.001.*rand(n,1),4000*rand(n,1),200*rand(n,1), 100*rand(n,1)]; for k=1:n X0(1,1)=M(k,1); X0(1,2)=M(k,2); X0(1,3)=M(k,3); X0(1,4)=M(k,4); X0(1,5)=M(k,5); ERR0(k)=toterrorpLBUFRR(X0(1,1),X0(1,2),X0(1,3),X0(1,4),X0(1,5)); rel(k)=X0(1,4)/X0(1,5); k=k+1; end [mindif i]=min(ERR0); XC(1,1)=M(i,1); XC(1,2)=M(i,2); XC(1,3)=M(i,3); XC(1,4)=M(i,4); XC(1,5)=M(i,5); XC(1,6)=mindif;