function [f,x]=myga(num,bounds,Myfun,N,CP,P)
%[f,x]=ga(num,bounds,fun,N,CP,P)
%該遺傳算法適用于：
%           目標(biāo)函數(shù)為求最大值，且解非負(fù)整數(shù)解
%bounds     邊界約束
%Myfun        為目標(biāo)函數(shù)
%num        初始種群數(shù)
%N          最大迭代次數(shù)
%CP         交叉概率
%P          突變概率
%f          目標(biāo)最優(yōu)解
%x          最優(yōu)解向量

m=nargin;
if m<6
disp('-_-  錯(cuò)誤!')
disp('>> 輸入變量太少')
disp('>>  按回車鍵查看幫助')
    pause
    help ga
    f='-_- ';
    x='沒(méi)有規(guī)矩不成方圓';
    break;
end
pop=INTinti(num,bounds);
fmax=pop(:,end);
endpop=pop;
n=size(endpop,2);
k=0;x=[];f=zeros(1,num);
while(k<N)
    pop=mutation(endpop);
    [cpop ,len,v]=cross(pop,bounds,CP);
    [pops]=changes(cpop,bounds,len,P);break;
    for i=1:num
        sol=pops(i,:);
        [f(i)]=Myfun(sol);
        if fmax(i)<f(i)
            fmax(i)=f(i);
            endpop(i,1:end-1)=pops(i,:);
        end
end
endpop(:,end)=fmax(:);
k=k+1;
end
[f,ii]=max(fmax);
x=endpop(ii,1:end-1);