|
本帖最后由 kyra 于 2020-3-10 09:30 编辑
integer,parameter :: N = floor(L/ls)
ls = 0.545, L= 100
所以 N = 183
real,parameter :: xx(4)=(/x00, x00+2*lt, x00+2*lc, x00+2*(lc+lt)/)
xx的大小是4
do I1 = 1,NM
do J1 = 1,N
Zpp(I1,J1)= ((2/(mr*L))**0.5)*sin(I1*pi*xx(J1)/L)!//此处J1循环到N,当循环到5时,xx越界(因为xx大小为4)
end do
end do |
|