Вычисление кривизны Гаусса - MathCAD - Ответ 1378335

function gc = gcurvature(x,y,z)
syms u v
xu   =   diff(x,u,1);
xuv  = diff(xu,v,1);
xv   =  diff(x,v,1);
xuu  =   diff(x,u,2);
xvv  = diff(x,v,2);
 
yu   =   diff(y,u,1);
yv   =   diff(y,v,1);
yuv  =   diff(yu,v,1);
yvv  = diff(y,v,2);
yuu  = diff(y,u,2);
 
zu   =   diff(z,u,1);
zv   = diff(z,v,1);
zuv  =   diff(zu,v,1);
zuu  =   diff(z,u,2);
zvv  = diff(z,v,2);
 
 for i=1:(size(z,1))
     for j=1:(size(z,2))
        Xu          =   [xu(i,j) yu(i,j) zu(i,j)];
        Xv          =   [xv yv zv];
        Xuu         =   [xuu yuu zuu];
        Xuv         =   [xuv yuv zuv];
        Xvv         =   [xvv yvv zvv];
        E           =   dot(Xu,Xu);
        F           =   dot(Xu,Xv);
        G           =   dot(Xv,Xv);
        m           =   cross(Xu,Xv);
        n           =   m/sqrt(sum(m.*m));
        L           =   dot(Xuu,n);
        M           =   dot(Xuv,n);
        N           =   dot(Xvv,n);
       
        %[U, V] = meshgrid(linspace(0, 2*pi, n), linspace(0, 2*pi, m));
        gc = ((L*N)-M^2)/((E*G)-F^2);
      %surfl(x, y, z );
     % hold on
     % plot(gc)
  % axis([-5 5 -5 5 -5 5])
   %axis off
   %colormap pink
  % hold off
      end
 end