subroutine doScatter
!! 
integer,parameter::N = 100
real(wp),dimension(N)::x,y,z
call random_number(x)
call random_number(y)
z = sqrt(x**2+y**2)
call figure()
call subplot(2,2,1)
call xylim([0.0_wp,1.0_wp],[0.0_wp,1.0_wp])
call scatter(x,y)
call ticks()
call labels('x','y','')
call subplot(2,2,2)
call xylim([0.0_wp,1.0_wp],[0.0_wp,1.0_wp])
call scatter(x,y,c=z)
call ticks()
call labels('x','y','')
call subplot(2,2,3)
call xylim([0.0_wp,1.0_wp],[0.0_wp,1.0_wp])
call scatter(x,y,s=(4.0_wp*z+1.0_wp),markColor='blue')
call ticks()
call labels('x','y','')
call subplot(2,2,4)
call xylim([0.0_wp,1.0_wp],[0.0_wp,1.0_wp])
call scatter(x,y,c=z,s=(4.0_wp*z+1.0_wp))
call ticks()
call labels('x','y','')
end subroutine doScatter