Помогите решить

← →
Profi © (2004-10-20 19:41) [0]

Есть дифур второго порядка:

q""=-2*E*q"-W*W*q-X0

E и W - const
X0 - задан таблично, причем шаг не известен и есть ли он вообще? Беру 1.
Пробывал Рунге-Кутом, изменяя уравнение до вида
Z"=-2*E*Z-W*W*Z*Z/2-X0
т.к. q0=0 Z0=0

Но, поскольку E порядка десятых, а W - сотых, то шаге на восьмом происходит переполненеи Z, даже если брать Extended.

Может у кого есть какие сображения, поделитесь!

← →
Profi © (2004-10-20 22:41) [1]

Рассмотрю любые предложения :)

← →
Maxi_L (2004-10-21 00:00) [2]

сори, но я только начинаю, пока ничем помочь не могу

← →
Profi © (2004-10-21 17:13) [3]

Либо орешек не прост, либо...
В общем я кое-что сделал, но возник вопрос: X0 - табличное значение и их 1024, какой брать шаг? 1 или 1024?

← →
Alx2 © (2004-10-21 17:21) [4]

q(x) = -X0/(W^2)+_C1*exp((-E+sqrt((E-W)*(E+W)))*x)+_C2*exp(-(E+sqrt((E-W)*(E+W)))*x); (с) Maple

← →
Profi © (2004-10-21 17:26) [5]

Alx2 © (21.10.04 17:21) [4]
:) В том то и дело, что найти надо q"", а не прсто q!

q" = _C1*(-E+sqrt((E-W)*(E+W)))^2*exp((-E+sqrt((E-W)*(E+W)))*x)+_C2*(-E-sqrt((E-W)*(E+W)))^2*exp(-(E+sqrt((E-W)*(E+W)))*x) (c) Maple

Alx2 © (21.10.04 18:07) [6]
Продифференцировать я и сам мог твою формулу! Ты ТЗ читал вообще? Нужен цыкл, а так ни C1 ни C2 мы не знаем!

Помогите решить Найти похожие ветки