non linear differential equation

please help me how to solve

 dq/dt + q =  a1 + b1*u(t) + 1/u(t)^2 

Answer 1

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)