n="justify">} out («fi.txt»); << nx << endl; (int i=0; i
{(int j=0; j
{<< u [j] [i] << « t»;// Запис значень функції у файл
} << endl;
}. close ();
/ * визначення змінних, за допомогою яких забезпечується вихід з програми і збіжність послідовних наближень * /
double sumst=0,sumst1=0,sumst2=0,sum0s=0,sum0s1=0,sum0s2=0;e=0.0001;dif=11;dif1=0;difeta=15;difeta1=0;dif2=0;(int n=0; n <5; n + +)
{(int n=0; n <5; n + +)
{= 0 ;/ / max raznost=0 ;/ / max raznost
/ / заповнення масиву eta ??p>
for (int k=0; k
{(int j=0; j
{(int i=0; i
{(i <= j)
{= 0;
/ / обчислення інтеграла
for (int s1=i; s1
{+=V [s1] [j] [k] * dt;
} [i] [j] [k]=X [k]-sumst;
} (i> j)
{= 0; (int s1=i; s1> j; s1 -)
{+=V [s1] [j] [k] * dt;
} [i] [j] [k]=X [k] + sumst;
}=fabs (eta1 [i] [j] [k]-eta [i] [j] [k]); (difeta1> difeta) difeta=difeta1;
}
}
} << «difeta =» << difeta << endl;
/ / заповнення масиву V [i] [j] [k] і W [i] [j] [k] (int k=0; k
{(int j=0; j
{(int i=0; i
{
[i] [j] [k]=eta1 [0] [j] [k] * exp (- 2.0 * i); s=0; (int s2=0; s2
{s +=exp (2.0 * s2) * V1 [s2] [j] [k] * V1 [s2] [j] [k] * dt;
} [i] [j] [k]=exp (- 2.0 * i) * (eta1 [0] [j] [k] * eta1 [0] [j] [k] * 0.5 + sum0s);=fabs (V1 [i] [j] [k]-V [i] [j] [k]); (dif1> dif) dif=dif1;
}
}
} << dif << endl; (fabs (dif2-dif)
{(int j=0; j
{(int i=0; i
{[i] [j] [k]=V1 [i] [j] [k];/* Останнє наближення стає попереднім і цикл продовжується * /
eta [i] [j] [k]=eta1 [i] [j] [k];// Аналогічно
}
}
}
} (int k=0; k
{(int j=0; j
{[j] [k]=W [j] [j] [k];/* Присвоюємо функції u значення масиву W [j] [j] [k] * /
}
}
} out2 («u.txt»); << nx << endl; (int i=0; i
{(int j=0; j
{<< u [j] [i] << « t»;// Запис результатів у файл
} << endl;
}. close ();
/ * запис файлу значень відомої функції u * /
for (int k=0; k