kaku=30; dt=2*M_PAI/kaku; for(i=0;i