| 查看: 900 | 回复: 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)))']; 这是代码，总是报错，帮忙看一下吧 |
回复此楼» 猜你喜欢
求调剂
已经有3人回复
-
265求调剂
已经有4人回复
-
085700资源与环境308求调剂
已经有6人回复
-
一志愿吉林大学材料学硕321求调剂
已经有12人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有3人回复
-
329求调剂
已经有5人回复
-
申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复
已经有4人回复
-
材料学硕318求调剂
已经有5人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有6人回复
-
081700化工学硕调剂
已经有3人回复
仙木映月
金虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 2257.6
- 散金: 51
- 红花: 17
- 帖子: 608
- 在线: 1094.4小时
- 虫号: 1929532
- 注册: 2012-08-08
- 性别: GG
- 专业: 机器人学及机器人技术
2楼2014-05-30 21:19:01
3楼2014-05-30 21:29:34
charmantang
禁虫 (正式写手)
★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2014-06-04 00:47:27
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2014-06-04 00:47:27
|
本帖内容被屏蔽 |
4楼2014-06-01 18:18:15
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
gzy_1035楼2014-06-02 14:30:13
6楼2014-06-03 09:53:35
方勤1
金虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 499.1
- 帖子: 45
- 在线: 24.7小时
- 虫号: 2326342
- 注册: 2013-03-07
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
7楼2014-06-03 22:13:42
8楼2014-06-04 11:28:32
方勤1
金虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 499.1
- 帖子: 45
- 在线: 24.7小时
- 虫号: 2326342
- 注册: 2013-03-07
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
9楼2014-06-06 08:30:14
