//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop
#include<iostream.h>
#include<conio.h>
#include<math.h>
//---------------------------------------------------------------------------
double h=0.1;
double e=0.01;
double f1(double x)
{
      return 0;
}
double f2(double x)
{
        return 50*x*(1-x);
}
double f3(double y)
{
        return 50*y*(1-y*y);
}
double f4(double x)
{
        return 50*x*(1-x);
}
 static void Derihle(double x, double y, double s, double f)
{
        s=1;//double x,y,s=1,f;
        int a=(1.0/h)+1;
        int b=(1.0/h)+1;
        double u[100][100];
        for (int i = 0; i < a; i++)
            {

                for (int j = 0; j < b; j++)
                {
                    u[i][j] = 0;
                }
            }
for (int i = 0; i < a; i++)
            {
                x = i * h;
                y = i * h;
                u[i][0] = f2(x);
                u[i][b - 1] = f4(x);
                u[0][i] = f1(y);
                u[a - 1][i] = f3(y);
            }
            int k = 0;
            while (s > e)
            {
                k++;
                for (int i = 1; i < a - 1; i++)
                {
                    for (int j = 1; j < b - 1; j++)
                    {
                        f = u[i][j];
                        u[i][j] = (u[i + 1][j] + u[i - 1][j] + u[i][j + 1] + u[i][j - 1]) / 4;
                        s = abs(f - u[i][j]);
                    }
                }
            }

                for (int i = 0; i < a; i++)
            {
                for (int j = 0; j < b; j++)
                {
                        cout<<u[j][i]<<"\t";
                }
                cout<<"\n";
            }

        }

        void main()
        {
            Derihle(x,y,s,f);
            getch();
        }




//---------------------------------------------------------------------------