задача с прямым течением времени, где ошибка?

[eaita]

исправьте пжл где моя ошибка,не могу понять где Вот код:

#include "stdafx.h"
#include <conio.h>
#include <math.h>
#include <iostream>
using namespace std;

const int B=-1;
float nu(float x) {
return 5*exp(B*x);
}

float m1(float t) {
return 5*exp(B*t);
}

float m2(float t) {
return 5*exp(B*(1+t));
}

void main() {
int i,j;
float x[10], t[10], y[10][10], tao, h=0.01, tr[10][10], ro[10][10], k[10][10], f[10][10], a[10], b[10], c[10], d[10], r, alfa[10], beta[10];
a[10]=0, b[0]=0, c[0]=1, c[10]=1;
t[0]=0;
tao=0.01;
t[1]=tao;
r=tao/(h*h);

for(i=0; i<10; i++){
x[i]=i*h;
y[i][0]=nu(x[i]);

for(j=0; j<10; j++){
t[j+1]=t[j]+tao;
y[0][j]=m1(t[j]);
y[10][j]=m2(t[j]);

ro[i][j+1]=1+x[i]+t[j+1];
k[i][j+1]=x[i]*x[i]+t[j+1]*t[j+1]+5;
tr[i][j]=5*exp(B*(x[i]+t[j]));
f[i][j+1]=5*B*exp(B*(x[i]+t[j+1]))*(ro[i][j+1]-2*x[i]-B*k[i][j+1]);
}}

for(i=0; i<10; i++){
	for(j=0; j<10; j++){
d[0]=m1(t[j+1]);
d[10]=m2(t[j+1]);

a[i]=r*k[i-1/2][j+1];
b[i]=r*k[i+1/2][j+1];
c[i]=ro[i][j+1]+(k[i-1/2][j+1]+k[i+1/2][j+1])*r;
d[i]=tao*f[i][j+1]+ro[i][j+1]*y[i][j];
	}
}

/*for(i=1; i<=9; i++) {
	alfa[0]=b[0]/c[0];
	beta[0]=d[0]/c[0];
	alfa[i]=b[i]/(c[i]-a[i]*alfa[i-1]);
	beta[i]=(d[i]+beta[i-1]*a[i])/(c[i]-a[i]*alfa[i-1]);
	for(j=0; j<9; j++){
		beta[10]=y[10][j+1];
	}}*/

for(i=9; i>=1; i--){
	alfa[0]=b[0]/c[0];
	alfa[i]=b[i]/(c[i]-a[i]*alfa[i-1]);
	beta[i]=(d[i]+beta[i-1]*a[i])/(c[i]-a[i]*alfa[i-1]);
	for(j=0; j<9; j++){
		beta[0]=m1(t[j+1]);
		beta[10]=m2(t[j+1]);
		y[i][j+1]=alfa[i]*y[i+1][j+1]+beta[i];
		printf("i=%1d, j+1=%1d, pogr=%0.5f",i,j+1,abs(y[i][j+1]-tr[i][j+1]));
		printf("\n");
	}
}
getch();
}

[Abstraction]

float x[10], t[10], y[10][10], tao, h=0.01, tr[10][10], ro[10][10], k[10][10], f[10][10], a[10], b[10], c[10], d[10], r, alfa[10], beta[10];

Мама... Я это не буду править, я этого боюсь. Не рассматривали для начала какие-нибудь способы, гм, некоторого уменьшения количества переменных?

[eaita]

пыталась сократить как могла, алгоритм дали, попорядку перевела на с++)

а на float обычно внимание не обращаю, походу дела добавляю и получаеца так(

[evg_m]

a[i]=r*k[i-1/2][j+1];
b[i]=r*k[i+1/2][j+1];
c[i]=ro[i][j+1]+(k[i-1/2][j+1]+k[i+1/2][j+1])*r;
что за индекс такой 1/2?

алгоритм дали, попорядку перевела на с++)

может будет лучше если мы увидим алгоритм.
поскольку

y[0][j]=m1(t[j]);
 y[10][j]=m2(t[j]);

вместе с определением m1 и m2 и желание симметрии наводит на мысль о том что надо было

y[i][0]=5 * exp(B*t[0]);
y[i][10] =5 * exp(B*t[10]);

и сделать это до цикла.

[eaita]

Постановка задачи и входные данные