Trying to solve motion equation using ODE45

Hello

I am trying to solve this motion equation , but i get no answer and MATLAB stays busy without converging . Is there a simplier way to solve this , Thank you

Knowing that Rini = R2.

function dydt = vdp1(t,y)
ep0=0.00000000000885;   %Vacuum Permittivity
epa=80;  %Area Permittivity               
epp=1000; %Particle Permittivity
fcm=(epp-epa)/(epp+2*epa); %Clausius–Mossotti-Factor
d=0.000001; %Particle Diameter
m=5.24*10^(-16); %Particle Mass
v=1000; %Applied Voltage
R2=3000; %External Radius
R1=1500; %Internal radius
vis=0.0013; %Dynamic Viscosity
dydt = [y(2);(4*pi*ep0*epa*fcm*((d/2)^3)*((v^2)/((log(R2/R1))^2))-(1/((y(1)-R2)^3))-3*pi*vis*d*y(2))/m];
end


[t,y] = ode45(@vdp1,[0 200],[3000 ; 0]);
plot(t,y(:,1),'-o',t,y(:,2),'-o')
title('Solution of the Equation with ODE45');
xlabel('Time t');
ylabel('Solution y');
legend('Y1 Position','Y2 Velocity')

Answer 1

[darova]

The problem

Try this modifications

[t,y] = ode45(@vdp1,[0 1e-6],[2900 ; 0]);