すいません、手に負えないのでお願いします。

① ラプラス方程式を解くプログラムを作れ。
② ①を用いて長方形領域 0≦x≦2、0≦y≦1の上でラプラス方程式を解け、
　　ただし境界条件（周上の関数値）はデータとして読み込むものとする。
③ 適当なxおよびyの値の組{(x1,y1),(x2,y2),...,(xm,ym)}（例えば，y=ax^3+bx^2+cx+d+ノイズ）
　　に対してべき多項式モデルf=Xcを当てはめた結果をグラフ表示せよ。なお、べき次数ｎ
　　およびデータ点数mを様々に変化させよ。