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

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

Yris1=AD-CB*Zris1;

Yris2=AD-CB*Zris2;

YY[schit]=Yris2;

ZZ[schit]=Zris2;

XX[schit]=gr1;

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 (Yris2>pol1 && Yris2<pol2)

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

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

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

} gr1=gr1+1;

}

}

}

if (smen==2)

{

ABD=(-stat2*gr1)/stat1-stat4/stat1;

AC=stat3/stat1;

ABX=ABD-Xc;

at=(AC*AC+1);

bt=(-2)*(ABX*AC+Zc);

ct=(ABX*ABX+Zc*Zc+(gr1-Yc)*(gr1-Yc)-Rc*Rc);

DISC=(bt/2)*(bt/2)-at*ct;

if (DISC>=0)

{