Математические модели технических объектов

//res = 3; //test1

//res = (2*xRight+3*y+5); //test2

res = 2*xRight*xRight+3*y*y*y+4; //test3

return res;

}

 

public static double V3(double x)

{

double res;

//res = 3; // test 1

//res = kappa3*(2*x+3*yLeft+5)-(x+2*yLeft+1)*3; //test2

res = kappa3*(2*x*x+3*yLeft*yLeft*yLeft+4) - (1+x+2*yLeft*yLeft)*(9*yLeft*yLeft); //test3

return res;

}

 

public static double V4(double x)

{

double res;

//res = 3; // test1

 

//res =2*x+3*yRight+5;// test2

res = 2*x*x+3*yRight*yRight*yRight+4; // test3;

return res;

}

// наш  функционал - правая часть уравнения

public static double fun(double x, double y)

{

double res;

//res = 0; // test 1;

//res = -10; //test2

res = -(24*x*x + 18*x*y + 22*y + 72*y*y*y + 4); // test3

return res;

}

 

}

}

 

 

 

 

 

 

 

Приложение 2

Анализ порядка  аппроксимации уравнения и граничных  условий, вывод выражения для  главного члена погрешности аппроксимации:

 

Анализ порядка аппроксимации  уравнения и вывод выражения  для главного члена погрешности аппроксимации.

 

Причем, коэффициент при h¹ равен нулю – разность левой и правой частей ( ), а коэффициент при h³ – главный член погрешности аппроксимации

В качестве порядка погрешности аппроксимации уравнения предлагается брать разность показателей степени h при главном члене и члене,  при котором стоит разность правой и левой частей.

В общем случае: получаем порядок аппроксимации p = (m+k)-m.

Здессь p = 2.

При граничных условиях:

 

Коэффициент при h⁰ равен нулю, т.к. он является разностью левой и правой частей гран. условия, а при h¹ стоит разность левой и правой частей уравнения.

В общем случае:   получаем порядок аппроксимации p=(l+k)-l.

2005