The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: help with ode45 function
Replies: 0  

Advanced Search

Back to Topic List Back to Topic List  

Posts: 3
Registered: 6/15/11
help with ode45 function
Posted: Feb 9, 2013 10:06 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

I'm trying to solve a first-order differential equation using 'ode45', and it is taking way too long. I'm wondering if I am doing something wrong.

Basically, the ODE I'm trying to solve is V' = constant*f(V)

So, I have a fun1.m file:

function f = fun1(t,Vm)

gk = 0.415;
gcl = 0.582;
gna = 0.01;
Ek = -74.7*10^-3;
Ena = 54.2*10^-3;
Ecl = -65.8*10^-3;
Cm = 1e-6;

f = (-1/Cm)*(gk*(Vm-Ek)+gna*(Vm-Ena)+gcl*(Vm-Ecl));

Then, I am doing:

clear all
Vm_initial = -60*10^-3;
[t1 f1] = ode23('fun1',[0 10],Vm_initial);

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.