Передача аргументов в ode45

Question

Я хочу передать простой аргумент в моей функции ode45. Моя функция заключается в следующем:

function dxdt = state(t,x,vgth,vgval) 
p=1; 
k=10^0.7; 
window1=1-((2*x)-1).^(2*p); 
dxdt=k*(vgval-vgth+1.2)*window1; 
end 

Основной сценарий:

clear all 
step=0.01; 
t = 0:step:10; 
f=2*0.157; 
vg = 5*sin(2*f*t); 
x0=0.01; 
vgth=1.9; 
[t,x] = ode45(@(t,x) state1 (t,x,vgth,vg(t)), t, x0); 
plot(t,x) 

Здесь Когда я исключить параметр vg(t) в аргументе мимоходом, все работает просто отлично. Но в аргументе это не работает с vg(t).

Answer 1

Сборка комментарии к ответу:

function dxdt = state(t,x,vgth,vgval) 
    p=1; 
    k=10^0.7; 
    window1=1-((2*x)-1).^(2*p); 
    dxdt=k*(vgval-vgth+1.2)*window1; 
end 

function vgval = vg(t) 
    f=2*0.157; 
    vgval = 5*sin(2*f*t); 
end 

step=0.01; 
t = 0:step:10; 
x0=0.01; 
vgth=1.9; 
[t,x] = ode45(@(t,x) state(t,x,vgth,vg(t)), t, x0); 

Это следует надеяться работать как задумано.