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

Reply via email to