在Catalysis.jl中使用addreaction！函数向反应网络添加反应的正确语法

问题描述

我最近开始考虑将模拟转换为Julia，并了解了Catalyst.jl程序包，该程序包具有用于创建和修改反应网络的便捷宏。
我想使用addreaction！（）函数向已经定义的反应网络中添加特定的反应，但是无法使用提到的语法来这样做。 这是我的网络：

@reaction_func Hsn(thr,N,lbd,cop) = ((thr^cop)/((thr^cop) + (N^cop))) - lbd*(1-((thr^cop)/((thr^cop) + (N^cop))))
@reaction_func Hsp(thr,cop) = (((thr^cop)/((thr^cop) + (N^cop))) - lbd*(1-((thr^cop)/((thr^cop) + (N^cop)))))/lbd
coupled_switch = @reaction_network begin
    (kA,kB,kC,kD),(A,B,C,D) → ∅
    gA*(Hsn(tBA,fBA,nBA)*Hsn(tCA,fCA,nCA)),∅ → A
    gB*(Hsn(tAB,A,fAB,nAB)*Hsn(tDB,fDB,nDB)),∅ → B
    gC*(Hsn(tAC,fAC,nAC)*Hsn(tDC,D,fDC,nDC)),∅ → C
    #gD*(Hsn(tBD,fBD,nBD)*Hsn(tCD,fCD,nCD)),∅ → D ==> This is the reaction I am trying to add as a test.
end gA gB gC gD kA kB kC kD tBA nBA fBA tCA nCA fCA tAB nAB fAB tDB nDB fDB tAC nAC fAC tDC nDC fDC tBD nBD fBD tCD nCD fCD

并添加我正在做的反应：

addreaction!(coupled_switch,gD*(Hsn(tBD,∅ → D )

我收到以下错误：

UndefVarError: Hsn not defined
in top-level scope at r.jl:45

我意识到这是完全错误的，但是我无法在文档中找到正确的方法。 该文档没有提及要传递给addreaction！（）函数的反应的格式。有人知道正确的方法吗？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）