Clear[lös,a,b,c,a0,b0,c0,n1,n2,k1,k2];
prec=16;<<Graphics`;

(*  Folgereaktionen		A -> B -> C		

 	Reaktionsordnungen				ni:
	Geschwindigkeitskonstanten		ki:
	Anfangskonzentrationen			i0:	*)
 	
	n1=0.0; n2=1.0;

	k1=0.4; k2=0.8;

	a0=1.0; b0=0.0; c0=0.0;


eqn ={a'[t]== -k1*(a[t])^n1,	a[0]==a0,
	  b'[t]== +k1*(a[t])^n1
	          -k2*(b[t])^n2,	b[0]==b0,
	  c'[t]==  k2*(b[t])^n2,	c[0]==c0};

st=0; en=100; (* start,ende der Reaktion  *)  
       
lös=NDSolve[eqn,{a[t],b[t],c[t]},{t,st,en},
				MaxSteps->1000,
				WorkingPrecision->prec];

a=a[t]/.lös[[1]];
b=b[t]/.lös[[1]];
c=c[t]/.lös[[1]];

Print["Die Rechnung ist beendet"];

Plot[{a,b,c},{t,st,5},
	 PlotStyle->{Red,Green,Blue},
	 PlotRange->{-0.6,(a0+b0+c0+0.1)},
	 PlotLabel->{"k1=",k1,"k2=",k2,
					 "a0=",a0,"b0=",b0},
	 AxesLabel->{"s","mol/l"}];