matlab - Matlab での関数の最大化 (fsolve と導関数の使用)

Question

関数を最大化するベクトルp（つまり）の値を見つけたいと思います。 p(1),p(2),p(3),...A(p)

私はそれを行うためにMATLABを使用していますが、fsolveそれが役立つと思われるものを見つけました。だから私は機能を作りましたA：

function A = myfun(p)

    R  = 0.1; 
    u1 = 500;
    u2 = 400;
    u3 = 300;

    A = ( (p(1)+p(2)+p(3)) * (1/u1+1/u2+1/u3)) * ...
        (1 + R*(p(1)^2+p(2)^2+p(3)^2) * (1/u1+1/u2+1/u3) );

そして、次のような連立方程式を解く必要があります。

diff(A,p(1))==0

diff(A,p(2))==0

diff(A,p(3))==0

結果のpベクトルが私の問題の解決策になります。

fsolveこの連立方程式 ( ) をどのように解くことができp0=[1 1 1]ますか?