function [sol]=fitness(sol,data)
% num min max a b c e f
Infoo=[1 36 114 94.705 6.73 0.00690 100 0.0840
2 36 114 94.705 6.73 0.00690 100 0.0840
3 60 120 309.540 7.07 0.02028 100 0.0840
4 80 190 369.030 8.18 0.00942 150 0.0630
5 47 97 148.890 5.35 0.01140 120 0.0770
6 68 140 222.330 8.05 0.01142 100 0.0840
7 110 300 287.710 8.03 0.00357 200 0.0420
8 135 300 391.980 6.99 0.00492 200 0.0420
9 135 300 455.760 6.6 0.00573 200 0.0420
10 130 300 722.820 12.9 0.00605 200 0.0420
11 94 375 635.200 12.9 0.00515 200 0.0420
12 94 375 654.690 12.8 0.00569 200 0.0420
13 125 500 913.400 12.5 0.00421 300 0.0350
14 125 500 1760.400 8.84 0.00752 300 0.0350
15 125 500 1728.300 9.15 0.00708 300 0.0350
16 125 500 1728.300 9.15 0.00708 300 0.0350
17 220 500 647.850 7.97 0.00313 300 0.0350
18 220 500 649.690 7.95 0.00313 300 0.0350
19 242 550 647.830 7.97 0.00313 300 0.0350
20 242 550 647.818 7.97 0.00313 300 0.0350
21 254 550 785.960 6.63 0.00298 300 0.0350
22 254 550 785.960 6.63 0.00298 300 0.0350
23 254 550 794.530 6.66 0.00284 300 0.0350
24 254 550 794.530 6.66 0.00284 300 0.0350
25 254 550 801.320 7.1 0.00277 300 0.0350
26 254 550 801.320 7.1 0.00277 300 0.0350
27 10 150 1055.100 3.33 0.52124 120 0.0770
28 10 150 1055.100 3.33 0.52124 120 0.0770
29 10 150 1055.100 3.33 0.52124 120 0.0770
30 47 97 148.890 5.35 0.01140 120 0.0770
31 60 190 222.920 6.43 0.00160 150 0.0630
32 60 190 222.920 6.43 0.00160 150 0.0630
33 60 190 222.920 6.43 0.00160 150 0.0630
34 90 200 107.870 8.95 0.00010 200 0.0420
35 90 200 116.580 8.62 0.00010 200 0.0420
36 90 200 116.580 8.62 0.00010 200 0.0420
37 25 110 304.450 5.88 0.01610 80 0.0980
38 25 110 304.450 5.88 0.01610 80 0.0980
39 25 110 304.450 5.88 0.01610 80 0.0980
40 242 550 647.830 7.97 0.00313 300 0.0350];
PD=9875
pmin=Infoo(:,2)';
pmax=Infoo(:,3)';
a=Infoo(:,4)';
b=Infoo(:,5)';
c=Infoo(:,6)';
e=Infoo(:,7)';
f=Infoo(:,8)';
Npar=data.Npar;
lb=data.lb;
ub=data.ub;
emp.x=[];
emp.v=[];
emp.fit=[];
par=repmat(emp,Npar,1);
for i=1:Npar
par(i).x=unifrnd(lb.x,ub.x);
par(i).v=0;
par(i)=fitness(par(i),data);
end
for j=1:40
cost=a(j)*x(i).^2+b(j)*x(i)+c(j)+abs(e(j)*sin(f(j)*(pmin(j)-x(i))));
sol.fit=sum(cost);
end