Контрольная работа: Сглаженная поверхность для границы трех атомов в пространстве

Zris2=((-bt/2)+sqrt(DISC))/at;

Zris1=((-bt/2)-sqrt(DISC))/at;

Xris1=ABD-AC*Zris1;

Xris2=ABD-AC*Zris2;

YY[schit]=gr1;

ZZ[schit]=Zris2;

XX[schit]=Xris2;

Rad[schit]=(sqrt((XX[schit]-Xvr)*(XX[schit]-Xvr)+(ZZ[schit]-Zvr)*(ZZ[schit]-Zvr)));

x1=XX[schit]-Xvr;x2=Rad[schit];y1=ZZ[schit]-Zvr;y2=0;

ugol[schit]=acos((x1*x2+y1*y2)/((sqrt(x1*x1+y1*y1)))/(sqrt(x2*x2+y2*y2)))*180/M_PI;

schit++;

{if (mon==0) {xp=Xris2;yp=gr1;mon++;}

else {line(xp+40,800-yp,40+Xris2,800-gr1);mon=0;}

circle(Xris2+40,800-gr1,2);

} gr1=gr1+1;

} }

}

}

double rk=0.0,pc;

if (Xkos1>rk) {rk=Xkos1;pc=Ykos1;}

if (Xkos1>rk) {rk=Xkos2;pc=Ykos2;}

if (Xkos1>rk) {rk=Xkos3;pc=Ykos3;}

setcolor(15);

setlinestyle(1,0,1);

line(Xc+40,800-Yc,x[1]+40,800-y[1]);

line(Xc+40,800-Yc,x[2]+40,800-y[2]);

line(Xc+40,800-Yc,x[3]+40,800-y[3]);

setlinestyle(0,0,1);

}