i + +)
{[0] [i]=fi (X [i]); [0] [0] [i]=X [i]; [0] [0] [i]=u [0] [i];
}
/ / V0 (s, t, x)=fi (x) (int i=0; i
{(int j=0; j
{(int k=0; k
{[i] [j] [k]=fi (X [k]);
/ / cout << V [i] [j] [k] << endl;
}
}
} (int j=0; j
{(int k=0; k
{[j] [k]=fi (X [k]);
}
} (int i=0; i
{(int j=0; j
{(int k=0; k
{[i] [j] [k]=fi1 (X [k]);
/ / cout << V [i] [j] [k] << endl;
}
}
} out («fi.txt»); << nx << endl; (int i=0; i
{(int j=0; j
{<< u [j] [i] << « t»;
} << endl;
}.close();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 (int k=0; k
{(int j=0; j
{(int i=0; i
{(i <= j)
{= 0; (int s1=i; s1
{+=V [s1] [j] [k] * dt;
} [i] [j] [k]=X [k]-sumst ;/ / sumst=integral ot s (i) do t (j) ot Vn - 1=sum (V * dt)
} (i> j)
{= 0; (int s1=i; s1> j; s1 -)
{+=V [s1] [j] [k] * dt;
} [i] [j] [k]=X [k] + sumst ;/ / sumst=integral ot s (i) do t (j) ot Vn - 1=sum (V * dt)
}=fabs (eta1 [i] [j] [k]-eta [i] [j] [k]);
if (difeta1> difeta) difeta=difeta1;// Знаходження максимальної різниці між наближеннями eta ??p>
}
}
} << «difeta =» << difeta << endl; (int k=0; k
{(int j=0; j
{(int i=0; i
{(int s2=0; s2
{+=(2 / ((1 + V [s2] [j] [k]) * (1 + V [s2] [j] [k])) - 2 / (( 1 + V [s2] [j] [k]) * (1 + V [s2] [j] [k]))) * dt;
} [i] [j] [k]=exp (-eta1 [0] [j] [k] * eta1 [0] [j] [k]) * (- 2 * ( eta1 [0] [j] [k])) * sumst1; s=0; (int s2=0; s2
{s +=(2 / (1 + V [s2] [j] [k]) - 2 / (1 + V [s2] [j] [k])) * dt;
}=0; (int s2=0; s2
{+=(W [s2] [j] [k] * W [s2] [j] [k]) * dt;
} [i] [j] [k]=exp (-eta1 [0] [j] [k] * eta1 [0] [j] [k]) + sum0s + sumst2;
//V1[i][j][k]=fi(X[k])*exp(-sum0s);=fabs(V1[i][j][k]-V[i][j][k]);<<«ijk»<<«t»<<i<<«t»<<j<<&la...