Решение нелинейного уравнения и построение графика

[xMass]

Всем привет
(язык Delphi)
У меня есть 3 вопроса,которые меня немного ставят в тупик.
Прежде расскажу задание начала задание:
"Разработать схему алгоритма и составить программу вычисления таблицы значений функции
Y=a*x^4+b*Ln(x^2+b^2)
известно XN (начальное значение икса), XK (конечное значение икса), DX (шаг, с которым топаем от начального до конечного значения)
a=7.72
b - решение уравнения tg(X)=X"

Ещё, прежде чем спросить выложу текст программы, что у меня получился

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, Grids, StdCtrls, ExtCtrls, TeeProcs, TeEngine, Chart;

type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Edit4: TEdit;
    Label4: TLabel;
    Button1: TButton;
    Button2: TButton;
    StringGrid1: TStringGrid;
    Label5: TLabel;
    Chart1: TChart;
    Button3: TButton;
    Series1: TLineSeries;
    TeeFunction1: TAddTeeFunction;
    Series2: TLineSeries;
    Series3: TLineSeries;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    {procedure FormCreate(Sender: TObject);}
  private
    { Private declarations }
  public
    { Public declarations }
  end;

Type W=array [1..10] of extended;

const a=7.72;
var
  Form1: TForm1;
  X, Y:W;
  K, N: integer;
  XN, XK, DX: extended;
  SK, SX, SY: string;
  b:real;
implementation

{$R *.dfm}

function F_func(z:real):real;
         begin F_Func:=sin(z)/cos(z)-z;
         end;

procedure iwem_b(var b:real{; str:byte});
          const km=150;
                eps=0.001;
          var i:byte;
              fm,fb,m,d:real;
          begin m:=0;
                d:=20;
                for i:=1 to km do
                    begin fm:=F_func(m);
                          b:=(m+d)/2;
                          fb:=F_func(b);
                          if abs(fb)<=eps then exit;
                          if (fm/100)*(fb/100)>0 then m:=b
                                                 else d:=b;
                    end;
                {str:=1}
                form1.label5.Caption:=floatToStr(b);
          end;

procedure TForm1.Button1Click(Sender: TObject);
var k:integer;
begin
XN:=StrToFloat(edit1.text);
XK:=StrToFloat(edit2.text);
DX:=StrToFloat(edit3.text);
{edit4.text:=FloatToStr(b);  }
{Eps:=StrToFloat(edit4.text);}

X[1]:=XN;
Y[1]:=X[1]*X[1]*X[1]*X[1]*a+Ln(X[1]*X[1]+b*b);

N:=trunc((XK-XN)/DX);


For K:=2 to N do
       begin
       X[K]:=X[K-1]+DX;
       Y[K]:=X[K]*X[K]*X[K]*X[K]*a+Ln(X[K]*X[K]+b*b);
       end;

       StringGrid1.Cells[0,0]:='N';
       StringGrid1.Cells[1,0]:='X';
       StringGrid1.Cells[2,0]:='Y';

For k:=1 to N do
       begin
            Str(K, SK);
            form1.StringGrid1.Cells[0,K]:=SK;

            Str(X[K]:10:4,SX);
            form1.StringGrid1.Cells[1,K]:=SX;

            Str(Y[K]:10:4,SY);
            form1.StringGrid1.Cells[2,K]:=SY;

          end;
   end;

procedure TForm1.Button2Click(Sender: TObject);
var k:integer;
begin
edit1.Clear;
edit2.Clear;
edit3.Clear;
edit4.Clear;
chart1.series[0].Clear;
   for K:=1 to N do
     form1.StringGrid1.Cells[0,K]:='';
     form1.StringGrid1.Cells[1,K]:='';
     form1.StringGrid1.Cells[2,K]:='';

end;

end.

Button1 - кнопа вычислить
Button2 - кнопа очистить

Непосредственно вопросы:
1) В label5 значение уравнения tgX=X у меня не выводит, в чем причина? Вообще каким путем можно определить считает ли у меня порграмма константу b?

2) Есть таблица со столбцами 1 - N, сколько получилось иксов, 2 - X, какие иксыб 3 - Y, какие получились игрик. При нажатии кнопки очистить у меня очищается только столбец с N, а с Y, X нет, почему?

3) Как построить график этого уравнения Y=a*x^4+b*Ln(x^2+b^2)
с помощью Chart?

И ещё, есть какая-нибудь удобная программа для написания блок схем ?
Заранее огромное спасибо!

[_Dmitry]

1. Этой строкой вы завершаете процедуру: if abs(fb)<=eps then exit;
поэтому строка form1.label5.Caption:=floatToStr(b) ; не будет выполнена никогда.
Завершать нужно не процедуру, а цикл if abs(fb)<=eps then break;
2. В данном цикле выполняется только первая строка, остальные строки вне цикла

   for K:=1 to N do   
     form1.StringGrid1.Cells[0,K]:='';   
     form1.StringGrid1.Cells[1,K]:='';   
     form1.StringGrid1.Cells[2,K]:='';

замените на

   for K:=1 to N do
   begin   
     form1.StringGrid1.Cells[0,K]:='';   
     form1.StringGrid1.Cells[1,K]:='';   
     form1.StringGrid1.Cells[2,K]:='';
   end;

3. Качайте с сайта клуба книгу Архангельского "100 компонентов Delphi", там все подробно с примерами.

[xMass]

_Dmitry
спасибо за ответ
2 - помогло, 3 - и вправду очень хорошая книга, спасибо
ПО поводу 1 - замена на

procedure iwem_b(var b:real{; str:byte});
          const km=150;
                eps=0.001;
          var i:byte;
              fm,fb,m,d:real;
          begin m:=0;
                d:=20;
                for i:=1 to km do
                    begin fm:=F_func(m);
                          b:=(m+d)/2;
                          fb:=F_func(b);
                          if abs(fb)<=eps then break;
                          if (fm/100)*(fb/100)>0 then m:=b
                                                 else d:=b;
                    form1.label5.Caption:=floatToStr(b);
                    end;
                {str:=1}

          end;

ни к чему не привела
значение констаны b по-прежнему не выводится в 5-й лэйбл

[_Dmitry]

значение констаны b по-прежнему не выводится в 5-й лэйбл

А где вы вызываете процедуру iwem_b, которая и пишет это значение? Я не нашел в вашей программе.

[xMass]

Сообщение от _Dmitry

А где вы вызываете процедуру iwem_b, которая и пишет это значение? Я не нашел в вашей программе.

Все, я нашел ошибку, спасибо огромное за помощь)

[DangerousV]

а что делает вот эта функция?
N:=trunc((XK-XN)/DX);