Есть прогресс.
f(x) = sin(x): http://s2.ipicture.ru/uploads/20130525/mTX3NE1K.png
f(x) = x: http://s2.ipicture.ru/uploads/20130525/01C5cQhu.png
На интервале 8 разбиений, но так как вертикальные черты прорисовываются и от точек для параболы, выглядит будто их много.
Собственно 2 вопроса:
1. Это ведь не верно? Будь верно, на f(x) = x, мы бы увидели не сходимость, а искривления параболы? Или я что-то не понимаю?
2. Как убрать вертикальные черты от точек для параболы, которые задаются во вложенном цикле на каждом отрезке разбиения?
Если не верно, подскажите в чём и я продолжу разбираться. Код ниже.
Код:
| Код |
procedure TForm1.simpson;
var
i,j:integer;
x1,x2,x3:real;
x,y,y1,y2,y3:real;
a,b,c:real;
begin
I_simpson:=0;
SeriesSimpson.Clear;
if N mod 2 <> 0 then
EditSimpson.Text:=trim(EditN.text)+'%2!=0'
else
begin
for i:=0 to N-1 do // Проходим по всем разбиениям
begin
i_simpson:=i_simpson + ((points[i].xR- points[i].xL)*(funcs[i].fL+4*funcs[i].fM+funcs[i].fR))/6; // Считаю.
SeriesSimpson.AddXY(points[i].xL, funcs[i].fL); // Рисую.
x := points[i].xL;
x1 := points[i].xL;
x3 := points[i].xR;
x2 := ( x1 + x3 ) / 2;
y1 := funcs[i].fL;
y3 := funcs[i].fR;
y2 := ( y1 + y3 ) / 2;
//Здесь наверное создавать вложенный цикл и в нём добавляем точки для рисования параболы.
while x < x3 do
begin
a := ( y3 - ( x3 * ( y2 - y1 ) + x2 * y1 - x1 * y2 ) / ( x2 - x1 ) ) / ( x3 * ( x3 - x1 - x2 ) + x1 * x2 );
b := ( ( y2 - y1 ) / ( x2 - x1 ) ) - a*( x1 + x2 );
c := ( ( x2*y1 - x1*y2 ) / x2 - x1 ) + a*x1*x2;
x := x+0.25;
y := a*sqr(x) + b*x + c;
SeriesSimpson.AddXY( x, Fx(x) );
end;
end;
SeriesSimpson.AddXY(points[N].xL, funcs[N].fL);
EditSimpson.Text:=floattostrf(i_simpson, ffGeneral, 4, 8);
SeriesSimpson.Active:=false;
end;
end;
|
С нетерпением жду ответа. |
)
И соответственно y2=(y3-y1)/2
Если сдвиг переменной не используется (а исходно предполагалось именно это), то там должна быть сумма.