Помогите дописать программу в Scilab

[Petr88]

Задача
Известно значение средней кривизны поверхности



Требуется найти значение функции f во всех узлах сетки и построить поверхность с помощью Scilab.

Я написал программу, только вот дальше не получается...

n=10; m=10;
x=linspace(0,1,n);
y=linspace(0,1,m);
u=zeros(n,m);
p=u;
q=u;
r=u;
s=u;
t=u;
v=u;
H=2*ones(n,m);
Hij=zeros(n*m,1);
Hij(:)=H(:,:);
for i=1:m
  u(1,i)=2*i;
  u(n,i)=u(1,i);
end
for i=1:n
  u(i,1)=4*i;
  u(i,m)=u(i,1);
end
C=splin2d(x,y,u);
[XX,YY] = ndgrid(x,y);
[u,p,q,r,s,t]=interp2d(XX,YY,x,y,C,"C0");


A=(ones(n,m)+q.*q)./((sqrt(ones(n,m)+p.*p+q.*q)).^3);
B=-2*p.*q./((sqrt(ones(n,m)+p.*p+q.*q)).^3);
Z=(ones(n,m)+p.*p)./((sqrt(ones(n,m)+p.*p+q.*q)).^3);

ZIJ=zeros(n*m,n*m);

Там дальше надо цикл написать, который значения вычислит и все... только не получается...

[AKM]

Пожалуйста, исправьте написание формул в соответствии с Правилами.
Здесь рассказано, как набирать формулы (здесь подробнее).
Используйте кнопку для редактирования своего сообщения.

Тема перемещена из "Помогите решить (М)" в карантин.
Как исправите - пишите сюда, чтобы тему вернули.