2015-03-22 3 views
0

У меня есть следующее уравнение:изменение переменной в Дифференциал уравнения с клена

eq1 := 2*diff(phi(r),r)/r+diff(phi(r),$(r,2)) + psi(r)^2*phi(r);

Я хочу, чтобы изменить независимую переменную г 1/г, то есть изменение переменной R = 1/г. Как написать это уравнение с использованием новой переменной

ответ

0

Вы хотите, чтобы команда dchange пакета PDEtools.

eq1 := 2*diff(phi(r),r)/r+diff(phi(r),$(r,2)) + psi(r)^2*phi(r); 

       /d  \         
       2 |--- phi(r)| /2  \     
       \ dr  / | d   |   2  
     eq1 := -------------- + |---- phi(r)| + psi(r) phi(r) 
        r   | 2  |     
           \ dr  /    

PDEtools[dchange](r=1/z, eq1); 

         /     /2  \\     
    3/d  \ 2 | /d  \ 2 | d   ||   2  
-2 z |--- phi(z)| - z |-2 z |--- phi(z)| - z |---- phi(z)|| + psi(z) phi(z) 
     \ dz  / |  \ dz  / | 2  ||     
         \      \ dz  //     

simplify(%); 
       /2  \      
        | d   | 4   2  
        |---- phi(z)| z + psi(z) phi(z) 
        | 2  |      
        \ dz  /     
+0

Он работает, спасибо – user55944

+0

Вы знаете, как изменить производную Diff (п, г) к производной Diff (ф (г), г) – user55944

+0

Параметр 'Команда value' изменяет инертную форму' Diff (f (r), r) 'в вызов активного' diff (f (r), r) '. Например, 'значение (Diff (f (r), r))'. Ваш комментарий странный, потому что вы использовали только 'f' и' f (r) ', и неясно, хотите ли вы использовать их как' f (r) '. Использование обоих из них, в частности, могло бы иметь смысл, если бы вы имели в виду «D (f) (r)» вместо этого, которые могут быть преобразованы с помощью «convert (D (f) (r), diff)». Вы должны были спросить об этом в совершенно новом Вопросе, а не в качестве комментария здесь. – acer

Смежные вопросы