Отправляет email-рассылки с помощью сервиса Sendsay

RFpro.ru: Консультации по информатике

Подписаться Бесплатная «Серебряная» новостная рассылка . Подписчиков 293 RSS
←   Май 2010   →
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31


За последние 60 дней ни разу не выходила


Сайт рассылки: http://rusfaq.ru/issues/8/19/228
Открыта: 17-04-2009

Автор
Калашников О.А.

Статистика

293 подписчиков
0 за неделю
←   Все выпуски   →

RFpro.ru: Информатика


Хостинг портала RFpro.ru:
Московский хостер
Профессиональный ХОСТИНГ на базе Linux x64 и Windows x64

РАССЫЛКИ ПОРТАЛА RFPRO.RU

Чемпионы рейтинга экспертов в этой рассылке

Гаряка Асмик
Статус: Специалист
Рейтинг: 3472
∙ повысить рейтинг » 		lamed
Статус: Бакалавр
Рейтинг: 2516
∙ повысить рейтинг » 		_Ayl_
Статус: Практикант
Рейтинг: 1833
∙ повысить рейтинг »

/ НАУКА И ОБРАЗОВАНИЕ / Точные и естественные науки / Информатика

Номер выпуска:188
Дата выхода:12.05.2010, 20:30
Администратор рассылки:Калашников О.А., Руководитель
Подписчиков / экспертов:154 / 137
Вопросов / ответов:1 / 2

Вопрос № 178272: Здравствуйте, уважаемые эксперты! Мне необходимо написать программу в qbasic, интегрирование ОДУ методом Рунге-Кутты, по данным алгоритмам. http://rfpro.ru/upload/2294 http://rfpro.ru/upload/2295 Индивидуально задание: дифференциальное ур...

Вопрос № 178272:

Здравствуйте, уважаемые эксперты!
Мне необходимо написать программу в qbasic, интегрирование ОДУ методом Рунге-Кутты, по данным алгоритмам.
http://rfpro.ru/upload/2294
http://rfpro.ru/upload/2295
Индивидуально задание: дифференциальное уравнение y’=y (y=y(x)), интервал [0,1], начальные условия x0 = 0, y0 = 1, шаг интегрирования h0 = 0,1.
Вот что у меня получилось:

Код:
DECLARE FUNCTION r! (x() AS SINGLE, y() AS SINGLE, h AS SINGLE, m AS SINGLE)
DECLARE FUNCTION f! (x AS SINGLE, y AS SINGLE)
DIM n AS INTEGER
DIM x0 AS SINGLE, y0 AS SINGLE, h AS SINGLE, m AS SINGLE, b AS SINGLE
CLS
INPUT "Nachalnoe uslovie x0:"; x0
INPUT "Nachalnoe uslovie y0:"; y0
INPUT "Pogreshnost:"; E
INPUT "Sha
 g pechati:"; h0
INPUT "Konec otrezka integrirovaniya [x0;b]:"; b
PRINT
n = FIX((b - x0) / h0) + 1
PRINT "Nachalnoe uslovie x0:"; x0
PRINT "Nachalnoe uslovie y0:"; y0
PRINT "Shag pechati:"; h0
FOR i = 1 TO n
h = h0
m = 1
y = r(x(), y(), h, m)
DO UNTIL ABS(y - y1) < E
y1 = y
h = h / 2
x = x0
y = y0
m = 2 * m
y = r(x, y, h, m)
LOOP
x0 = x0 + h0
y0 = y
PRINT "x0 ="; x0,"y0 ="; y0, "h0
="; h0, "m ="; m
NEXT i
END

FUNCTION f (x AS SINGLE, y AS SINGLE)
f = y(x)
END FUNCTION

FUNCTION r! (x() AS SINGLE, y() AS SINGLE, h AS SINGLE, m AS SINGLE)
FOR j = 1 TO m
k1 = f(x(j), y(j))
k2 = f(x(j) + h / 2, y(j) + h * k1 / 2)
k3 = f(x(j) + h / 2, y(j) + h * k2 / 2)
k4 = f(x(j) + h, y(j) + h * k3)
y(j + 1) = y(j) + h / 6 * (k1 + 2 * k2 + 2 * k3 + k4)
x = x + h
NEXT j
r = y(j + 1)
END FUNCTIO
 N

Программа выдает ошибку, связанную с переменной X. Видимо я что-то не поняла с процедурами-функциями. Помогите, пожалуйста, исправить программу.
В результате должны получиться значения x и y(x) из приведенной таблицы.
http://rfpro.ru/upload/2296
Заранее спасибо!

Отправлен: 07.05.2010, 20:23
Вопрос задал: Ushastik1985, Посетитель
Всего ответов: 2
Страница вопроса »

Отвечает lamed, Бакалавр :
Здравствуйте, Ushastik1985. Ответ в приложении. QBasic
Внесены некоторые исправления.
1. Исправлена передача параметров. Надо считать x[i+1] и x[i] итерациями, а не элементами массива.
2. Исправлена функция f=y.
3. Исправлен вызов функции y = r(x0, y0, h, m) в первом цикле.
Удачи!

Приложение:

Ответ отправил: lamed, Бакалавр
Ответ отправлен: 07.05.2010, 23:01
Номер ответа: 261274

Оценка ответа: 5
Комментарий к оценке:
Спасибо за полный и основательный ответ.

Вам помог ответ? Пожалуйста, поблагодарите эксперта за это!
Как сказать этому эксперту "спасибо"?
  • Отправить SMS #thank 261274 на номер 1151 (Россия) | Еще номера »
  • Отправить WebMoney:

    • Отвечает Лиджи-Гаряев Владимир, Студент :
    Здравствуйте, Ushastik1985.

    В предыдущем ответе, x выбегает из отрезка [0;1]
    Предлагаю свой вариант исправлений.
    Код:
    DECLARE FUNCTION r! (x AS SINGLE, y AS SINGLE, h AS SINGLE, m AS SINGLE)
    DECLARE FUNCTION f! (x AS SINGLE, y AS SINGLE)
    CLS
    INPUT "Nachalnoe uslovie
x0:"; x0
    INPUT "Nachalnoe uslovie y0:"; y0
    INPUT "Pogreshnost:"; E
    INPUT "Shag pechati:"; h0
    INPUT "Konec otrezka integrirovaniya [x0;b]:"; b
    PRINT
    n = FIX((b - x0) / h0) + 1
    PRINT "Nachalnoe uslovie x0:"; x0
    PRINT "Nachalnoe uslovie y0:"; y0
    PRINT "Shag pechati:"; h0
    y1 = y0
    FOR i = 1 TO n
    h = h0
    m = 1
    y = r(x, y, h, m)
    DO UNTIL ABS(y - y1) < E<
 br>y1 = y
    h = h / 2
    x = x0
    y = y0
    m = 2 * m
    y = r(x, y, h, m)
    LOOP
    x0 = x0 + h0
    y0 = y
    PRINT "x ="; x0, "y ="; y0, "h ="; h0, "m ="; m
    NEXT i
    END

    FUNCTION f (x, y)
    f = y
    END FUNCTION

    FUNCTION r! (x, y, h, m)
    FOR j = 1 TO m
    k1 = f(x, y)
    k2 = f(x + h / 2, y + h * k1 / 2)
    k3 = f(x + h / 2, y + h * k2 / 2)
    k4 = f(x + h, y + h * k3)
    y = y + h * (k1 + 2 * k2 + 2 * k3 + k4) / 6
    x = x + h
    NEXT j
    r
= y
    END FUNCTION

    Ответ отправил: Лиджи-Гаряев Владимир, Студент
    Ответ отправлен: 07.05.2010, 23:28
    Номер ответа: 261275

    Оценка ответа: 5

    Вам помог ответ? Пожалуйста, поблагодарите эксперта за это!
    Как сказать этому эксперту "спасибо"?
  • Отправить SMS #thank 261275 на номер 1151 (Россия) | Еще номера »
  • Отправить WebMoney:

    • Оценить выпуск »
    Нам очень важно Ваше мнение об этом выпуске рассылки!

    Задать вопрос экспертам этой рассылки »

    Скажите "спасибо" эксперту, который помог Вам!

    Отправьте СМС-сообщение с тестом #thank НОМЕР_ОТВЕТА
    на короткий номер 1151 (Россия)
    Номер ответа и конкретный текст СМС указан внизу каждого ответа.

    Полный список номеров »

    * Стоимость одного СМС-сообщения от 7.15 руб. и зависит от оператора сотовой связи. (полный список тарифов)
    ** При ошибочном вводе номера ответа или текста #thank услуга считается оказанной, денежные средства не возвращаются.
    *** Сумма выплаты эксперту-автору ответа расчитывается из суммы перечислений на портал от биллинговой компании.

    © 2001-2010, Портал RFpro.ru, Россия
    Авторское право: ООО "Мастер-Эксперт Про"
    Автор: Калашников О.А. | Программирование: Гладенюк А.Г.
    Хостинг: Компания "Московский хостер"
    Версия системы: 2010.6.14 от 03.03.2010
    В избранное