FlexPDE (Проблемы с расчет. обл.)Срочно

[Sintanial]

Ребят очень срочно помогите пожалуйста.
Я уже создавал подобную тем в ветке математики...но увы без кода программы слишком не понятно.
И так моя проблема
Мне необходимо нарисовать область так что бы :угол отмеченный красным был такой же как углы отмеченные синим- при этом можно изменять всё что угодно но так что бы угол отмеченный черным был всегда 0 градусов ( то бишь горизонтален- как сейчас на рисунке). И еще одна вещь....надо изменять всё так что бы не поменялся параметр $S$, то бишь площадь в моём случае(это показывается на графике когда запускаешь программу.


Мучаюсь уже неделю....... я конечно понимаю что разбираться в чужом коде это гиблое дело, но другого варианта как вам объяснить мою проблему не знаю
Нажав на надпись "Оффтопик" появится код программы(Просто не хотел что бы много занимало место на странице....а кроме как в оформить ввиде оффтопика больше не знаю как )

(Оффтоп)

title 'Staged Geometry'
VARIABLES {Peremennie}
C(threshold=0.1)
coordinates{Koordinatnai setka v kotoroi nahodits'a oblast'}
  ycylinder('R','z')
 select
     stages=50
 definitions{Constanti and Zavisimosti}
 k=stage
 d=0.2
 C_0=1
 deltad=d/10
 d_t=d+deltad*k
 theta_1=2*Pi/3
 theta_2=2*Pi/3
 deltaR=d/10
 R_t=deltaR*k
 R_tx=R_t*sin(theta_1)
 R_tc=d_t+0.75*R_t*cos(theta_1)
 R_ty= d_t-R_t*(cos(theta_1)+1.5)
 R_m=d_t+R_t*cos(theta_2){d_t/cos(theta_2)/2}
 l=1
 deltal=-1/8*(8*cos(theta_2)^2*deltad*l-2*d*deltad*Pi+4*d*deltad*theta_2-deltad^2*k*Pi+
 2*deltad^2*k*theta_2+2*sin(2*theta_2)*d*deltad+sin(2*theta_2)*deltad^2*k-2*Pi*deltaR^2*k*cos(theta_2)^2-
 4*deltaR^2*k*cos(theta_2)^2*theta_1+4*deltaR^2*k*sin(2*theta_1)*cos(theta_2)^2)/(cos(theta_2)^2*(d+deltad*k))
 l_t=l+deltal*k
 S = d_t*l_t-1/8*(d+deltad*k)^2*(Pi-2*theta_2)/cos(theta_2)^2+1/8*(d+deltad*k)^2*sin(2*theta_2)/cos(theta_2)^2
-1/2*Pi*deltaR^2*k^2+1/2*deltaR^2*k^2*(1/2*Pi-theta_1)+1/2*deltaR^2*k^2*sin(2*theta_1);
 
EQUATIONS {Uravnenie}
(1/R)*dR(C)+dRR(C)+dZZ(C)=0
 
 
 boundaries{Raschetnai oblast' and granichnie uslovia }
 
   region 1
     start"outer"(l_t,d_t) Natural(C)=0
     line to (R_tx,d_t) value(C)=C_0
     ARC (radius=-R_tx) to (0.,R_ty) Natural(C)=0
     line to (0,0) Natural(C)=0
     line to (l_t,0) Value(C)=C_0 ARC(radius=-R_m)
     TO FINISH
 
  monitors
       grid(R,z)
 
  plots
       grid(R,z)  report(AREA_integral(1))
For cycle=1 contour(C) {grafik raschetnoi oblasti C}
 
  histories
     history (l_t) versus k
 
     history (S) versus k
  end