| 24小時(shí)熱門版塊排行榜 |
| 9 | 1/1 | 返回列表 |
| 查看: 911 | 回復(fù): 8 | ||
[求助]
求解方程組 已有2人參與
|
|
function y=broyden(x0) a=eye(length(x0)); x1=x0-myfun2(x0)/a; n=1; while(norm(x1-x0)>=1.0e-6)&(n<=100000000) x0=x1; x1=x0-myfun2(x0)/a; p=x1-x0; q=myfun2(x1)-myfun(x0); a=a+(q-p*a)'*p/norm(p); n=n+1; end y=x1; n function y= myfun2(x) a=2.6; y=['x(1)*0.4-0.21+x(2)*0.6','x(2)-x(1)^(1/(1-a))*(2/3)*((1-x(1))^(a/(a-1))*(0.21/0.79)^(a/(a-1))-x(1)^(1/(1-a)))']; 這是代碼,總是報(bào)錯(cuò),幫忙看一下吧 |
金蟲(chóng) (正式寫手)

禁蟲(chóng) (正式寫手)
|
本帖內(nèi)容被屏蔽 |
鐵桿木蟲(chóng) (著名寫手)
方丈大師
金蟲(chóng) (初入文壇)

金蟲(chóng) (初入文壇)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|