dt =(2.0d0*dx)/((3.0d0**0.5d0)*pi*cmax) 自動化