No BSD License
function f=shooting(N) global N1 N2 N3 N1=N(1); N2=N(2); [t,y]=ode45('MaxwellStefan',[0 1],[0.319 0.528]); f(1)=y(end,1); f(2)=y(end,2);
Contact us at files@mathworks.com