Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- public class testing_programm {
- static double epsilon = Math.pow(10,-5);
- static double lyambda = 1/5;
- public static void main(String[] args) {
- // ебучий первый пункт
- double a = 0;
- double b = 1;
- double n = 10;
- double h = (b - a) / (n - 1);
- double[][] A = new double[10][10];
- A[0][0] = A[9][9] = h / 2;
- for (int j = 1; j <= n - 2; j++) {
- A[j][j] = h;
- }
- double[][] s = new double[10][10];
- double[] F = new double[10];
- for (int j = 0; j < n; j++) {
- s[j][j] = a + j * h;
- double x = a + (j + 1) * h;
- F[j] = (x * x - 1);
- }
- double[][] K = new double[10][10];
- for(int i = 0; i < n; i++) {
- for(int j = 0; j < n; j++) {
- K[i][j] = 1 / (10 - (a + (i + 1) * h) * (a + (j + 1) * h));
- }
- }
- double[][] A_f = new double[10][10];
- for (int i = 0; i < n; i++) {
- for (int j = 0; j < n; j++) {
- A_f[i][j] = 0 - lyambda * A[j][j] * K[i][j];
- if(i == j) {
- A_f[i][j] += 1;
- }
- }
- }
- inversion(A_f,(int)n);
- double[] U_x = new double[10];
- double sum = 0;
- for(int i = 0; i < n; i++) {
- for (int j = 0; j < n; j++) {
- U_x[i] += A_f[i][j] * F[j];
- }
- System.out.println(U_x[i]);
- }
- //ебучий второй пункт
- n = 3;
- double[] X_i = {0,0.5,1};
- }
- static void inversion(double[][] A, int N)
- {
- double temp;
- double [][] E = new double [N][N];
- for (int i = 0; i < N; i++)
- for (int j = 0; j < N; j++)
- {
- E[i][j] = 0f;
- if (i == j)
- E[i][j] = 1f;
- }
- for (int k = 0; k < N; k++)
- {
- temp = A[k][k];
- for (int j = 0; j < N; j++)
- {
- A[k][j] /= temp;
- E[k][j] /= temp;
- }
- for (int i = k + 1; i < N; i++)
- {
- temp = A[i][k];
- for (int j = 0; j < N; j++)
- {
- A[i][j] -= A[k][j] * temp;
- E[i][j] -= E[k][j] * temp;
- }
- }
- }
- for (int k = N - 1; k > 0; k--)
- {
- for (int i = k - 1; i >= 0; i--)
- {
- temp = A[i][k];
- for (int j = 0; j < N; j++)
- {
- A[i][j] -= A[k][j] * temp;
- E[i][j] -= E[k][j] * temp;
- }
- }
- }
- for (int i = 0; i < N; i++)
- for (int j = 0; j < N; j++)
- A[i][j] = E[i][j];
- }
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement