コメントありがとうございます。 調査した所、この手法で限界効果を直接求めることは難しいかもしれません・・・ ただ、導出された関数を効用を見たい変数で偏微分すれば良いということならば、例えばPythonのsympyライブラリを用いれば計算可能です。