|
Здравствуйте, Валерий Васильевич.
Я неправильно переписал два уравнения - первое и третье и поэтому решал не ту систему. Видимо, эти описки из-за моей болезни.
Система, которую я решал: FindRoot[{(1 + r) (26 pi + 56 w) == 56 pi, (1 + r) (16 pi + 16 w) == 48, (1 + r) (12 pi + 9 w) == 8 pc, 80 w == 5 pc}, {pi, 1.7}, {w, 0.27}, {r, 0.52}, {pc, 4.3}]
Система Стидмана: FindRoot[{(1 + r) (28 pi + 56 w) == 56 pi, (1 + r) (16 pi + 16 w) == 48, (1 + r) (12 pi + 8 w) == 8 pc, 80 w == 5 pc}, {pi, 1.7}, {w, 0.27}, {r, 0.52}, {pc, 4.3}]
Поэтому точное решение по системе Стидмана немного отличается.
РЕШЕНИЕ СИСТЕМЫ СТИДМАНА. {pi -> 1.70416, w -> 0.26849, r -> 0.520797, pc -> 4.29584}
С точностью до 16 знаков после запятой оно правильное.
pi = 1.7041594578792296 w = 0.26849003388254816 r = 0.5207972893961477 pc = 4.295840542120771
Разность между левой и правой частями всех четырёх уравнений равна нулю. Если взять Ваши проверочные формулы, то получим следующее.
КОД ПРОГРАММЫ.
pix = 1.7041594578792296 wx = 0.26849003388254816 rx = 0.5207972893961477 pcx = 4.295840542120771
Print["Proverka"] (1 + rx) (28 pix + 56 wx) - 56 pix (1 + rx) (16 pix + 16 wx) - 48 (1 + rx) (12 pix + 8 wx) - 8 pcx 80 wx - 5 pcx
Print["Proverka po Valeriu"] Print["pi"] (1 + rx) (28 pix + 56 wx)/56 pix
Print["w"] 5 pcx/80 wx
Print["r"] 8 pcx/(12 pix + 8 wx) - 1 rx
Print["pc"] 80 wx/5 pcx ===================
Результат. Все четыре разности левой и правой части равны нулю.
Проверка по Валерию. Результат всех четырёх формул проверки по Валерию совпадает с решением с точностью до 5 знаков после запятой.
Проверочные формулы Валерия дают следующие результаты, которые совпадают с решением. pi = 1.70416 w = 0.26849 r = 0.520797 pc = 4.29584 ==============================================================================
P.S. Рад, что у Вас появилась возможность использовать Mathematica. Программа эта очень обширная и содержит в себе множество возможностей математического изучения разных конструкций. Но не всегда бывает ясно, как на её языке правильно составить код чтобы получить, то, что ищешь. Я до сих пор сталкиваюсь со сложностями при решении ряда сложных задач.
Так, например, в предпоследней главе своей будущей книги я изучаю математические свойства модели, которая воспроизводила бы основные закономерности развития реальной экономики за период с 1400 по 1830 годы. Этот период включает в себя две революции цен и состоит из трёх стадий рыночной экономики. Благодаря исследованиям экономических историков мы знаем сейчас общую динамику цен за весь этот период, имеются также свидетельства о норме прибыли и нормах прибавочной стоимости в 18 и 19 веках. Модель может воспроизвести все эти исторические факты. Но решение уравнений для этой модели "в лоб" не даётся, так как априори не известны начальные (в смысле последующих итераций) значения параметров модели. А при неверном их выборе алгоритм поиска решения даёт либо почти нулевые либо очень большие значения, то есть не сходится. Поэтому решение удалось найти методом многократной прогонки через ряд укороченных систем уравнений.
Это - моё изобретение в области Mathematica. Решать задачу со множеством (в моей реалистичной модели - 8 уравнений) сложных уравнений не напрямую, а сначала разбить её на несколько систем уравнений с меньшим числом уравнений и потом прогонять все эти укороченные системы, уточняя параметры раз за разом, пока не достигнешь решения с нужной точностью.
Задачу эту я решил. Дописываю последнюю главу. Надеюсь, успею.
Доброго Вам здоровья и успехов в исследованиях. С уважением, Григорий.
|