#include"const.h"
#include<iostream>
#include<cmath>
using namespace std;

double ***junc(int nx, int ny, int p)
{
    int i,j,n,k,ii;
    double hx,hy, ***a,hxr,hyr,xc,yc,xw,yw;
    k=p-1;
    n=nx+ny+1;
    hy=SW/ny;
    hx=WF/nx;
    a= new double**[n];
    for (i=0;i<=n;i++)
    {
        a[i]=new double*[p];
        for (j=0;j<=p;j++)
            a[i][j]=new double[2];

    }


    for (j=0;j<=p;j++)
    {
        a[0][j][0]=0;
        a[0][j][1]=j*(SW-R)/p;
        cout<<"0 "<<j<<'\n';
    }

    for(i=1;i<=nx;i++)
    {
        xc=R/sqrt(SW*SW/(i*i*hx*hx)+1); //x_circle
        yc=SW-sqrt(R*R-xc*xc);//y_circle
//        cout<<yc<<'\n';
        xw=i*hx; //x_wall
        yw=0; //y_wall
        hxr=(xw-xc)/(k+1);
        hyr=(yc-yw)/(k+1);
  //      cout<<"hyr "<<hyr<<'\n';

        for(j=0;j<=p;j++)
        {
            a[i][j][0]=xc+hxr*j;
            a[i][j][1]=yc-hyr*j;

            cout<< i <<' '<<j<<'\n';
        }
    }

    for(i=nx+1;i<n;i++)
    {
        ii=i-nx;
        xc=R/sqrt((SW-ii*hy)*(SW-ii*hy)/(WF*WF)+1); //x_circle
        yc=SW-sqrt(R*R-xc*xc);//y_circle
//        cout<<yc<<'\n';
        xw=WF; //x_wall
        yw=hy*(ii); //y_wall
        hxr=(xw-xc)/(k+1);
        hyr=(yc-yw)/(k+1);
//        cout<<"hyr "<<hyr<<'\n';

        for(j=0;j<=p;j++)
        {
            a[i][j][0]=xc+hxr*j;
            a[i][j][1]=yc-hyr*j;
            cout<< i <<' '<<j<<'\n';
        }
    }

    return a;
}