residual.m

来自「MDPSAS工具箱是马里兰大学开发的」· M 代码 · 共 9 行

function A = residual(A)

unpack(A)

% dySiH4/dt
dydt = (Q*C - (Q*C + kg*C*V*y)*y - kg*C*V*y) / (C*V);

A = set(A,'resid',{dydt});