r=dsolve('Da=-k(1)*a','Db=k(2)*(((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2))^2-(k(3)+k(4))*b-k(5)*b*((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2)+k(6)*k(5)/(k(6)+k(7))*b*((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2)-(k(9)+k(11)+k(13)+k(16))*b^2+k(10)*d+k(12)*e+k(14)*f','Dc=k(7)*k(5)/(k(6)+k(7))*b*((k(5)^2*k(7)^2/(k(6)+k(7))^2*b^2+4*k(2)*(k(1)*a+k(3)*b))^0.5-k(5)*k(6)/(k(6)+k(7)))/2/k(2)-k(8)*c','Dd=k(9)*b^2-k(10)*d','De=k(11)*b^2-k(12)*e','Df=k(13)*b^2-(k(14)+k(15))*f','Dg=k(15)*f','Dh=k(16)*b^2','Di=k(4)*b+k(8)*c','a(0)=100','b(0)=0','c(0)=0','d(0)=0','e(0)=0','f(0)=0','g(0)=0','h(0)=0','i(0)=0','t')
??? Error using ==> mupadmex
Error in MuPAD command: Unknowns must be functions [ode::new]
Error in ==> sym.sym>sym.mupadmexnout at 2018
out = mupadmex(fcn,args{:});
Error in ==> dsolve>mupadDsolve at 190
[var_list,R] = mupadmexnout('symobj::dsolve',sys,x,ignoreConstraints);
Error in ==> dsolve at 97
[R,vars] = mupadDsolve(ignoreConstraints,varargin{1:narg});