second order partial differential element in an ODE

Hi, Does anyone know who to write a script for the below equation?
z' = 1/L*[Ω Pctr *Vs *cos(Ωt)−(R +2∂L/∂u*v)*z−(1/c+ (∂^2 L)/(∂u^2 )v^2+ ∂L/∂u *v')*y].
The main concern is the second order partial differential element.
Thanks a lot.

4 Comments

Additionally, L, u, v, y and z are all functions of t.
Davis sent me a copy of the paper, but I do not know how to deal with the diff(L(t),u(t)) part.
There must be something more relating z to some of u, v and y (L?)? Or some of them to make this a second order ODE.
There is; Davis has shown only one of the equations here.
L(t) has a known exponential formula, but u(t) is defined through another differential equation.
There are two particular difficulties:
1) I do not know how to differentiate a function with respect to another function, diff(L(t),u(t))
2) because u(t) is being worked out numerically by an ode*() call (as u' is one of the equations), the formula for u(t) is not initially known, which makes differentiating with respect to u(t) suspect.
Ok. I would (simplistically) start trying to make functions for dL/du and d^2L/du^2 - which should work if the OP just want to solve the ODE for one given L (or if L, dL/du, d^2L/du^2 is called with additional parameters). This might work if L doesn't show up in the LHS of the ODE - which would turn it into something else than an ODE...

Sign in to comment.

Answers (0)

Tags

Asked:

on 3 Sep 2011

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!