function diffusin(h)
	x=0:h:10;
	n=length(x);
	alea=sqrt(h)*rand(2,n,"normal");
	res=zeros(2,n);
	res(1,1)=0.1;
	res(2,1)=0.1;
	for j=1:(n-1)
		res(2,j+1)=res(2,j)+1*h+sin(res(2,j))*alea(2,j);
	end
	plot2d(x,res(2,:))
	plot2d(x,x)
	a=3.1415*ones(1,n);
	b=6.2830*ones(1,n);
	c=9.4245*ones(1,n);
	plot2d(x,a)
	plot2d(x,b)
	plot2d(x,c)
endfunction