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) endfunction
q=1; y=URQ(q,eta(q),zeta(q)) Rgds, Rafael
_______________________________________________ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users