for jj=1:jjend %你的程序 for j=3:k-1 i2(j)=(vU2(j+1)-vU2(j))*12*C/(23*h)+16*i2(j-1)/23-5*i2(j-2)/23;%%%%%未知数在这里 是C end i2(k)=(vU2(k)-vU2(k-1))*C/h;%%%%%未知数在这里 是C i2=i2'; i4=-i2; u9(1)=0;u9(2)=0; for j=3:k-1 u9(j+1)=(i4(j+1)-i4(j))*12*L2/(23*h)+16*u9(j-1)/23-5*u9(j-2)/23; end u9(k)=(i4(k)-i4(k-1))*L2/h; u9=u9'; u4=vU2-u9-R2*i4; %比较 if norm(u4-vU1)<errr break else c=c+cinc end