function out = Ghmapout{1} = @init;out{2} = @fun_eval;out{3} = @jacobian;out{4} = @jacobianp;out{5} = @hessians;out{6} = @hessiansp;out{7} = @der3;out{8} = @der4;out{9} = @der5;% --------------------------------------------------------------------------function dydt = fun_eval(t,kmrgd,alpha,beta,r,s)dydt=[kmrgd(2),;alpha-beta*kmrgd(1)-kmrgd(2)^2+r*kmrgd(1)*kmrgd(2)+s*kmrgd(2)^3,;];% --------------------------------------------------------------------------function [tspan,y0,options] = inithandles = feval(Ghmap);y0=[0,0];options = odeset('Jacobian',handles(3),'JacobianP',handles(4),'Hessians',handles(5),'HessiansP',handles(6));tspan = [0 10];% --------------------------------------------------------------------------function jac = jacobian(t,kmrgd,alpha,beta,r,s)jac=[[0,1];[-beta+r*kmrgd(2),-2*kmrgd(2)+r*kmrgd(1)+3*s*kmrgd(2)^2]];% --------------------------------------------------------------------------function jacp = jacobianp(t,kmrgd,alpha,beta,r,s)jacp=[[0,0,0,0];[1,-kmrgd(1),kmrgd(1)*kmrgd(2),kmrgd(2)^3]];% --------------------------------------------------------------------------function hess = hessians(t,kmrgd,alpha,beta,r,s)hess1=[[0,0];[0,r]];hess2=[[0,0];[r,-2+6*s*kmrgd(2)]];hess(:,:,1) =hess1;hess(:,:,2) =hess2;% --------------------------------------------------------------------------function hessp = hessiansp(t,kmrgd,alpha,beta,r,s)hessp1=[[0,0];[0,0]];hessp2=[[0,0];[-1,0]];hessp3=[[0,0];[kmrgd(2),kmrgd(1)]];hessp4=[[0,0];[0,3*kmrgd(2)^2]];hessp(:,:,1) =hessp1;hessp(:,:,2) =hessp2;hessp(:,:,3) =hessp3;hessp(:,:,4) =hessp4;%---------------------------------------------------------------------------function tens3  = der3(t,kmrgd,alpha,beta,r,s)tens31=[[0,0];[0,0]];tens32=[[0,0];[0,0]];tens33=[[0,0];[0,0]];tens34=[[0,0];[0,6*s]];tens3(:,:,1,1) =tens31;tens3(:,:,1,2) =tens32;tens3(:,:,2,1) =tens33;tens3(:,:,2,2) =tens34;%---------------------------------------------------------------------------function tens4  = der4(t,kmrgd,alpha,beta,r,s)tens41=[[0,0];[0,0]];tens42=[[0,0];[0,0]];tens43=[[0,0];[0,0]];tens44=[[0,0];[0,0]];tens45=[[0,0];[0,0]];tens46=[[0,0];[0,0]];tens47=[[0,0];[0,0]];tens48=[[0,0];[0,0]];tens4(:,:,1,1,1) =tens41;tens4(:,:,1,1,2) =tens42;tens4(:,:,1,2,1) =tens43;tens4(:,:,1,2,2) =tens44;tens4(:,:,2,1,1) =tens45;tens4(:,:,2,1,2) =tens46;tens4(:,:,2,2,1) =tens47;tens4(:,:,2,2,2) =tens48;%---------------------------------------------------------------------------function tens5  = der5(t,kmrgd,alpha,beta,r,s)tens51=[[0,0];[0,0]];tens52=[[0,0];[0,0]];tens53=[[0,0];[0,0]];tens54=[[0,0];[0,0]];tens55=[[0,0];[0,0]];tens56=[[0,0];[0,0]];tens57=[[0,0];[0,0]];tens58=[[0,0];[0,0]];tens59=[[0,0];[0,0]];tens510=[[0,0];[0,0]];tens511=[[0,0];[0,0]];tens512=[[0,0];[0,0]];tens513=[[0,0];[0,0]];tens514=[[0,0];[0,0]];tens515=[[0,0];[0,0]];tens516=[[0,0];[0,0]];tens5(:,:,1,1,1,1) =tens51;tens5(:,:,1,1,1,2) =tens52;tens5(:,:,1,1,2,1) =tens53;tens5(:,:,1,1,2,2) =tens54;tens5(:,:,1,2,1,1) =tens55;tens5(:,:,1,2,1,2) =tens56;tens5(:,:,1,2,2,1) =tens57;tens5(:,:,1,2,2,2) =tens58;tens5(:,:,2,1,1,1) =tens59;tens5(:,:,2,1,1,2) =tens510;tens5(:,:,2,1,2,1) =tens511;tens5(:,:,2,1,2,2) =tens512;tens5(:,:,2,2,1,1) =tens513;tens5(:,:,2,2,1,2) =tens514;tens5(:,:,2,2,2,1) =tens515;tens5(:,:,2,2,2,2) =tens516;