H=input('pleasse input the value of H:');
M=input('pleasse input the value of M:');
disp(H);
disp(M);
if H>0.5 && H<1
Q=0
else if H>0 && H<0.5
U=0;
for u=1:Inf
U=U+u^(H-1.5)
end;
Q=(0.5-H)*U;
end;
end;
disp(Q);
G(u:1)=randn(M-1);
for u=(t-M)
t-1)
T(1,
=(t-u)^(H-1.5);
end;
G(t)=randn(1);
f=(H-0.5)*(T*G(u))+Q*G(t);
plot(t,f,'b*');
grid on;