C语言一维和二维混用问题

[show off]

问下我第一张后面那个for循环是对的吗？？运行结果不出错，但是没数据，，

[赵4老师]

请贴代码，不要贴截图。

[show off]

#include "stdio.h"
#include "math.h"
#include "4trmul.c"
#include "6gaus.c"
#define NJ 120
#define N1 10000
main()
{
int i, j, p = 7850;
double dt = 0.00001, t = 1, s = 3, em = 2.0e11, dy2 = 0.100, dy1 = 0.050, g = 9.8, d, nl = 120;
d = s / nl;
double bb = 0.5, rr = 0.25, ss, me, ge, a0, a1, a2, a3, a4, a5, a6, a7, pi = 3.1415926525;
a0 = 1 / (rr*dt*dt);
a1 = bb / (rr*dt);
a2 = 1 / (rr*dt);
a3 = 1 / (rr * 2) - 1;
a4 = bb / rr - 1;
a5 = (bb / rr - 2)*(dt / 2.0);
a6 = dt*(1 - bb);
a7 = bb*dt;
pi = 3.1415926525;
ss = 0.25*pi*(dy2*dy2 - dy1*dy1);
me = p*ss*d / 6.0;
ge = em*ss / d;
printf("%13.3lf%13.8lf%13.3d\n", d, me, p);
double m[NJ][NJ] = { 0.0 }, sk[NJ][NJ] = { 0.0 }, su[NJ] = { 0.0 }, skk[NJ][NJ] = { 0.0 }, ff[NJ] = { 0.0 }, ui[NJ] = { 0.0 }, f, u[NJ][N1] = { 0.0 }, uu[N1] = { 0.0 }, u1[NJ][N1] = { 0.0 }, u2[NJ][N1] = { 0.0 };
for (i = 0; i < NJ; i++)
for (j = 0; j < NJ; j++)
{
if (i == j) sk[i][j] = 2 * ge;
if (i == j - 1 || i - 1 == j) sk[i][j] = -1 * ge;
}
sk[NJ - 1][NJ - 1] = 1 * ge;
for (i = 0; i < NJ; i++)
{
for (j = 0; j < NJ; j++)
{
if (i == j) m[i][j] = 4 * me;
if (i == j - 1 || i - 1 == j) m[i][j] = 1 * me;
}
}
m[NJ - 1][NJ - 1] = 2 * me;
for (i = 0; i < NJ; i++)
{
for (j = 0; j < NJ; j++)
{
skk[i][j] = sk[i][j] + a0*m[i][j];
}
}
f = 200000 * sin(30 * pi*dt*(i + 1));
for (i = 1; i < N1 ; i++)
{
for (j = 0; j < NJ; j++)
{
uu[j] = a0*u[i-1][j] + a2*u1[i-1][j] + a3*u2[i-1][j];
}
trmul(m, uu, 5, ui);
for (j = 0; j < N1; j++)
{
ff[j] = f + ui[j];
}
gaus(skk, ff, 5);
for (j = 0; j < NJ; j++)
{
u[i][j] = ff[j];
u2[i][j] = a0*u[i][j] - a0*u[i-1][j] - a2*u1[i-1][j] - a3*u2[i-1][j];
u1[i][j] = u1[i-1][j] + a6*u2[i][j] + a7*u2[i][j];

}
}
printf("%13.1lf\n", u[i+4]；
}

就这个头文件没问题

[zhouqunhai]

这图看的眼睛疼