Discussion:GTP

Autres discussions [liste]

NB : GTP est une abréviation aussi utilisée par

corrige moi ce code fin je puisse tracer l courbes de resonnance[modifier le code]

Implicit none

Integer,parameter::n=100000

double precision:: f,omega, gama, mu, al1, al2, al3,eta ,al4, al5,al6,ha

double precision:: fs, fp, omegamax, omega0, amax, a, db,lamda

double precision::a11,b11,c11,x1,x2,pmt,delta

open(1,file='data1.dat')

!!!!!! Enter the parameters !!!!!!!!!!!!!!

! eta=0.; gama=3.; f=0.5; lamda=0.5; mu=0.1;omega=2.

eta=-0.30; gama=3.; f=0.5; lamda=0.5; mu=0.1;omega=2.

! !!!!!!!!!!!!!!!!!!! Initialisation of the parameters !!!!!!!!!!!!!!!!!!!!!!!!!!!!

omega0=0.00; a=0.; omegamax=2.; amax=1.; db=0.00001

!Pause

do pmt=omega0, omegamax, 0.00001

omega=pmt

a11=((9*eta**2)/16)-((9*omega**9+gama**2)/(2)+(f**2)/64)

b11=-(((3*omega**2*eta)/(2))+((3*gama*omega**4*lamda)/(2))+((f**2)/(2))+(3*eta)/(2))

c11=((omega**4)+1+(lamda**2*omega**2)-(2*omega**2)-(f**2))

if(a11.ne.0.) then

delta=b11**2-4*a11*c11

if(delta .ge.0.) then

x1=(-b11+sqrt(delta))/(2*a11)

x2=(-b11-sqrt(delta))/(2*a11)

if (x1.ge.0.) write(1,*) pmt,sqrt(x1)

if (x2.ge.0.) write(1,*) pmt,sqrt(x2)

end if

else if (b11.ne.0.)then

x1=-c11/b11

if(x1.ge.0.) write(1,*) pmt,sqrt(x1)

end if