global mg h l L0;
mg = 0.1; h=0.1; l=0.1; L0=0.1;
function e=eta(q)
e= mg./(2*q)
endfunction
function z=zeta(q)
z= asinh(mg.*h./(2*q.*l.*sinh(eta(q))))+eta(q);
endfunction
function y=URQ(q, e, z)
y= mg.*L0./(q.*l)-sinh(2*e-z)-sinh(z)
endfunctionq=1; y=URQ(q,eta(q),zeta(q)) Rgds, Rafael
_______________________________________________ users mailing list [email protected] http://lists.scilab.org/mailman/listinfo/users
