Модераторы: Snowy, MetalFan, bems, Poseidon
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> задать цикл 
:(
    Опции темы
serj07
Дата 27.1.2010, 21:31 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 5
Регистрация: 27.1.2010

Репутация: нет
Всего: нет
помогите разобраться, а то сам запутался.
Есть значение например Х= 0.999542.
задаем значение n=5
в процедуре идет расчет площади S
потом находим погрешность X-S=G
если G>0.0001 то n умножаем на 2 до тех пор пока G<= 0.0001
Всем спасибо!!!
PM MAIL   Вверх
~FoX~
Дата 27.1.2010, 21:36 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


НЕ рыжий!!!
****


Профиль
Группа: Участник Клуба
Сообщений: 2819
Регистрация: 8.10.2003
Где: Зеленоград

Репутация: 5
Всего: 68
Код

while g <=0.0001 do begin
  g := g*2;
end



--------------------
user posted image
…множественность никогда не следует полагать без необходимости…
PM MAIL WWW ICQ Jabber   Вверх
Frees
Дата 27.1.2010, 21:36 (ссылка) | (голосов:1) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Завсегдатай
Сообщений: 2233
Регистрация: 2.12.2005
Где: Екатеринбург

Репутация: 16
Всего: 54
я один ничего не понял?

Добавлено через 1 минуту и 5 секунд
похоже один.


--------------------
Кольцов Виктор Владимирович
PM MAIL ICQ   Вверх
chip_and_dayl
Дата 27.1.2010, 21:44 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 603
Регистрация: 22.4.2007
Где: Украина, Киев

Репутация: 1
Всего: 9
~FoX~
n умножаем на 2. 

Правда не понятно чему присваивать результат

Добавлено через 1 минуту и 4 секунды
Frees
Ну, если взять автора темы и меня, то нас трое smile


--------------------
Фсё будет хорошо!!!
PM MAIL Skype   Вверх
~FoX~
Дата 27.1.2010, 21:54 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


НЕ рыжий!!!
****


Профиль
Группа: Участник Клуба
Сообщений: 2819
Регистрация: 8.10.2003
Где: Зеленоград

Репутация: 5
Всего: 68
Оппа, это я пропустил...
Код

while g <=0.0001 do begin
  n := n*2;
  s := //считаем площадь
  g := X-S;
end



--------------------
user posted image
…множественность никогда не следует полагать без необходимости…
PM MAIL WWW ICQ Jabber   Вверх
serj07
Дата 27.1.2010, 22:12 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 5
Регистрация: 27.1.2010

Репутация: нет
Всего: нет
задание:

Составить программу для вычисления интеграла - I  
с заданной точностью а по формуле трапеций - S1 
где h=1/n. Для достижения необходимой точности число n следует удваивать, пока не будет достигнуто такое n, для которого 
Начальное число узлов n = 5. Требуемая точность q = 0,0001.
Вывести значения n, Ỉn и погрешности q .



интеграл I = 0.9547

Это расчет Трапеции S1


Код

begin
    n:= 5;
    h:=0;
    s1:=0;
    step:=1/n;
      begin
         while h<1 do
            begin
              h:=h+step;
              f1:=sin(h)/h;
              f2:=sin(h+step)/(h+step);
              s1:=s1+(((f1+f2)/2)*step);

            end;

      end;
    end;   



надо вычислить погрешность q<=0.0001

умножая n на 2 в преведущем коде.

этот код вставляю не получается.

Код


   begin
        while q<=0.0001 do
          begin
            q:= I-S1;
            n:=n*2
PM MAIL   Вверх
serj07
Дата 27.1.2010, 23:48 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 5
Регистрация: 27.1.2010

Репутация: нет
Всего: нет
погрешность счетает один раз, n не перемножает до указоной погрешности

Код


procedure TForm1.Button1Click(Sender: TObject);
 var
   q,I, h,f1,f2,s1,step: real;
    n, cnt: word;


begin

    I:= 0.945967;
                       try
    n:=5;
                       EXCEPT
                         MessageDlg('Íåâåðíûå äàííûå', mtError, [mbOK],0);
                             Exit
                             END;

    q:=0;
    h:=0;
    s1:=0;
    cnt:=0;
    step:=1/n;
    chart1.Series[0].Clear;

   begin
        while q<=0.0001 do

       begin

            n:=n*2;

          begin
              while h<1 do

           begin
              h:=h+step;
              f1:=sin(h)/h;
              f2:=sin(h+step)/(h+step);
              Chart1.Series[0].AddXY(h,f1);
              s1:=s1+(((f1+f2)/2)*step);
              inc(cnt)
           end;
         end;
       

 q:= I-s1;
     

  end;
   end;
     begin
        caption:=IntToStr(cnt);
        Label1.Caption:= FloatToStr(s1);
        Label5.Caption:= FloatToStr(n);
        Label4.Caption:= FloatToStr(q);
      end;
 end;
end.
PM MAIL   Вверх
serj07
Дата 28.1.2010, 23:40 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 5
Регистрация: 27.1.2010

Репутация: нет
Всего: нет
немогу разобраться пруга виснит, помогите найти ошибку

Код


procedure TForm1.Button1Click(Sender: TObject);
 var

   q:currency;
   I, h,f1,f2,s1,step: real;
    n, cnt: word;

begin

    I:= 0.004267;  //  это значение интеграла


                     try
    n:= StrToInt(Edit1.Text); // это  количество узлов


                     EXCEPT
                         MessageDlg('Неверные данные', mtError, [mbOK],0);
                        Exit
                      END;

    q:=0;
    h:=0;
    s1:=0;
    cnt:=0;
    step:=1/n;
    chart1.Series[0].Clear;

       begin
              while h<1 do   // здесь ищем второе значение для сравнения

           begin
              h:=h+step;
              f1:=sin(h)/h;
              f2:=sin(h+step)/(h+step);
              Chart1.Series[0].AddXY(h,f1);
              s1:=((f1+f2)/2)*step;
              inc(cnt);
           end;
        end;

        begin

            q:= s1-I;  // находим погрешность

            repeat   // умножаем n на 2 пока недобьемся погрешности 0,001
             n:=n*2;
            until q<=0.001 ;
         end;

      begin
        caption:=IntToStr(cnt);
        Label1.Caption:= FloatToStr(s1);
        Label5.Caption:= FloatToStr(n);
        Label4.Caption:= FloatToStr(q);
      end;
end;
end.
PM MAIL   Вверх
chip_and_dayl
Дата 28.1.2010, 23:56 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 603
Регистрация: 22.4.2007
Где: Украина, Киев

Репутация: 1
Всего: 9
Код

var

   q:currency;
   I, h,f1,f2,s1,step: real;
    n, cnt: word;

begin

    I:= 0.004267;  //  это значение интеграла


                     try
    n:= StrToInt(Edit1.Text); // это  количество узлов


                     EXCEPT
                         MessageDlg('Неверные данные', mtError, [mbOK],0);
                        Exit
                      END;

    q:=0;
    h:=0;
    s1:=0;
    cnt:=0;
    step:=1/n;
    chart1.Series[0].Clear;


              while h<1 do   // здесь ищем второе значение для сравнения

           begin
              repeat
              h:=h+step;
              f1:=sin(h)/h;
              f2:=sin(h+step)/(h+step);
              Chart1.Series[0].AddXY(h,f1);
              s1:=((f1+f2)/2)*step;
              inc(cnt);


            q:= s1-I;  // находим погрешность

            // умножаем n на 2 пока недобьемся погрешности 0,001
             n:=n*2;
            until q<=0.001 ;
         end;
      begin
        caption:=IntToStr(cnt);
        Label1.Caption:= FloatToStr(s1);
        Label5.Caption:= FloatToStr(n);
        Label4.Caption:= FloatToStr(q);
      end;
end;


Добавлено через 27 секунд
Работает, но за корректность данных не ручаюсь


--------------------
Фсё будет хорошо!!!
PM MAIL Skype   Вверх
  
 Ответ в темуСоздание новой темы Создание опроса
Правила форума "Delphi: Для новичков"
SnowyMetalFan
bemsPoseidon
Rrader

Запрещается!

1. Публиковать ссылки на вскрытые компоненты

2. Обсуждать взлом компонентов и делиться вскрытыми компонентами

  • Литературу по Дельфи обсуждаем здесь
  • Действия модераторов можно обсудить здесь
  • С просьбами о написании курсовой, реферата и т.п. обращаться сюда
  • Вопросы по реализации алгоритмов рассматриваются здесь
  • 90% ответов на свои вопросы можно найти в DRKB (Delphi Russian Knowledge Base) - крупнейшем в рунете сборнике материалов по Дельфи

Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, Snowy, MetalFan, bems, Poseidon, Rrader.
 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Delphi: Для новичков | Следующая тема »

Подписаться на тему | Подписка на этот форум | Скачать/Распечатать тему


 




[ Время генерации скрипта: 0.0565 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
 По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.