?? lesliegower.m
字號:
function out = LeslieGowerout{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;out{10} = @userf1;% --------------------------------------------------------------------------function dydt = fun_eval(t,kmrgd,b1,b2,ma,mj,my,mz,cjj,cja,cyy,cyz,cjy,cjz,cyj,cya)dydt=[(b1*kmrgd(2))/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4)),;(1-mj)*kmrgd(1)+(1-ma)*kmrgd(2),;(b2*kmrgd(4))/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4)),;(1-my)*kmrgd(3)+(1-mz)*kmrgd(4),;];% --------------------------------------------------------------------------function [tspan,y0,options] = inithandles = feval(LeslieGower);y0=[0,0,0,0];options = odeset('Jacobian',handles(3),'JacobianP',handles(4),'Hessians',handles(5),'HessiansP',handles(6));tspan = [0 10];% --------------------------------------------------------------------------function jac = jacobian(t,kmrgd,b1,b2,ma,mj,my,mz,cjj,cja,cyy,cyz,cjy,cjz,cyj,cya)jac=[[-b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^2*cjj,b1/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))-b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^2*cja,-b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^2*cjy,-b1*kmrgd(2)/(1+cjj*kmrgd(1)+cja*kmrgd(2)+cjy*kmrgd(3)+cjz*kmrgd(4))^2*cjz];[1-mj,1-ma,0,0];[-b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^2*cyj,-b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^2*cya,-b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^2*cyy,b2/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))-b2*kmrgd(4)/(1+cyj*kmrgd(1)+cya*kmrgd(2)+cyy*kmrgd(3)+cyz*kmrgd(4))^2*cyz];[0,0,1-my,1-mz]];?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -