clear all;
clc;
format short;

cc(1,:)=[20 24 25 22 23 28 23 24 25 27 20 26 20 28 24 26 21 23 19 22 23 23 21 22 24 26 24 27 25 21 24 21 23 22 22 24 26 27 28 24 21 20 22 25 24 26 24 25 24 21 20 23 22 22 22 22 23 22 24 25 27 21 26 28 24 25 24 24 27 27 29 23 24 25 21 24 21 23 26 25 19 28 24 26 27 27 29 22 24 18 23 24 26 25 26 27];
cc(2,:)=[4 3 2 3 4 3 2 3 4 3 2 3 4 3 2 1 2 2 3 4 2 3 5 3 2 1 2 2 4 3 1 2 4 3 2 3 4 3 3 2 1 3 3 4 2 1 2 3 3 2 3 4 4 3 1 1 3 2 1 2 3 3 4 3 2 2 1 2 3 4 3 2 1 2 3 3 2 4 2 3 4 2 3 2 1 2 3 4 3 3 4 2 2 3 3 4];
cc(3,:)=[2 1 2 3 3 4 3 2 2 1 2 3 4 3 2 1 2 3 3 2 4 2 3 4 2 3 2 1 2 3 4 3 3 4 2 2 3 3 4 4 3 2 3 4 3 2 3 4 3 2 3 4 3 2 1 2 2 3 4 2 3 4 3 2 1 2 2 4 3 1 2 4 3 2 3 4 3 3 2 1 3 3 4 2 1 2 3 3 2 3 4 4 3 1 1 3];

v(1,1:9)=[466100 466100 466100 466100 466100 1372500 1455700 1491400 1504300];
h(1,1:9)=[1166 1166 1166 1166 1166 1236 1240 1241.87 1242.51];

v(2,1:9)=[66290 70140 74140 78310 82820 85120 87420 89720 92020];
v(3,1:9)=[33800 37000 40800 46100 504000 52500 74000 88400 106200];
v(3,200:208)=[33800 37000 40800 46100 504000 52500 74000 88400 106200];

h(2,1:9)=[982 984 986 988 990 991 992 993 994];
h(3,1:9)=[880 882 884 887 889 890 899 904.2 910];

n=3;

a=polyfit(v(3,1:9),h(3,1:9),n);

for tt=1:0
    a(tt)
end

p(1:96)=[2213.182 2248.814 1951.312 1698.104 1116.228 1029.144 1005.312 1019.088 1017.408 1020.768 1021.104 1032.192 1023.12 1000.272 1014.048 1011.696 999.6 1013.376 1015.392 1021.776 1037.232 1046.304 1039.584 1038.912 1039.248 1037.568 1081.14 1418.366 1680.136 1730.278 1872.286 2192.83 2566.202 3185.16 3414.376 3917.122 3911.432 4442.352 4373.522 4485.124 4524.184 4585.402 4421.82 4446.464 4463.276 4513.4 4433.202 4631.592 4441.916 4248.732 4246.436 4119.032 4166.74 4298.232 4275.5 4375.492 4437.34 4488.23 4484.714 4482.696 4427.826 4419.074 4443.042 4407.496 3899.074 3928.91 4121.99 4291.738 4356.966 4390.434 4029.394 4058.2 3921.746 3764.098 3503.244 3605.764 4065.02 4571.61 4655.234 4712.454 4721.958 4814.112 4817.884 4798.02 4697.89 4647.854 3962.894 3965.158 3533.862 3117.786 3057.966 3088.324 3285.392 2663.634 1454.554 1249.32];

nnn=1;
nnn=10.9218;

q(1,193:200)=[1 1 1 1 1 1 1 1];
q(2,195:200)=[1 1 1 1 1 1];

opt(1:3)=[0 0 0];

qmax1=200;
qmax2=200;
qmax3=200;
qmin1=1;
qmin2=1;
qmin3=1;

q(1,1)=20;q(2,1)=2;

for jj=1:9

    %q(1,jj)=(p(jj)*nnn-(36.9919+0.0007061*v(3,jj)-2.3679*10^(-12)*v(3,jj)^2+2.1528*10^(-15)*v(3,jj)^3)*q(3,jj)-(5.2975+0.002023802*v(2,jj)-1.760249*10^(-8)*v(2,jj)^2+6.4832*10^(-14)*v(2,jj)^3)*q(2,jj))/(64.6+0.000311319*v(1,jj+199)-2.0832*10^(-10)*v(1,jj+199)^2+5.4337*10^(-17)*v(1,jj+199)^3);
    %*v(1,jj)^2)*q(1,jj)-(46.10616+0.626623*10^(-6)*v(2,jj)-0.156644*10^(-13)*v(2,jj)^2+0.19337127*10^(-22)*v(2,jj)^3-0.118*10^(-31)*v(2,jj)^4)*q(2,jj))/(26.651562+0.119706*10^(-6)*v(3,jj+199)-0.122998*10^(-15)*v(3,jj+199)^2+0.84559*10^(-25)*v(3,jj+199)^3);

    fitness_function=64.6*v(1,jj)+0.00015566*v(1,jj)^2-0.6944*10^(-10)*v(1,jj)^3+1.3594*10^(-17)*v(1,jj)^4+5.2975*v(2,jj)+0.001011901*v(2,jj)^2-0.58675*10^(-8)*v(2,jj)^3+1.6208*10^(-14)*v(2,jj)^4+36.9919*v(3,jj)+0.000353*v(3,jj)^2-0.7893*10^(-9)*v(3,jj)^3+0.5382*10^(-15)*v(3,jj)^4

    p(jj)*nnn;

    d=64.6+0.00311319*v(1,jj)-2.0832*10^(-8)*v(1,jj)^2+5.4337*10^(-14)*v(1,jj)^3;
    b=5.2975+0.02023802*v(2,jj)-1.760249*10^(-6)*v(2,jj)^2+6.4832*10^(-11)*v(2,jj)^3;
    c=36.9919+0.0071*v(3,jj+199)-2.3679*10^(-7)*v(3,jj+199)^2+2.1528*10^(-12)*v(3,jj+199)^3;
    
end