function Simulation1
date=load('D:\plog.txt');
x=date;
for t=1:481
y(t)=x(t,:);
end
s(0)=0;
y(0)=0;
for t=1:481
s(t)=y(t)-c-a*y(t-1)-b*s(t-1);
end
d=g^(2);
f=1/(sqrt(2*pi*d))*exp{-s(t)^(2)/(2*d)};
h=cumprod(f);
L=log(h);
A=[a,b,c,d];
d1=diff(L,A);
fun=infine('d1');
[a,b,c,d]=ode23(fun,0.00,0.00,0.00,1);
end