Ловлю Matlab, Кто и что умеет... Опции

[vikusin]

По причине своих огромных долгов перед наукой плтно занимаюсь МатЛабом. Вот тут задалась целью собрать коллекцию его жучков, чтобы разуверить граждан студентов во всесильности данной программы ( а то становится прямо болезнью -все решит MatLab или Maple, а сами даже задачу проанализировать перед тем, как ее в машину заталкивать не удасуживаются).

Раньше спокойно некорректные задачи подсовывала и чаами любовалась на мучения. А новые версии уже на такой мелочи на сыпятся. Думаю, что делать.

Если кто-нибудь налетал, подскажите где копать.

[Elena]

Ага, я уже успела поковырять на оф. сайте MathWorks. Придется вычислять в мапле и передавать в матлаб. Лень мне писать программу по вычислению дробных производных...

[lehtym]

Ага, я уже успела поковырять на оф. сайте MathWorks. Придется вычислять в мапле и передавать в матлаб. Лень мне писать программу по вычислению дробных производных...

Можно в такие выси не лезть, а решать численно дифференциальное уравнение Лежандра.
Тут есть два пути:
- вычесть логарифмическую особенность в единице (если есть) и поставить краевую задачу для регулярной части.
- решать начальную задачу: значения Лежандров в нуле явно выписываются через Гамма-функцию.

Первый путь я сам использовал - получается не очень быстро, но в целом терпимо.
На втором пути надо грамотно выбрать численную схему: при комплексных индексах задача жесткая получается.

[Elena]

Можно в такие выси не лезть, а решать численно дифференциальное уравнение Лежандра.
Тут есть два пути:
- вычесть логарифмическую особенность в единице (если есть) и поставить краевую задачу для регулярной части.
- решать начальную задачу: значения Лежандров в нуле явно выписываются через Гамма-функцию.

Первый путь я сам использовал - получается не очень быстро, но в целом терпимо.
На втором пути надо грамотно выбрать численную схему: при комплексных индексах задача жесткая получается.

Если б это была основная задача, я бы, может, и решала её с помощью первого пути, а т.к. эта задача была вспомогательной, то решалась она наиболее быстрым (для решающего) способом.

Кстати, еще один глюк Матлаба:
Если написать в условии (if или while) что-нибудь в духе a>0.0001, забыв по рассеянности, что число a комплексное, и не поставив модуль, то Матлаб даже не вякнет, что что-то в условии не так, а спокойно сравнит только действительную часть этого числа.