1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
| for j:= 1 to 60 do begin
for k:=1 to 48 do begin
for i:=1 to nbratome do begin
r[i]:=sqrt(abs(sqr(xgh[j]-xb[i])+sqr(ygh[k]-yb[i]))); {jai du rajouter abs car il refusait de faire la racine du nb !}
if r[i]<rayonVdW[i] then
affichpixel[j,k]:=false { affichpixel est un tableau de booleen qui me servira par la suite }
else
affichpixel[j,k]:=true;
V:=charge[i]/r[i]; { la division par zero pose pb a ce nivo la ! }
Vtot[j,k]:=Vtot[j,k]+V; { vtot est le potentiel sur lensemble de ma molecule dapres les pixels de coordonées j,k de mon ecran avec une certaine resolution predefinie ! }
end;
end;
end;
for j:=1 to 60 do begin
for k:=1 to 48 do begin
writeln(Vtot[j,k]); { je voulais verifier le calcul de mes potentiels , mais bon rien ne saffiche !! LOL }
end;
end; |
Partager