请有时间的高手帮忙运行一下代码，thanx！

kevinhcj

%subrutine<BR>function C_b1 = element_Cb1(z0,z1,t,bb)<BR>syms Xi <BR>psi = bb;
<P>N1 = 1-10*Xi^3+15*Xi^4-6*Xi^5;<BR>N2 = 10*Xi^3-15*Xi^4+6*Xi^5;<BR>N3 = Xi-6*Xi^3+8*Xi^4-3*Xi^5;<BR>N4 = -4*Xi^3+7*Xi^4-3*Xi^5;<BR>N5 = 1/2*(Xi^2-3*Xi^3+3*Xi^4-Xi^5);<BR>N6 = 1/2*(1*Xi^3-2*Xi^4+Xi^5);</P>
<P>N  = [N1,N3,N5,N2,N4,N6];</P>
<P>NN1 = diff(N,Xi,1);<BR>cc = (psi'*(NN1'*NN1)*psi)*(N'*N);<BR>C_b1 = int(cc,Xi,0,1);<BR>C_b1 = C_b1;<BR>C_b1 = double(C_b1);</P>
<P>return</P>
<P>% main<BR>clear all; clc;<BR>format long;</P>
<P>n = 10<BR>m = 100;<BR>zz = [0,1:n]';<BR>delta_t = 0.1;</P>
<P>C_b1_u = zeros(3*n+3);<BR>psi_q_try = ones(3*(n+1),1);</P>
<P>for j = 1:m<BR>    <BR>    t = delta_t*j;<BR>       <BR>        for k = 1:n<BR>            bb = psi_q_try(3*k-2:3*(k+1));<BR>            C_b1 = element_Cb1(zz(k),zz(k+1),t,bb);<BR>            C_b1_u(3*k-2:3*(k+1),3*k-2:3*(k+1)) = C_b1_u(3*k-2:3*(k+1),3*k-2:3*(k+1)) + C_b1;<BR>        end <BR>end<BR><BR>在同一设置下运行多次，所能执行的循环数不一致！</P>

MVH

运行到j和k都等于9的时候停止

LILI-C

不好意思楼主 我还想问一下 你说的那个矩阵积分的方法 算4阶矩阵的时候 10分钟都积不出来 是这样的么

kevinhcj

int(FUN,X,a,b)<BR>其中FUN可以是矩阵，使用就跟单个函数一样。至于所需的时间，看矩阵的内容而定了，要不你重新发个帖子，把代码帖出来，让我看看。

happy

<DIV class=quote><B>以下是引用<I>kevinhcj</I>在2006-3-15 18:00:47的发言：</B><BR>
<P><BR>在同一设置下运行多次，所能执行的循环数不一致！</P></DIV>
<P><BR>试了两次，结果都和2楼一样，不存在循环数不一致问题</P>

happy

另外有一个疑问，你的子程序中根本和z0,z1,t这几个变量无关，为什么还要传递呢？