Большие рассчеты

#1 Опубликовано: 26.05.2010 22:45:10

0 сообщений из 5 понравились пользователям.

Группа: User

Сижу считаю физику. Сдавать завтра. Не считает! Описывать словами долго, файлик в аттаче.
Программка классная! Спасибо Вам за нее =)
клик

#2 Опубликовано: 26.05.2010 23:04:39

Andrey Ivashov

2270 сообщений из 3734 понравились пользователям.

Группа: Super Administrator

Чтобы получить символьный ответ dhN1 нажмите Ctrl+. (или кнопку "->" на панели Арифметика) после ввода имени переменной.

#3 Опубликовано: 26.05.2010 23:17:34

Odin

0 сообщений из 5 понравились пользователям.

Группа: User

Пардоньте(С) не указал что именно не считает.
цикл, в котором считается dhN2 в качестве результата возвращает пустой массив из 0. Пробовал выбирать различные типы оптимизации в ответ получал: значение либо не поместилось в Int32, либо деление на 0(?!) при различных типах оптимизации.
аналогично для dhN1.
попробовал подставить вручную значения в формулу (дифференциал h(...) по N1 ) и всё нормально вычислилось (переменная q).
мне надо именно вычисление этих массивов. уже самые различные варианты перепробовал.

#4 Опубликовано: 26.05.2010 23:22:38

Andrey Ivashov

2270 сообщений из 3734 понравились пользователям.

Группа: Super Administrator

Несколько смущает задание функции
[math]qwe(N1;N2;R1;R2)←-{w*Bv1*N2^{5/4}*ln(R2/R2)^4}/{S*A0*v1^4*(N1^{1/4}-N2^{1/4})^5}[/math]
С таким аргументом логарифма не удивительно, что функция всегда возвращает ноль

#5 Опубликовано: 26.05.2010 23:28:26

Odin

0 сообщений из 5 понравились пользователям.

Группа: User

ладно. это дифференциал по N1 функции h(...).
пытался даже напрямую таким образом: dhN1:=d/dN1 h(...). аналогичная ситуация.
а если в эту функцию вручную подставить значения все, то она возвращает нормальное значение -- видно по переменной q

#6 Опубликовано: 26.05.2010 23:38:10

Andrey Ivashov

2270 сообщений из 3734 понравились пользователям.

Группа: Super Administrator

И ещё раз, настоятельно прошу обратить внимание на аргумент при натуральном логарифме [math]ln(R2/R2)[/math] - ответ всегда будет равен нулю! У вас ошибка в формуле.

#7 Опубликовано: 26.05.2010 23:45:45

Odin

0 сообщений из 5 понравились пользователям.

Группа: User

хорошо. тогда получается, что производная по N1 функции h(...) будет всегда равна нулю при R2>R1 ? в таком случае почему значение функции не равно нулю если подставлять туда все значения вручную?

#8 Опубликовано: 26.05.2010 23:51:50

Andrey Ivashov

2270 сообщений из 3734 понравились пользователям.

Группа: Super Administrator

Ну почему Вы не хотите проверить то, что сами написали?

Почувствуйте разницу (посмотрите на аргумент логарифма!): вот производная:
[math]diff(h(N1;N2;R1;R2);N1)—-{w*Bv1*N2^{5/4}*ln(R2/R1)^4}/{S*A0*v1^4*(nthroot(N1;4)-nthroot(N2;4))^5}[/math]

а вот Ваша функция:
[math]qwe(N1;N2;R1;R2)←-{w*Bv1*N2^{5/4}*ln(R2/R2)^4}/{S*A0*v1^4*(N1^{1/4}-N2^{1/4})^5}[/math]

Программа всё вычислила верно, а вот Вы зачем-то после копирования ответа его модифицировали и при этом допустили ошибку.

#9 Опубликовано: 26.05.2010 23:53:52

Odin

0 сообщений из 5 понравились пользователям.

Группа: User

эх блин... извините.. время позднее. мозги не работают толком