Помогите решить задачу и дописать программу в 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);

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

-- Ср май 18, 2011 09:50:10 --

Или может ее как-то по-другому можно решить, иначе, чем я делаю...

[Petr88]

Тут же получается метод конечных разностей?

[AKM]

PAV в сообщении #171140 писал(а):

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