关于程序的调试

jgchen1973

我想编一段程序,求半径为从1到r的,高度从１到z变化的圆柱体的体积，输出的时候按照第一行为半径，第一列为高度，其余的地方为体积，我编的程序如下，总是不行，请高手指点，谢谢！

function [v,b]=fmianji(r,z)
fid1=fopen('mianji1.dat','a+');
%fid2=fopen('mianji2.dat','a+');
c=zeros(r+1,z+1);
for k=1:r;
    c(1,k+1)=k;
    s(k)=pi*k^2;
    l(k)=2*pi*k;
       for h=1:z;
          c(h+1,1)=h;
          v(k)=s(k)*h;
          b(k)=2*s(k)+l(k)*h;
          c(h+1,k+1)=v(k);
          fprintf(fid1,'% f\n',c);
          %fprintf(fid2,'% f % f\n',y,b(k));
    end
end
fclose(fid1);
%fclose(fid2);

[ 本帖最后由 sigma665 于 2008-5-17 16:37 编辑 ]

sigma665

有什么错误提示吗
或者运行的结果

ch_j1985

原帖由 jgchen1973 于 2008-5-17 15:15 发表

                               
登录/注册后可看大图

我想编一段程序,求半径为从1到r的,高度从１到z变化的圆柱体的体积，输出的时候按照第一行为半径，第一列为高度，其余的地方为体积，我编的程序如下，总是不行，请高手指点，谢谢！

function [v,b]=fmianji(r,z)
...

参考此帖可以解决您的问题：
http://www.chinavib.com/forum/viewthread.php?tid=54241&page=1#pid340884

happy

输出的格式游戏，呵呵

1. function [v,b]=fmianji(r,z)
   
2. fid1=fopen('mianji2.dat','a+');
   
3. for i=1:r+1
   
4.     fprintf(fid1,' %f',i-1);
   
5. end
   
6. fprintf(fid1,'\n');
   
7. for h=1:z;
   
8.     fprintf(fid1,' %f',h);
   
9.     for k=1:r
   
10.         v=pi*k^2*h;
    
11.         fprintf(fid1,' %f',v);
    
12.     end
    
13.     fprintf(fid1,'\n');
    
14. end
    
15. fclose(fid1);

复制代码

jgchen1973

谢谢教授,可以了而且完全符合我的要求!!!!