That would just be 6 lines of code in the function defining your DE:
function dxy = deFcn(t,xy)
dxy(1) = C1 * (xy0 - xy(1)) - C2*(xy(1)-xy(4)/C3)
:
:
dxy(6) = C4 * (xy7 - xy(6)) - C2*(xy(3)-xy(6)/C3)
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!