Математический форум Math Help Planet
Обсуждение и решение задач по математике, физике, химии, экономике Теоретический раздел |
Часовой пояс: UTC + 3 часа [ Летнее время ] |
новый онлайн-сервис число, сумма и дата прописью |
|
Часовой пояс: UTC + 3 часа [ Летнее время ] |
Страница 1 из 1 |
[ Сообщений: 5 ] |
|
Автор | Сообщение | |
---|---|---|
Alena164 |
|
|
q[1]:P*diff(f(r),r[1])=R[1]; q[2]:P*diff(f(r),r[2])=R[2]; функция следующая: f(r):=B1*r[1]^(b1)*r[2]^(b2); пробовала через солв, но выдает ошибку: solve([q[1],q[2]],[r[1],r[2]]); algsys: tried and failed to reduce system to a polynomial in one variable; give up. -- an error. To debug this try: debugmode(true); затем решила попробовать метод ньютона, но опять выдал ошибку mnewton([q[1],q[2]],[r[1],r[2]],[1,1]); Maxima encountered a Lisp error: Error in LAMBDA [or a callee]: LAMBDA [or a callee] requires more than zero arguments. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. не подскажите, в чем проблема?) Заранее спасибо) |
||
Вернуться к началу | ||
erjoma |
|
|
Вернуться к началу | ||
Alena164 |
|
|
Спасибо, но максима оказывается только через ньютона такие вещи решает, оказывается надо было просто сначала все неизвестные задать, это ведь численный метод) и все заработало)
у меня возник другой вопрос.. мне нужно сделать это в общем виде, чтобы я задала функцию и входные параметры и на выходе получила мои r[i] я сделала так: B1:1;u1:(1/2); u2:(1/3);P:10;R[1]:1; R[2]:2; f(r):= B1*r[1]^(u1)*r[2]^(u2); for k:1 thru n step 1 do (q[k]:P*diff(f(r),r[k])=R[k]); mnewton([q[1],q[2]],[r[1],r[2]],[5,5]); но теперь мне нужно в mnewton записать не конкретно q1 и q2, а q[k] и r[k], чтобы я могла изменить только функцию, но не менять саму "программу". но к сожалению максима запоминает только последнее значение, в данном случае 2. и даже если я прошу сделать список из q[k] ( а список в максиме это все что в []) он все равно запоминает послед значение.. теперь я думаю, что в максиме этого сделать нельзя, может есть какие то мысли?) |
||
Вернуться к началу | ||
WWWorms |
|
||
Помогите решить !!!! Мой е-mail Senya-ac@mail.ru
|
|||
Вернуться к началу | |||
Avgust |
|
|
Главное - получить верный ответ с использованием Maple.
Я Вашу систему решал так. Из первого уравнения получим: [math]y=-\frac 79 x[/math] Подставим во второе и получим [math]{\frac {32}{9}}\,x+15\,\ln \left( \frac 29\,x \right) -5=0[/math] Составим итерацию Ньютона: [math]x_{n+1}=x_n-\frac{{\frac {32}{9}}\,x_n+15\,\ln \left( \frac 29\,x_n \right) -5}{\frac {32}{9}+\frac{15}{x_n}}[/math] Теперь составим программу в Maple при начальном [math]x_0=2[/math]: restart; x := 2: for n to 6 do x1 := evalf(x-((32/9)*x+15*ln((2/9)*x)-5)/(32/9+15/x)): print(n, x1, -(7/9)*x1); x := x1 end do: Точность в последней строке больше, чем требовал препод. Можно было бы ограничиться и четырьмя циклами. В принципе и самое начало моего поста можно было бы запустить в прогу. При таком подходе начальное [math]y_0=-0.5[/math] оказывается лишним. |
||
Вернуться к началу | ||
[ Сообщений: 5 ] |
Часовой пояс: UTC + 3 часа [ Летнее время ] |
Кто сейчас на конференции |
Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 4 |
Вы не можете начинать темы Вы не можете отвечать на сообщения Вы не можете редактировать свои сообщения Вы не можете удалять свои сообщения Вы не можете добавлять вложения |