| 编程中国 | 业界新闻 | 技术文章 | 视频教程 | 下载频道 | 程序源码 | 个人空间 | 编程论坛
全能ASP/PHP/ASP.NET主机,支持月付专业 MSSQL 数据库空间,支持月付专业 MySQL 数据库空间,支持月付买域名,送MP3、MP4
高端软件开发 = 年薪十万不是梦赛孚耐:软件保护加密专家身份认证令牌USB KEY买空间,免费送域名(厦门中资源)
共有 354 人关注过本帖
标题:求助 M 文件编写错误
收藏  订阅  推荐  打印 
hebutyll
Rank: 1
等级:新手上路
帖子:8
积分:196
注册:2008-7-29
求助 M 文件编写错误

这是我编写的M文件,但是不知道哪里错了,还请帮助!
%jacobi.m
function s=jacobi (A, b, X0, eps)
if nargin==3
    eps=1.0e-6
elseif nargin<3
    error
    return
end
% 求对角阵D和
D=diag(diag(A));
D=inv(D);
L=tril(A, -1);
U=triu(A, 1);
B=-D*(L+U);
f=D*b;
%进行迭代求解
while norm (s-X0) >=eps
    X0=s;
    s=B*X0+f;
end
return
当在 matlab中输入
a=[10 -1 2;2 5 -1;-1 2 5];
b=[5 10 15];
x0=[0 0 0];
eps=0.0001;
s=jacobi(a, b, x0, eps)
就会出现错误提示。
搜索更多相关主题的帖子: eps  文件  编写  jacobi  nargin  
2008-7-29 11:26
关于我们 | 广告合作 | 编程中国 | 清除Cookies | Archiver | WAP | TOP

编程中国 版权所有,并保留所有权利。鲁ICP备08000592号
Powered by Discuz, Processed in 0.048331 second(s), 9 queries.
Copyright©2004-2008, BCCN.NET, All Rights Reserved