Курсовая работа: Решение параболических уравнений

{

clrscr();

FILE *myfile;

FILE *plotter;

float a[120] [120];

float b[120] [120];

float u[120] [120];

float T = 0.05;

float l = 1;

float h = 0.1;

float tau = 0.01;

int n, i, j, k;

float s = pow (h, 2) / tau;

n = ceil (l / h);

for (i = 0; i <= 119; i++)

{

for (j = 1; j <= 119; j++)

{

u[i] [j] = 0;

a[i] [j] = 0;

b[i] [j] = 0;

}

}

for (i = 0; i <= n; i++)

{

u[i] [0] = phi (i * h);

}

for (j = 0; j <= floor (T /tau); j++)

{