molecule m;
float xyz[4000],fret;
file file1;
file1=fopen("outx", "w" );

m = getpdb ("1h22_10.pdb");
readparm(m,"1h22_10.prmtop");
mm_options( "cut=999., ntpr=10, nsnb=99999, diel=C, gb=0 , dielc=4.0 ,wcons=200.0,maxarcdot=1500" );
setxyz_from_mol( m, NULL, xyz );
mme_init( m, "::H*", "::!@H=", xyz, file1);

// conjugate gradient minimization
conjgrad(xyz, 3*m.natoms, fret, mme, 0.00001, 0.000001, 1000 );

// Newton-Raphson minimization\fP
mm_options( "ntpr=1" );
newton( xyz, 3*m.natoms, fret, mme, mme2, 0.00000001, 0.0, 20 );

setmol_from_xyz( m,NULL, xyz );
putpdb("1h22_10_md.pdb",m);

// get the normal modes:
nmode( xyz, 3*m.natoms, mme2, 10, 0, 0.0, 0.0, 0);

fclose(file1);