Planteo

> restart;

> u:=L/n;

> with(linalg):

Warning, new definition for norm

Warning, new definition for trace

> n:=5;

> L:=5;u;

>

> h:=array(0..n,0..n,[]);

> hx[i,j]:=array(1..n-1,1..n-1,[]);hy[i,j]:=array(1..n-1,1..n-1,[]);hxx[i,j]:=array(1..n-1,1..n-1,[]);

> hyy[i,j]:=array(1..n-1,1..n-1,[]);hxy[i,j]:=array(1..n-1,1..n-1,[]);

generación de expresiones diferenciales discretizadas

> for i from 1 by 1 to n-1 do

> for j from 1 by 1 to n-1 do

> hx[i,j]:=(h[i+1,j]-h[i-1,j])/(2*u);

> hy[i,j]:=(h[i,j+1]-h[i,j-1])/(2*u);

> hxx[i,j]:=(h[i+1,j]+h[i-1,j]-2*h[i,j])/(u^2);

> hyy[i,j]:=(h[i,j+1]+h[i,j-1]-2*h[i,j])/(u^2);

> hxy[i,j]:=((h[i+1,j+1]-h[i-1,j+1])-(h[i+1,j-1]-h[i-1,j-1]))/(4*u^2);

> term[i,j]:=(1+hy[i,j]^2)*hxx[i,j]+(1+hx[i,j]^2)*hyy[i,j]-2*hx[i,j]*hy[i,j]*hxy[i,j];

> eq[i,j]:=term[i,j]=0;print(eq[i,j]);

> poli[i,j]:=term[i,j]:

> od:od: