global A, global B, global epsr
B=1.36e-03;
epsr=1.0e-06;

A=0.2377;
[z0 tz0 k0 zz0]=rk4stb("fsurf2",[0;0;0;0;0;0],0,0.015,80000,2);
A=0.303;
[z1 tz1 k1 zz1]=rk4stb("fsurf2",[0;0;0;0;0;0],0,0.015,80000,4);
A=0.3514;
[z2 tz2 k2 zz2]=rk4stb("fsurf2",[0;0;0;0;0;0],0,0.015,80000,6);
A=0.39;
[z3 tz3 k3 zz3]=rk4stb("fsurf2",[0;0;0;0;0;0],0,0.015,80000,8);

off0=max(z0(2,:));
off1=max(z1(2,:));
off2=max(z2(2,:));
off3=max(z3(2,:));
gg=plot(z0(1,:),z0(2,:)-off0,z1(1,:),z1(2,:)-off1,z2(1,:),z2(2,:)-off2,z3(1,:),z3(2,:)-off3);

surfaceenergy0=zz0(11)
gravenergy0=B*(zz0(12)-zz0(8)*zz0(10))