@skatersgitt

Помогите, пожалуйста, разобраться со сложившейся ситуацией. Есть проблемный кусок программы, ума не приложу, что маткаду не нравится. При определенных значениях (L<=1) всё считается. При бОльших значениях маткад выдаёт ошибку. Притом, если менять значения ρs или/и μs, расчёт продолжается. Пересчитывал это уравнение вручную, пересчитывал в вольфраме: при сколь угодно маленьких/больших значениях ls/L корни для скорости (W) есть. Полагаю, проблема в логике этой программы. Помогите, пожалуйста, переписать, господа. Неделю бьюсь...
P.S. Это вырезка из большой программы. Есть еще похожее место с подобной проблемой. С ним постараюсь сам разобраться, коль скоро добрые люди дадут понять, в чём же загвоздка.
Заранее спасибо!

Вложения

problem.zip (57.6 Кб, 9 просмотров)

0

@Symon

1. Сравнение z=1 внутри цикла for выглядит несуразно. Оно выполняется только при первом заходе в цикл. Лучше этот случай
вынести из цикла и при выполнении условия в цикл уже не заходить.
2. Имеется лишняя строка внутри цикла, где есть сравнение =. Это нигде не используется.
3. Использование root здесь ни к чему. Решается уравнение простым извлечением корня.
4. Исправление всего сказанного вас не спасет.
Основная ошибка в строке с Cd. Это значение можно вычислить только в случае, если Re попадает в два интервала. Если нет, маткад попадает в прострацию, что и происходит при некоторых значениях параметров.
5. Не ясен смысл цикла for, где W обновляется несколько раз (по сути это рекурсия). Чем больше z, тем больше повторений. Понять это можно только если знать, какую задачу решаете.
6. Реальные константы в программе выглядят некрасиво. Лучше заменить их буквенными обозначениями перед программой.
7. Для вычисления Сd лучше определить функцию (вне программы).
8. На рабочем листе неплохо навести порядок. Все данные разбросаны как попало.

0