| 查看: 373 | 回复: 0 | ||
[求助]
自定义子程序调用问题求解,求助大神
|
|
这个是主程序 N=129; X0=-4; XE=1.4; X=zeros(1,N); P=zeros(1,N); H=zeros(1,N); V=zeros(1,N); DX=(XE-X0)/(N-1); for i=1:N X(i)=X0+(i-1)*DX; H(i)=0.5*X(i)*Xi); if abs(X(i))>1 P(i)=0; else P(i)=sqrt(1-(X(i))^2); end end Ak=Ak(N); V=V(Ak,X0,XE,N,P); Ak的子程序昨天求贴运行过,单独计算是可以的function Ak=Ak(N) i=1:N; Ak(i)=(i+0.5).*(log(abs(i+0.5))-1)-(i-0.5).*(log(abs(i-0.5))-1); 下面是计算v的子程序 function V=V(Ak,X0,XE,N,P) for i=1:N; PAI1=0.318309886;%pi的倒数% C=log((XE-X0)/(N-1)); V(1)=0; j=1; for j=1:N; IJ=abs(i-j); if IJ>0; V(i)=V(i)+(Ak(IJ)+C)*((XE-X0)/(N-1))*P(j); end end end 每次一运行就提示 Subscript indices must either be real positive integers or logicals. 实在不知道什么原因,希望有大神帮忙看看,不胜感激! @月只蓝 发自小木虫IOS客户端 |
» 猜你喜欢
什么时候开奖?
已经有13人回复
散金币祈福
已经有94人回复
物理学I论文润色/翻译怎么收费?
已经有55人回复
青基已中
已经有1人回复
散金币祈福
已经有186人回复
求助,如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
求助中医院临床中医中标的青年标书学习,谢谢
已经有2人回复










回复此楼