输入矩阵到ode45的功能文件

Question

I have a code(which requires a lot input be given by the user) which will give me anxn matrix(say A), which I have to use to solve a system of ODEs X'=AX. How do I include this matrix A in the function file(.m file) of ode45. If I include the code in the function file as below:

function xp=g_test(t,x);
k=input('something');
A=some manipulation of inputs;
xp=A*x;
end

Matlab asks for input at each timestep(typically my problem has 30k timesteps). So how do I include/ pass the matrix A to the function?

Answer 1

You can create a function that returns a function handle :

function odeFcn = makeODE(a,b,c)  
    A = some_function(a, b, c);
    odeFcn = @(t,x) A*x;
end

Now you can call ode45 with input matrices a, b, c :

outputVector = ode45(makeODE(a,b,c), [t0, t1], x0);

---

The inspiration is taken from gnovice's answer here .