Matlab - Ecuaciones diferenciales

clear all;clc

syms u(t) %declara simbólica la variable 

intervalo=0:10; %defines tu intervalo de 0 a 10  de uno en uno 

uo=[10, 20, 50, 60, 70, 80 , 100]; %aquí las diferentes condiciones iniciales 

hold on %aquí para que las gráficas queden en una misma figura

for k=1:length(uo) % para las primera condiciones inicial hasta la última 

sol(t)=dsolve(diff(u,t)==(0.7-0.01*u)*u,u(0)==uo(k)); %resolver ecuación 

f=matlabFunction(sol); %convertir a función anónima 

y=f(intervalo);

if isequal([1,1],size(y))

    y=repmat(y,1,length(intervalo));

end

plot(intervalo,y) %representar

end