линейная двумерная интерполяция

[армфайт]

написать подпрограмму линейной двумерной интерполяции функци z(x,y)алгоритм метода даны три точки (x1,y1,z1)(x1,y2,z2)(x2,y1,z3)прово дим через них плоскость z=ax+by+c даны точкf(x,y)в качестве интерполированого значения в этой точке берем z=ax+by+c входные параметры (x1,y1,z1,x1,y2,z2,x2,y1,z3)точки(x 1,y1,z1)(x1,y2,z2)и(x2,y1,z3) x,y-точки(x,y) выходные параметры z-ответ(=ax+by+c)si-индикатор ошибки=0 найдено z в режиме интерполяции точка(x,y)лежит внутри треугольника образованого точками (x1,y1)(x1,y2)(x2,y1)=1 найденого z в режиме экстрополяции. точка (x,y)лежит вне треугольника образованого точками (x1,y1)(x1,y2) и (x2.y1)=2 эти три точки не определяют плоскость однозначно. даже если примерно можете решить буду очень признателен