Форум программистов
 

Восстановите пароль или Зарегистрируйтесь на форуме, о проблемах и с заказом рекламы пишите сюда - alarforum@yandex.ru, проверяйте папку спам!

Вернуться   Форум программистов > Delphi программирование > Общие вопросы Delphi
Регистрация

Восстановить пароль

Купить рекламу на форуме - 42 тыс руб за месяц

Ответ
 
Опции темы Поиск в этой теме
Старый 20.07.2012, 17:00   #1
EvilHero
Пользователь
 
Регистрация: 19.11.2010
Сообщений: 12
По умолчанию Решение уравнения методом половинного деления и постойка его графика

Уважаемые форумчане ! Помогите пожалуйста разобраться в своих ошибках. Имеется уравнение (x^3)(x-a)^(1/3)-b=0, a и b задаются пользователем, как и точность вычисления корней e, а так же вводятся интервал для постойки графика и нахождения корней в этом интервале. Задача частично у меня решена, но график строится только если я задаю левую границу интервала больше, чем a.
вот собственно код:
Код:
unit Unit1;

interface

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

type
  TForm1 = class(TForm)
    Label1: TLabel;
    Edit1: TEdit;
    Label2: TLabel;
    Edit2: TEdit;
    Label3: TLabel;
    Edit3: TEdit;
    Label4: TLabel;
    Label5: TLabel;
    Edit4: TEdit;
    Edit5: TEdit;
    Button1: TButton;
    Label6: TLabel;
    Label7: TLabel;
    Chart1: TChart;
    Series1: TLineSeries;
    Button2: TButton;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation

{$R *.dfm}

procedure TForm1.Button1Click(Sender: TObject);
function f(x:real):real;

begin
f:=exp(3*ln(x))*(exp((1/3)*ln(x-(strtoint(Edit1.Text)))))-(strtoint(Edit2.Text))
end;
var z,a,b,c,e:real;

begin
e:=(strtoint(edit3.Text));
a:=(strtoint(edit4.Text));
b:=(strtoint(edit5.Text));
if f(a)*f(b)<0 then
begin
Label6.Caption:='Числа введены неверно!';
exit;
end;
repeat
c:=(a+b)/2;
if f(a)*f(c)<0 then
b:=c
else a:=c;
until abs(b-a)<e;
z:=(a+b)/2;
Label6.Caption:=floattostr(z);
end;

procedure TForm1.Button2Click(Sender: TObject);
function f(x:real):real;
begin
f:=exp(3*ln(x))*(exp((1/3)*ln(x-(strtoint(Edit1.Text)))))-(strtoint(Edit2.Text))
end;
var
y:integer;
k:real;
begin
Series1.Clear;
for y:= -100 to 100 do
begin
k:=f(y);
Series1.AddXY(y,k,'',clRed);
end;
end;

end.
EvilHero вне форума Ответить с цитированием
Ответ


Купить рекламу на форуме - 42 тыс руб за месяц



Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
решение методом половинного деления MeVer1c Паскаль, Turbo Pascal, PascalABC.NET 0 01.03.2012 14:21
Решение задачи методом половинного деления Pascal Oolen'ka Помощь студентам 0 26.02.2012 20:16
Нахождение корней уравнения методом половинного деления в ИСР Turbo Pascal 72nikola Фриланс 2 14.06.2011 16:06
Численные методы. Решение нелинейного уравнения методом половинного деления. gree Помощь студентам 1 11.11.2009 18:36
Решение уравнения методом деления отрезка пополам. Методом секущей. Panda196 Паскаль, Turbo Pascal, PascalABC.NET 3 25.11.2008 09:06