私はこの方程式を持っています:f(a,b,x)=t0-a+(a^2*(1+((x-x0)^2/b^2)))^0.5
一次導関数を取得したい場合:df(a,b,x)/d(a,b)
for x0,t0
=constantおよびa,b,x
=array。
このスクリプトを使用する場合、これは本当ですか
> syms f; syms t0; syms x; syms x0; syms
> a; syms b;
>
>
> f=t0-a+(a^2*(1+((x-x0)^2/b^2)))^0.5;
> f_1=diff(f,a)+diff(f,b)
この方法を使用すると、正しいプロットが得られないため、質問しています...何か提案はありますか?