Matlab code for a diffyq. In file "system.m": function [dstatedt] = system(t,state,p1,p2) dAdt=...; dBdt=...; dstatedt=[dAdt;dBdt]; From the command line: [t,y]=ode45(@system,[tstart tend],[A0 B0],p1,p2);