%fitting parameters
a=1.016;
b=-7.943;
c=0.02556;
n=1.07;
p=1.2343;
f=0.9225;
parm0=[a b c n p f];

%equation
%K=a+b(c*vp.^1./vin.^n*(0.75*10^(-9))^(n-1)*1.46^(p-n+1)/0.00075^(p))^f

options = optimset('Display','iter','TolFun',10^(-10),'TolX',10^(-10),'MaxIter',10000,'MaxFunEvals',10000);
fitfun = @(x) ErrorEstimation(x);
[parm1,fval,exitflag,output]=fminsearch(fitfun,parm0,options);

a=parm1(1);
b=parm1(2);
c=parm1(3);
n=parm1(4);
p=parm1(5);
f=parm1(6);

%% Plot
%vin
vin=zeros(6,15);
vp=zeros(6,15);
Phi=zeros(6,15);
vin(1,:)=0.1;
vin(2,:)=0.2;
vin(3,:)=0.3;
vin(4,:)=0.4;
vin(5,:)=0.5;
vin(6,:)=0.6;
vp(1,:)=[3.12234E-07	7.31008E-07	1.17054E-06	1.61585E-06	2.06239E-06	2.50831E-06	2.95258E-06	3.39472E-06	3.83396E-06	4.27029E-06	4.70285E-06	5.13133E-06	5.55574E-06	5.97521E-06	6.38916E-06];
vp(2,:)=[3.21914E-07	6.91474E-07	1.08E-06	1.46943E-06	1.86364E-06	2.25875E-06	2.65404E-06	3.04909E-06	3.44356E-06	3.83774E-06	4.23076E-06	4.62291E-06	5.01418E-06	5.40458E-06	5.79353E-06];
vp(3,:)=[3.68E-07	7.42345E-07	1.13106E-06	1.52524E-06	1.92181E-06	2.31939E-06	2.71735E-06	3.11537E-06	3.51304E-06	3.91071E-06	4.3075E-06	4.70401E-06	5.09993E-06	5.49528E-06	5.88975E-06];
vp(4,:)=[4.15634E-07	7.93158E-07	1.18379E-06	1.57954E-06	1.97753E-06	2.3766E-06	2.77613E-06	3.17583E-06	3.57554E-06	3.97466E-06	4.37378E-06	4.77261E-06	5.17086E-06	5.56882E-06	5.9662E-06];
vp(5,:)=[4.6357E-07	8.4403E-07	1.23623E-06	1.63312E-06	2.03218E-06	2.43227E-06	2.83287E-06	3.23368E-06	3.63455E-06	4.03512E-06	4.4357E-06	4.83569E-06	5.23569E-06	5.6351E-06	6.03423E-06];
vp(6,:)=[5.12057E-07	8.95047E-07	1.28856E-06	1.68637E-06	2.08622E-06	2.48712E-06	2.88854E-06	3.29007E-06	3.69181E-06	4.09355E-06	4.49471E-06	4.89616E-06	5.29703E-06	5.6976E-06	6.09789E-06];
Phi(1,:)=[1.095140173	1.055699536	1.032037149	1.015753886	1.003934898	0.995082525	0.98822036	0.98269316	0.978037998	0.974092731	0.970580126	0.967409684	0.964531826	0.961846349	0.959312757];
Phi(2,:)=[1.055591544	1.048904013	1.03724734	1.027212151	1.019072605	1.012451139	1.007263689	1.003061992	0.99960162	0.996945534	0.994667808	0.992754629	0.99118528	0.989868892	0.988667791];
Phi(3,:)=[1.032968305	1.040204959	1.034118857	1.027131587	1.02096664	1.015554349	1.011100446	1.007608062	1.004601255	1.002306909	1.000373087	0.998791788	0.997523486	0.996518631	0.995642928];
Phi(4,:)=[1.016843865	1.032713459	1.030377655	1.025710446	1.02096876	1.016395694	1.012715658	1.009551127	1.006948715	1.004738428	1.002944106	1.001588985	1.000421025	0.999526262	0.998770972];
Phi(5,:)=[1.006906314	1.026107903	1.026935916	1.023895877	1.020070993	1.016464628	1.013126699	1.010298886	1.007978038	1.0059482	1.004392317	1.003011346	1.002002435	1.001126504	1.000417647];
Phi(6,:)=[0.997531006	1.020125319	1.023208667	1.02197346	1.018916328	1.015889123	1.013138462	1.01055906	1.008449864	1.006674429	1.00508788	1.003895976	1.002833779	1.002038753	1.001410777];
Phi=Phi(:,3:end);
vp=vp(:,3:end);
vin=vin(:,3:end);

figure;
plot(vp.^1./vin.^n*(0.75*10^(-9))^(n-1)*1.46^(p-n+1)/0.00075^(p),Phi,'.')
hold on;

%d
d=zeros(9,15);
vp=zeros(9,15);
Phi=zeros(9,15); 
d(1,:)=0.00075;
d(2,:)=0.00085;
d(3,:)=0.001;
d(4,:)=0.00125;
d(5,:)=0.0015;
d(6,:)=0.00175;
d(7,:)=0.002;
d(8,:)=0.004;
d(9,:)=0.008;
vp(1,:)=[3.12147E-07	7.30863E-07	1.17033E-06	1.61556E-06	2.06204E-06	2.5079E-06	2.95229E-06	3.39414E-06	3.83338E-06	4.26942E-06	4.70197E-06	5.13046E-06	5.55458E-06	5.97405E-06	6.388E-06];
vp(2,:)=[3.010734E-07	7.188776E-07	1.158176E-06	1.603399E-06	2.049931E-06	2.495847E-06	2.940198E-06	3.382395E-06	3.821771E-06	4.258325E-06	4.691288E-06	5.120660E-06	5.545928E-06	5.966579E-06	6.382613E-06];
vp(3,:)=[2.90272E-07	7.07E-07	1.14572E-06	1.59059E-06	2.04E-06	2.48238E-06	2.92627E-06	3.37E-06	3.80729E-06	4.24355E-06	4.68E-06	5.1056E-06	5.53118E-06	5.95E-06	6.36925E-06];
vp(4,:)=[2.801999E-07	6.955027E-07	1.133532E-06	1.577648E-06	2.023056E-06	2.467643E-06	2.910660E-06	3.351411E-06	3.789370E-06	4.224365E-06	4.656045E-06	5.084237E-06	5.508767E-06	5.929110E-06	6.345268E-06];
vp(5,:)=[2.746331E-07	6.889302E-07	1.126134E-06	1.569311E-06	2.013637E-06	2.457091E-06	2.898800E-06	3.338184E-06	3.774661E-06	4.208086E-06	4.638167E-06	5.064616E-06	5.487139E-06	5.905738E-06	6.320123E-06];
vp(6,:)=[2.71317E-07	6.847462E-07	1.121074E-06	1.563269E-06	2.006411E-06	2.448681E-06	2.888958E-06	3.326619E-06	3.761414E-06	4.192970E-06	4.621037E-06	5.045367E-06	5.465835E-06	5.882068E-06	6.294189E-06];
vp(7,:)=[2.69059E-07	6.817391E-07	1.117246E-06	1.558410E-06	2.000447E-06	2.441393E-06	2.880159E-06	3.316309E-06	3.749407E-06	4.179016E-06	4.605028E-06	5.027225E-06	5.445388E-06	5.859300E-06	6.268852E-06];
vp(8,:)=[2.634069E-07	6.726477E-07	1.103347E-06	1.538625E-06	1.973521E-06	2.406237E-06	2.835683E-06	3.261314E-06	3.682693E-06	4.099602E-06	4.511770E-06	4.918978E-06	5.321226E-06	5.71813E-06	6.11002E-06];
vp(9,:)=[2.614147E-07	6.665703E-07	1.090811E-06	1.517695E-06	1.942508E-06	2.363315E-06	2.779216E-06	3.189367E-06	3.593795E-06	3.992227E-06	4.384119E-06	4.769743E-06	5.149098E-06	5.521913E-06	5.888733E-06];
Phi(1,:)=[1.095140173	1.055700E+00	1.032037E+00	1.015754E+00	1.003935E+00	9.950825E-01	9.882204E-01	9.826932E-01	9.780380E-01	9.740927E-01	9.705801E-01	9.674097E-01	9.645318E-01	9.618463E-01	9.593128E-01];
Phi(2,:)=[5.539450E-01	1.056425E+00	1.032349E+00	1.016474E+00	1.005435E+00	9.974272E-01	9.913802E-01	9.866515E-01	9.828253E-01	9.796641E-01	9.769163E-01	9.744982E-01	0.972298829	0.970274968	0.968401913];
Phi(3,:)=[1.102908174	1.056333E+00	1.032170E+00	1.016886E+00	1.006752E+00	9.996870E-01	9.945541E-01	9.907976E-01	9.878425E-01	9.854903E-01	9.834967E-01	9.817860E-01	9.802767E-01	9.788802E-01	9.775939E-01];
Phi(4,:)=[1.104946E+00	1.054815E+00	1.030861E+00	1.016544E+00	1.007504E+00	1.001553E+00	9.975767E-01	9.948158E-01	9.928193E-01	9.913146E-01	9.901179E-01	9.891128E-01	9.882412E-01	9.874220E-01	9.866561E-01];
Phi(5,:)=[1.065936E+00	1.045374E+00	1.025142E+00	1.013327E+00	1.005468E+00	1.000860E+00	9.979917E-01	9.961581E-01	9.949503E-01	9.942173E-01	9.935876E-01	9.930632E-01	9.925565E-01	9.920685E-01	0.991620443];
Phi(6,:)=[1.10325287	1.05111062	1.027710778	1.014724731	1.007125636	1.002649781	0.999966359	0.998356329	0.997388687	0.996781804	0.996369283	0.995997682	0.995694245	0.99534514	0.995004933];
Phi(7,:)=[1.102144E+00	1.049449E+00	1.026363E+00	1.013818E+00	1.006732E+00	1.002711E+00	1.000443E+00	9.992161E-01	9.985393E-01	9.982133E-01	9.979520E-01	9.977707E-01	9.975886E-01	9.973628E-01	9.970968E-01];
Phi(8,:)=[1.093472845	1.04021298	1.018680707	1.008806997	1.004141732	1.00222237	1.001658105	1.001758844	1.002089602	1.002441487	1.002719658	1.002867087	1.002897356	1.00280074	1.002633801];
Phi(9,:)=[1.084619956	1.030996616	1.01156887	1.003998715	1.001438071	1.00106761	1.001521795	1.002157306	1.002730537	1.003137136	1.003302875	1.003297407	1.003159499	1.002896386	1.002576827];
Phi=Phi(:,3:end);
vp=vp(:,3:end);
d=d(:,3:end);
plot(vp.^1./0.1.^n*(0.75*10^(-9))^(n-1)*1.46^(p-n+1)./d.^(p),Phi,'.')

%D
D=zeros(5,15); 
vp=zeros(5,15); 
Phi=zeros(5,15); 
D(1,:)=0.55*10^(-9);
D(2,:)=0.75*10^(-9);
D(3,:)=0.85*10^(-9);
D(4,:)=0.95*10^(-9);
D(5,:)=1.15*10^(-9);
vp(1,:)=[3.11479E-07	7.27607E-07	1.16315E-06	1.60347E-06	2.04407E-06	2.48311E-06	2.91944E-06	3.35257E-06	3.78193E-06	4.20663E-06	4.62669E-06	5.04121E-06	5.45022E-06	5.85283E-06	6.24905E-06];
vp(2,:)=[3.12089E-07	7.30746E-07	1.17019E-06	1.61536E-06	2.06181E-06	2.50761E-06	2.95171E-06	3.39385E-06	3.83309E-06	4.26884E-06	4.70139E-06	5.12988E-06	5.554E-06	5.97318E-06	6.38713E-06];
vp(3,:)=[3.12234E-07	7.31735E-07	1.17243E-06	1.61922E-06	2.06756E-06	2.51558E-06	2.96217E-06	3.40722E-06	3.84966E-06	4.28919E-06	4.72581E-06	5.15895E-06	5.58801E-06	6.01271E-06	6.43277E-06];
vp(4,:)=[3.12467E-07	7.32781E-07	1.17464E-06	1.62289E-06	2.07291E-06	2.52285E-06	2.97176E-06	3.41914E-06	3.86419E-06	4.30721E-06	4.74703E-06	5.18395E-06	5.61708E-06	6.04643E-06	6.47143E-06];
vp(5,:)=[3.12729E-07	7.34177E-07	1.17778E-06	1.62812E-06	2.08061E-06	2.5334E-06	2.98543E-06	3.43629E-06	3.88571E-06	4.33308E-06	4.77785E-06	5.22028E-06	5.65981E-06	6.09585E-06	6.5284E-06];
Phi(1,:)=[0.471467635	1.046284743	1.023237255	1.008408251	0.998464114	0.99140006	0.986154067	0.982096925	0.978807543	0.975993119	0.973543503	0.971319105	0.969302824	0.967420334	0.96568654];
Phi(2,:)=[1.096709306	1.055922382	1.031798391	1.015520386	1.003726413	0.994886402	0.987921011	0.982415193	0.977756305	0.973738647	0.970249793	0.967104588	0.964206532	0.961512262	0.958992465];
Phi(3,:)=[1.100556722	1.059811764	1.03540187	1.018549196	1.006138233	0.996505893	0.988901896	0.982753065	0.977512845	0.972960714	0.968987284	0.965380342	0.962061111	0.958975439	0.956078273];
Phi(4,:)=[1.102821961	1.063429917	1.039008422	1.021619217	1.008518857	0.998413797	0.990191752	0.983426705	0.977646829	0.972670563	0.968194964	0.964154165	0.960433192	0.956956883	0.953682624];
Phi(5,:)=[1.108264085	1.069184776	1.044750361	1.026743064	1.012843017	1.001734605	0.992548866	0.984761554	0.978073469	0.972216281	0.966857765	0.96203801	0.957571132	0.953365315	0.949394332];
Phi=Phi(:,3:end);
vp=vp(:,3:end);
D=D(:,3:end);
plot(vp.^1./0.1.^n.*(D).^(n-1)*1.46^(p-n+1)./0.00075.^(p),Phi,'.')

%L
L=zeros(7,15); 
vp=zeros(7,15); 
Phi=zeros(7,15); 
L(1,:)=0.26;
L(2,:)=0.41;
L(3,:)=0.56;
L(4,:)=0.71;
L(5,:)=1.46;
L(6,:)=1;
L(7,:)=1.25;
vp(1,:)=[3.3385E-07	7.62197E-07	1.21317E-06	1.67186E-06	2.13366E-06	2.59692E-06	3.06084E-06	3.52475E-06	3.98867E-06	4.45242E-06	4.91585E-06	5.37862E-06	5.84107E-06	6.30286E-06	6.764E-06];
vp(2,:)=[3.31001E-07	7.57733E-07	1.20668E-06	1.66308E-06	2.12227E-06	2.58271E-06	3.04346E-06	3.5041E-06	3.96443E-06	4.42415E-06	4.88324E-06	5.3415E-06	5.79893E-06	6.25533E-06	6.7108E-06];
vp(3,:)=[3.28215E-07	7.53553E-07	1.20086E-06	1.65521E-06	2.11229E-06	2.57028E-06	3.02834E-06	3.48595E-06	3.94303E-06	4.39912E-06	4.8543E-06	5.30827E-06	5.76103E-06	6.21242E-06	6.66245E-06];
vp(4,:)=[3.25465134E-07	7.49537589E-07	1.19529098E-06	1.64797845E-06	2.10299720E-06	2.55873327E-06	3.01429001E-06	3.46918921E-06	3.92307221E-06	4.37581946E-06	4.82719184E-06	5.27706981E-06	5.72527404E-06	6.17192408E-06	6.61606351E-06];
vp(5,:)=[2.7631624E-07	6.3973020E-07	1.0211788E-06	1.4083250E-06	1.7972735E-06	2.1863965E-06	2.5749090E-06	2.9621715E-06	3.3482130E-06	3.7322198E-06	4.1141916E-06	4.4938379E-06	4.8708680E-06	5.2452818E-06	5.6164980E-06];
vp(6,:)=[3.2039E-07	7.42384E-07	1.18564E-06	1.63535E-06	2.08701E-06	2.53884E-06	2.98999E-06	3.43991E-06	3.88826E-06	4.33453E-06	4.77889E-06	5.22113E-06	5.66082E-06	6.09754E-06	6.53129E-06];
vp(7,:)=[3.15899E-07	7.36068E-07	1.17725E-06	1.62448E-06	2.07328E-06	2.52186E-06	2.96923E-06	3.415E-06	3.85843E-06	4.29914E-06	4.73713E-06	5.17207E-06	5.60361E-06	6.03108E-06	6.45482E-06];
Phi(1,:)=[1.03982183	1.048396281	1.041260834	1.03319755	1.026427841	1.020831657	1.01630133	1.012531596	1.009588756	1.007290087	1.005428078	1.003943057	1.002802438	1.001920384	1.001227768];
Phi(2,:)=[1.0538068E+00	1.0509658E+00	1.0398594E+00	1.0303427E+00	1.0225673E+00	1.0165070E+00	1.0117273E+00	1.0079518E+00	1.0050288E+00	1.0027209E+00	1.0009138E+00	9.9948929E-01	9.9836728E-01	9.9747891E-01	9.9676766E-01];
Phi(3,:)=[1.061209215	1.051816379	1.03837664	1.027607741	1.019300241	1.012842617	1.007850354	1.003967927	1.000931049	0.998622978	0.99665482	0.995178648	0.993970038	0.992945028	0.992109415];
Phi(4,:)=[1.0665528E+00	1.0520430E+00	1.0368965E+00	1.0251319E+00	1.0163688E+00	1.0094893E+00	1.0041779E+00	1.0002704E+00	9.9703400E-01	9.9450107E-01	9.9246119E-01	9.9084458E-01	9.8939519E-01	9.8823576E-01	0.987123314];
Phi(5,:)=[1.085017446	1.057476273	1.036175875	1.020511024	1.009004025	0.999984548	0.992958554	0.987185415	0.982430128	0.978351816	0.974777484	0.971609698	0.968719691	0.966068551	0.963594139];
Phi(6,:)=[1.088608324	1.05666306	1.036860859	1.022661213	1.012360058	1.004628348	9.99E-01	9.94E-01	9.90E-01	9.87E-01	9.85E-01	9.83E-01	9.81E-01	0.979033888	0.977423176];
Phi(7,:)=[1.093789549	1.056091319	1.034096473	1.01888539	1.007762026	0.99937799	0.992975365	0.987902343	0.98373267	0.980205082	0.977178566	0.974505784	0.972084962	0.969836604	0.96778241];
Phi=Phi(:,3:end);
vp=vp(:,3:end);
L=L(:,3:end);
plot(vp.^1./0.1.^n*(0.75*10^(-9))^(n-1).*L.^(p-n+1)./0.00075.^(p),Phi,'.')

%Correction Function
x=0:0.001:0.2;
plot(x,a+b*(c*x).^f,'--');

function Result = ErrorEstimation(x) %%variables x
a=x(1);
b=x(2);
c=x(3);
n=x(4);
p=x(5);
f=x(6);

Error=0;

%%vin,vp
vin=zeros(6,15); 
vp=zeros(6,15); 
Phi=zeros(6,15); 
vin(1,:)=0.1;
vin(2,:)=0.2;
vin(3,:)=0.3;
vin(4,:)=0.4;
vin(5,:)=0.5;
vin(6,:)=0.6;
vp(1,:)=[3.12234E-07	7.31008E-07	1.17054E-06	1.61585E-06	2.06239E-06	2.50831E-06	2.95258E-06	3.39472E-06	3.83396E-06	4.27029E-06	4.70285E-06	5.13133E-06	5.55574E-06	5.97521E-06	6.38916E-06];
vp(2,:)=[3.21914E-07	6.91474E-07	1.08E-06	1.46943E-06	1.86364E-06	2.25875E-06	2.65404E-06	3.04909E-06	3.44356E-06	3.83774E-06	4.23076E-06	4.62291E-06	5.01418E-06	5.40458E-06	5.79353E-06];
vp(3,:)=[3.68E-07	7.42345E-07	1.13106E-06	1.52524E-06	1.92181E-06	2.31939E-06	2.71735E-06	3.11537E-06	3.51304E-06	3.91071E-06	4.3075E-06	4.70401E-06	5.09993E-06	5.49528E-06	5.88975E-06];
vp(4,:)=[4.15634E-07	7.93158E-07	1.18379E-06	1.57954E-06	1.97753E-06	2.3766E-06	2.77613E-06	3.17583E-06	3.57554E-06	3.97466E-06	4.37378E-06	4.77261E-06	5.17086E-06	5.56882E-06	5.9662E-06];
vp(5,:)=[4.6357E-07	8.4403E-07	1.23623E-06	1.63312E-06	2.03218E-06	2.43227E-06	2.83287E-06	3.23368E-06	3.63455E-06	4.03512E-06	4.4357E-06	4.83569E-06	5.23569E-06	5.6351E-06	6.03423E-06];
vp(6,:)=[5.12057E-07	8.95047E-07	1.28856E-06	1.68637E-06	2.08622E-06	2.48712E-06	2.88854E-06	3.29007E-06	3.69181E-06	4.09355E-06	4.49471E-06	4.89616E-06	5.29703E-06	5.6976E-06	6.09789E-06];
Phi(1,:)=[1.095140173	1.055699536	1.032037149	1.015753886	1.003934898	0.995082525	0.98822036	0.98269316	0.978037998	0.974092731	0.970580126	0.967409684	0.964531826	0.961846349	0.959312757];
Phi(2,:)=[1.055591544	1.048904013	1.03724734	1.027212151	1.019072605	1.012451139	1.007263689	1.003061992	0.99960162	0.996945534	0.994667808	0.992754629	0.99118528	0.989868892	0.988667791];
Phi(3,:)=[1.032968305	1.040204959	1.034118857	1.027131587	1.02096664	1.015554349	1.011100446	1.007608062	1.004601255	1.002306909	1.000373087	0.998791788	0.997523486	0.996518631	0.995642928];
Phi(4,:)=[1.016843865	1.032713459	1.030377655	1.025710446	1.02096876	1.016395694	1.012715658	1.009551127	1.006948715	1.004738428	1.002944106	1.001588985	1.000421025	0.999526262	0.998770972];
Phi(5,:)=[1.006906314	1.026107903	1.026935916	1.023895877	1.020070993	1.016464628	1.013126699	1.010298886	1.007978038	1.0059482	1.004392317	1.003011346	1.002002435	1.001126504	1.000417647];
Phi(6,:)=[0.997531006	1.020125319	1.023208667	1.02197346	1.018916328	1.015889123	1.013138462	1.01055906	1.008449864	1.006674429	1.00508788	1.003895976	1.002833779	1.002038753	1.001410777];

Phi=Phi(:,3:end);
vp=vp(:,3:end);
vin=vin(:,3:end);

Recovery=((4*vp*1.46)./(0.1*0.00075)).^2;

PhiM = a+b*(c./vin.^n.*vp.^1*(0.75*10^(-9))^(n-1)*1.46^(p-n+1)/0.00075^(p)).^(f);
Error = Error + sum(sum(((Phi-PhiM)).^2.*Recovery,1),2); %error weight: recovery

%%d,vp
d=zeros(9,15); 
vp=zeros(9,15); 
Phi=zeros(9,15); 
d(1,:)=0.00075;
d(2,:)=0.00085;
d(3,:)=0.001;
d(4,:)=0.00125;
d(5,:)=0.0015;
d(6,:)=0.00175;
d(7,:)=0.002;
d(8,:)=0.004;
d(9,:)=0.008;
vp(1,:)=[3.12147E-07	7.30863E-07	1.17033E-06	1.61556E-06	2.06204E-06	2.5079E-06	2.95229E-06	3.39414E-06	3.83338E-06	4.26942E-06	4.70197E-06	5.13046E-06	5.55458E-06	5.97405E-06	6.388E-06];
vp(2,:)=[3.010734E-07	7.188776E-07	1.158176E-06	1.603399E-06	2.049931E-06	2.495847E-06	2.940198E-06	3.382395E-06	3.821771E-06	4.258325E-06	4.691288E-06	5.120660E-06	5.545928E-06	5.966579E-06	6.382613E-06];
vp(3,:)=[2.90272E-07	7.07E-07	1.14572E-06	1.59059E-06	2.04E-06	2.48238E-06	2.92627E-06	3.37E-06	3.80729E-06	4.24355E-06	4.68E-06	5.1056E-06	5.53118E-06	5.95E-06	6.36925E-06];
vp(4,:)=[2.801999E-07	6.955027E-07	1.133532E-06	1.577648E-06	2.023056E-06	2.467643E-06	2.910660E-06	3.351411E-06	3.789370E-06	4.224365E-06	4.656045E-06	5.084237E-06	5.508767E-06	5.929110E-06	6.345268E-06];
vp(5,:)=[2.746331E-07	6.889302E-07	1.126134E-06	1.569311E-06	2.013637E-06	2.457091E-06	2.898800E-06	3.338184E-06	3.774661E-06	4.208086E-06	4.638167E-06	5.064616E-06	5.487139E-06	5.905738E-06	6.320123E-06];
vp(6,:)=[2.71317E-07	6.847462E-07	1.121074E-06	1.563269E-06	2.006411E-06	2.448681E-06	2.888958E-06	3.326619E-06	3.761414E-06	4.192970E-06	4.621037E-06	5.045367E-06	5.465835E-06	5.882068E-06	6.294189E-06];
vp(7,:)=[2.69059E-07	6.817391E-07	1.117246E-06	1.558410E-06	2.000447E-06	2.441393E-06	2.880159E-06	3.316309E-06	3.749407E-06	4.179016E-06	4.605028E-06	5.027225E-06	5.445388E-06	5.859300E-06	6.268852E-06];
vp(8,:)=[2.634069E-07	6.726477E-07	1.103347E-06	1.538625E-06	1.973521E-06	2.406237E-06	2.835683E-06	3.261314E-06	3.682693E-06	4.099602E-06	4.511770E-06	4.918978E-06	5.321226E-06	5.71813E-06	6.11002E-06];
vp(9,:)=[2.614147E-07	6.665703E-07	1.090811E-06	1.517695E-06	1.942508E-06	2.363315E-06	2.779216E-06	3.189367E-06	3.593795E-06	3.992227E-06	4.384119E-06	4.769743E-06	5.149098E-06	5.521913E-06	5.888733E-06];
Phi(1,:)=[1.095140173	1.055700E+00	1.032037E+00	1.015754E+00	1.003935E+00	9.950825E-01	9.882204E-01	9.826932E-01	9.780380E-01	9.740927E-01	9.705801E-01	9.674097E-01	9.645318E-01	9.618463E-01	9.593128E-01];
Phi(2,:)=[5.539450E-01	1.056425E+00	1.032349E+00	1.016474E+00	1.005435E+00	9.974272E-01	9.913802E-01	9.866515E-01	9.828253E-01	9.796641E-01	9.769163E-01	9.744982E-01	0.972298829	0.970274968	0.968401913];
Phi(3,:)=[1.102908174	1.056333E+00	1.032170E+00	1.016886E+00	1.006752E+00	9.996870E-01	9.945541E-01	9.907976E-01	9.878425E-01	9.854903E-01	9.834967E-01	9.817860E-01	9.802767E-01	9.788802E-01	9.775939E-01];
Phi(4,:)=[1.104946E+00	1.054815E+00	1.030861E+00	1.016544E+00	1.007504E+00	1.001553E+00	9.975767E-01	9.948158E-01	9.928193E-01	9.913146E-01	9.901179E-01	9.891128E-01	9.882412E-01	9.874220E-01	9.866561E-01];
Phi(5,:)=[1.065936E+00	1.045374E+00	1.025142E+00	1.013327E+00	1.005468E+00	1.000860E+00	9.979917E-01	9.961581E-01	9.949503E-01	9.942173E-01	9.935876E-01	9.930632E-01	9.925565E-01	9.920685E-01	0.991620443];
Phi(6,:)=[1.10325287	1.05111062	1.027710778	1.014724731	1.007125636	1.002649781	0.999966359	0.998356329	0.997388687	0.996781804	0.996369283	0.995997682	0.995694245	0.99534514	0.995004933];
Phi(7,:)=[1.102144E+00	1.049449E+00	1.026363E+00	1.013818E+00	1.006732E+00	1.002711E+00	1.000443E+00	9.992161E-01	9.985393E-01	9.982133E-01	9.979520E-01	9.977707E-01	9.975886E-01	9.973628E-01	9.970968E-01];
Phi(8,:)=[1.093472845	1.04021298	1.018680707	1.008806997	1.004141732	1.00222237	1.001658105	1.001758844	1.002089602	1.002441487	1.002719658	1.002867087	1.002897356	1.00280074	1.002633801];
Phi(9,:)=[1.084619956	1.030996616	1.01156887	1.003998715	1.001438071	1.00106761	1.001521795	1.002157306	1.002730537	1.003137136	1.003302875	1.003297407	1.003159499	1.002896386	1.002576827];
 
Phi(:,1:2)=NaN;
Phi=Phi(:,3:end);
vp=vp(:,3:end);
d=d(:,3:end);

Recovery=((4*vp*1.46)./(0.1*d)).^2;

PhiM = a+b*(c./0.1.^n.*vp.^1*(0.75*10^(-9))^(n-1)*1.46^(p-n+1)./d.^(p)).^(f);
Error = Error + sum(sum(((Phi-PhiM)).^2.*Recovery,1),2); %error weight: recovery

%%D,vp
D=zeros(5,15); 
vp=zeros(5,15); 
Phi=zeros(5,15); 
D(1,:)=0.55*10^(-9);
D(2,:)=0.75*10^(-9);
D(3,:)=0.85*10^(-9);
D(4,:)=0.95*10^(-9);
D(5,:)=1.15*10^(-9);
vp(1,:)=[3.11479E-07	7.27607E-07	1.16315E-06	1.60347E-06	2.04407E-06	2.48311E-06	2.91944E-06	3.35257E-06	3.78193E-06	4.20663E-06	4.62669E-06	5.04121E-06	5.45022E-06	5.85283E-06	6.24905E-06];
vp(2,:)=[3.12089E-07	7.30746E-07	1.17019E-06	1.61536E-06	2.06181E-06	2.50761E-06	2.95171E-06	3.39385E-06	3.83309E-06	4.26884E-06	4.70139E-06	5.12988E-06	5.554E-06	5.97318E-06	6.38713E-06];
vp(3,:)=[3.12234E-07	7.31735E-07	1.17243E-06	1.61922E-06	2.06756E-06	2.51558E-06	2.96217E-06	3.40722E-06	3.84966E-06	4.28919E-06	4.72581E-06	5.15895E-06	5.58801E-06	6.01271E-06	6.43277E-06];
vp(4,:)=[3.12467E-07	7.32781E-07	1.17464E-06	1.62289E-06	2.07291E-06	2.52285E-06	2.97176E-06	3.41914E-06	3.86419E-06	4.30721E-06	4.74703E-06	5.18395E-06	5.61708E-06	6.04643E-06	6.47143E-06];
vp(5,:)=[3.12729E-07	7.34177E-07	1.17778E-06	1.62812E-06	2.08061E-06	2.5334E-06	2.98543E-06	3.43629E-06	3.88571E-06	4.33308E-06	4.77785E-06	5.22028E-06	5.65981E-06	6.09585E-06	6.5284E-06];
Phi(1,:)=[0.471467635	1.046284743	1.023237255	1.008408251	0.998464114	0.99140006	0.986154067	0.982096925	0.978807543	0.975993119	0.973543503	0.971319105	0.969302824	0.967420334	0.96568654];
Phi(2,:)=[1.096709306	1.055922382	1.031798391	1.015520386	1.003726413	0.994886402	0.987921011	0.982415193	0.977756305	0.973738647	0.970249793	0.967104588	0.964206532	0.961512262	0.958992465];
Phi(3,:)=[1.100556722	1.059811764	1.03540187	1.018549196	1.006138233	0.996505893	0.988901896	0.982753065	0.977512845	0.972960714	0.968987284	0.965380342	0.962061111	0.958975439	0.956078273];
Phi(4,:)=[1.102821961	1.063429917	1.039008422	1.021619217	1.008518857	0.998413797	0.990191752	0.983426705	0.977646829	0.972670563	0.968194964	0.964154165	0.960433192	0.956956883	0.953682624];
Phi(5,:)=[1.108264085	1.069184776	1.044750361	1.026743064	1.012843017	1.001734605	0.992548866	0.984761554	0.978073469	0.972216281	0.966857765	0.96203801	0.957571132	0.953365315	0.949394332];
 
Phi(:,1:2)=NaN;
Phi=Phi(:,3:end);
vp=vp(:,3:end);
D=D(:,3:end);

Recovery=((4*vp*1.46)./(0.1*0.00075)).^2;

PhiM = a+b*(c./0.1.^n.*vp.^1.*(D).^(n-1)*1.46^(p-n+1)./0.00075.^(p)).^(f);
Error = Error + sum(sum(((Phi-PhiM)).^2.*Recovery,1),2); %error weight: recovery

%%L,vp
L=zeros(7,15); 
vp=zeros(7,15); 
Phi=zeros(7,15); 
L(1,:)=0.26;
L(2,:)=0.41;
L(3,:)=0.56;
L(4,:)=0.71;
L(5,:)=1.46;
L(6,:)=1;
L(7,:)=1.25;
vp(1,:)=[3.3385E-07	7.62197E-07	1.21317E-06	1.67186E-06	2.13366E-06	2.59692E-06	3.06084E-06	3.52475E-06	3.98867E-06	4.45242E-06	4.91585E-06	5.37862E-06	5.84107E-06	6.30286E-06	6.764E-06];
vp(2,:)=[3.31001E-07	7.57733E-07	1.20668E-06	1.66308E-06	2.12227E-06	2.58271E-06	3.04346E-06	3.5041E-06	3.96443E-06	4.42415E-06	4.88324E-06	5.3415E-06	5.79893E-06	6.25533E-06	6.7108E-06];
vp(3,:)=[3.28215E-07	7.53553E-07	1.20086E-06	1.65521E-06	2.11229E-06	2.57028E-06	3.02834E-06	3.48595E-06	3.94303E-06	4.39912E-06	4.8543E-06	5.30827E-06	5.76103E-06	6.21242E-06	6.66245E-06];
vp(4,:)=[3.25465134E-07	7.49537589E-07	1.19529098E-06	1.64797845E-06	2.10299720E-06	2.55873327E-06	3.01429001E-06	3.46918921E-06	3.92307221E-06	4.37581946E-06	4.82719184E-06	5.27706981E-06	5.72527404E-06	6.17192408E-06	6.61606351E-06];
vp(5,:)=[2.7631624E-07	6.3973020E-07	1.0211788E-06	1.4083250E-06	1.7972735E-06	2.1863965E-06	2.5749090E-06	2.9621715E-06	3.3482130E-06	3.7322198E-06	4.1141916E-06	4.4938379E-06	4.8708680E-06	5.2452818E-06	5.6164980E-06];
vp(6,:)=[3.2039E-07	7.42384E-07	1.18564E-06	1.63535E-06	2.08701E-06	2.53884E-06	2.98999E-06	3.43991E-06	3.88826E-06	4.33453E-06	4.77889E-06	5.22113E-06	5.66082E-06	6.09754E-06	6.53129E-06];
vp(7,:)=[3.15899E-07	7.36068E-07	1.17725E-06	1.62448E-06	2.07328E-06	2.52186E-06	2.96923E-06	3.415E-06	3.85843E-06	4.29914E-06	4.73713E-06	5.17207E-06	5.60361E-06	6.03108E-06	6.45482E-06];
Phi(1,:)=[1.03982183	1.048396281	1.041260834	1.03319755	1.026427841	1.020831657	1.01630133	1.012531596	1.009588756	1.007290087	1.005428078	1.003943057	1.002802438	1.001920384	1.001227768];
Phi(2,:)=[1.0538068E+00	1.0509658E+00	1.0398594E+00	1.0303427E+00	1.0225673E+00	1.0165070E+00	1.0117273E+00	1.0079518E+00	1.0050288E+00	1.0027209E+00	1.0009138E+00	9.9948929E-01	9.9836728E-01	9.9747891E-01	9.9676766E-01];
Phi(3,:)=[1.061209215	1.051816379	1.03837664	1.027607741	1.019300241	1.012842617	1.007850354	1.003967927	1.000931049	0.998622978	0.99665482	0.995178648	0.993970038	0.992945028	0.992109415];
Phi(4,:)=[1.0665528E+00	1.0520430E+00	1.0368965E+00	1.0251319E+00	1.0163688E+00	1.0094893E+00	1.0041779E+00	1.0002704E+00	9.9703400E-01	9.9450107E-01	9.9246119E-01	9.9084458E-01	9.8939519E-01	9.8823576E-01	0.987123314];
Phi(5,:)=[1.085017446	1.057476273	1.036175875	1.020511024	1.009004025	0.999984548	0.992958554	0.987185415	0.982430128	0.978351816	0.974777484	0.971609698	0.968719691	0.966068551	0.963594139];
Phi(6,:)=[1.088608324	1.05666306	1.036860859	1.022661213	1.012360058	1.004628348	9.99E-01	9.94E-01	9.90E-01	9.87E-01	9.85E-01	9.83E-01	9.81E-01	0.979033888	0.977423176];
Phi(7,:)=[1.093789549	1.056091319	1.034096473	1.01888539	1.007762026	0.99937799	0.992975365	0.987902343	0.98373267	0.980205082	0.977178566	0.974505784	0.972084962	0.969836604	0.96778241];

Phi(:,1:2)=NaN;
Phi=Phi(:,3:end);
vp=vp(:,3:end);
L=L(:,3:end);

Recovery=((4*vp.*L)./(0.1*0.00075)).^2;

PhiM = a+b*(c./0.1.^n.*vp.^1*(0.75*10^(-9))^(n-1).*L.^(p-n+1)./0.00075.^(p)).^(f);
Error = Error + sum(sum(((Phi-PhiM)).^2.*Recovery,1),2); %error weight: recovery

Result=Error;
end