non linear differential equation

1 view (last 30 days)
mohan
mohan on 30 Nov 2012
please help me how to solve
dq/dt + q = a1 + b1*u(t) + 1/u(t)^2

Answers (1)

Azzi Abdelmalek
Azzi Abdelmalek on 30 Nov 2012
Edited: Azzi Abdelmalek on 30 Nov 2012
function dq=myeq(t,q,t1,u1,a1,b1)
u=interp1(t1,u1,t)
dq=-q+a1+b1+1/u^2;
%Then call myeq
a1=1;
b1=2;
q0=0;
t1=1:0.1:10; % vector time
u1=t1.^2 % your vector u
[t q]=ode45(@(t,q) myeq(t,q,t1,u1,a1,b1),t1,q0)
plot(t,q)

Categories

Find more on Programming in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!