function [mil,x,parstart] = milstei2(m,dt)
% Simulates a process on [0,T=m*dt] with drift function  
% and diffusion defined in mildiff.  Initial value ini.
par=inputpar(4);
parstart=par([3 4]);
dw=normrnd(0,sqrt(dt),1,m+20);
%w1 = [0,cumsum(dw)];
ini=-par(1)/par(2);
mil = ini*ones(1,m+20);
for i=1:(m+20)
d=mildiff(mil(i),dt,dw(i),par);
mil(i+1) = mil(i) + d;
end
x=dt*[0,1:m];
mil=mil(20:(m+20))';
%plot(x',mil,'r')
