func.m

这次上传的绝对是好东西

% func 
% ---- 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% function [fvec,J]=func(n,m,x,nprob,option)           
% Selects the appropriate test function based on nprob.
%                                                      
% Created on 10/22/94 by Madhu Lamba                   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


function [fvec,J] = func(n,m,x,nprob,opt)

if (nprob==1)
	funcname = 'rose';
    [fvec,J]=rose(n,m,x,opt);
elseif (nprob==2)
        funcname = 'froth';[fvec,J]=froth(n,m,x,opt);
elseif (nprob==3)
        funcname = 'badscp';[fvec,J]=badscp(n,m,x,opt);
elseif (nprob==4)
        funcname = 'badscb';[fvec,J]=badscb(n,m,x,opt);
elseif (nprob==5)
        funcname = 'beale';[fvec,J]=beale(n,m,x,opt);
elseif (nprob==6)
        funcname = 'jensam';[fvec,J]=jensam(n,m,x,opt);
elseif (nprob==7)
        funcname = 'helix';[fvec,J]=helix(n,m,x,opt);
elseif (nprob==8)
        funcname = 'bard';[fvec,J]=bard(n,m,x,opt);
elseif (nprob==9)
        funcname = 'gauss';[fvec,J]=gauss(n,m,x,opt);
elseif (nprob==10)
        funcname = 'meyer';[fvec,J]=meyer(n,m,x,opt);
elseif (nprob==11)        
        funcname = 'gulf';[fvec,J]=gulf(n,m,x,opt);
elseif (nprob==12)   
        funcname = 'box';[fvec,J]=box(n,m,x,opt);
elseif (nprob==13)
        funcname = 'sing';[fvec,J]=sing(n,m,x,opt);
elseif (nprob==14)
        funcname = 'wood';[fvec,J]=wood(n,m,x,opt);
elseif (nprob==15)
        funcname = 'kowosb';[fvec,J]=kowosb(n,m,x,opt);
elseif (nprob==16)       
        funcname = 'bd';[fvec,J]=bd(n,m,x,opt);
elseif (nprob==17)
        funcname = 'osb1';[fvec,J]=osb1(n,m,x,opt);
elseif (nprob==18)
        funcname = 'biggs';[fvec,J]=biggs(n,m,x,opt);
elseif (nprob==19)
        funcname = 'osb2';[fvec,J]=osb2(n,m,x,opt);
elseif (nprob==20)
        funcname = 'watson';[fvec,J]=watson(n,m,x,opt);
elseif (nprob==21)
        funcname = 'rosex';[fvec,J]=rosex(n,m,x,opt);
elseif (nprob==22)
        funcname = 'singx';[fvec,J]=singx(n,m,x,opt);
elseif (nprob==23)
        funcname = 'pen1';[fvec,J]=pen1(n,m,x,opt);
elseif (nprob==24)
        funcname = 'pen2';[fvec,J]=pen2(n,m,x,opt);
elseif (nprob==25)
        funcname = 'vardim';[fvec,J]=vardim(n,m,x,opt);
elseif (nprob==26)
        funcname = 'trig';[fvec,J]=trig(n,m,x,opt);
        %elseif (nprob==27)
       % funcname = 'almost';[fvec,J]=almost(n,m,x,opt);
elseif (nprob==28)
        funcname = 'bv';[fvec,J]=bv(n,m,x,opt);
elseif (nprob==29)
        funcname = 'ie';[fvec,J]=ie(n,m,x,opt);
elseif (nprob==30)
        funcname = 'trid';[fvec,J]=trid(n,m,x,opt);
elseif (nprob==31)
        funcname = 'band';[fvec,J]=band(n,m,x,opt);
elseif (nprob==32)        
        funcname = 'lin';[fvec,J]=lin(n,m,x,opt);
elseif (nprob==33)        
        funcname = 'lin1';[fvec,J]=lin1(n,m,x,opt);
elseif (nprob==34)
        funcname = 'lin0';[fvec,J]=lin0(n,m,x,opt);
end;
%[fvec,J] = feval(funcname,n,m,x,opt);
       
%