c7chua.m

来自「高等引用数学的matlab求解的配套程序」· M 代码 · 共 5 行

function dx=c7chua(t,x,flag,a,b,alpha,beta,gamma)
dx=[alpha*(x(2)-x(1)-nlfun(x(1),a,b));
    x(1)-x(2)+x(3); -beta*x(2)-gamma*x(3)];
function f0=nlfun(x,a,b)
f0=b*x+(a-b)*(abs(x+1)-abs(x-1))/2;