syms u v W P real;
x=cosh(u)*cos(v);
y=cosh(u)*sin(v);
z=u;
Ru=[diff(x,u) diff(y,u) diff(z,u)];
Rv=[diff(x,v) diff(y,v) diff(z,v)];
g11=Ru*Ru'
g12=Ru*Rv'
g22=Rv*Rv'
Ruu=[diff(x,u,2) diff(y,u,2) diff(z,u,2)];
Rvv=[diff(x,v,2) diff(y,v,2) diff(z,v,2)];
Ruv=[diff(diff(x,u),v) diff(diff(y,u),v) diff(diff(z,u),v)];
N=simplify(cross(Ru,Rv));
n=simplify(N/sqrt(N*N'));
Norman=simplify(n*n')
b11=simplify(Ruu*n')
b12=simplify(Ruv*n')
b22=simplify(Rvv*n')
dudvK=simplify(solve((b11*g12-b12*g11)*W^2+(b11*g22-b22*g11)*W+b12*g22-b22*g12,W))
dudvA=simplify(solve(b11*P^2+2*b12*P+b22,P))
[v,u]=ode45('difeq',[0,2*pi],0)
