how to solve ode if value of constant is vector

function [dydt] = diffvar(t,y)

dydt=-k*y+a;

end

end

in this if a is vector how to solve

Answer 1

[Ameer Hamza]

If k and a are constant vectors of equal length the following will work

k = rand(10, 1);
a = rand(10, 1);
IC = zeros(10, 1);
tspan = [0 10];
[t, y] = ode45(@(t,y) diffvar(t,y,k,a), tspan, IC)
plot(t, y)

function [dydt] = diffvar(t,y,k,a)
    dydt=-k.*y+a;
end