function [y] = funcion_dist0(u1,u2,u3,u4,u5,u6) %u1=i1 %u2=i2 %u3=i3 %u4=i4 %u5=R %u6=flag imax=15; limsup=20; liminf=10; lim=10; R=sqrt(imax/u5); angvert=((u1+u2)*R^2+limsup+liminf)/2; d=R*sin((angvert-imax)*pi/180); z=R*sin(angvert*pi/180); l=sqrt(R^2-z^2); anghor=(u3+u4)*R^2/2; x=l*sin(anghor*pi/180); yin=-x/tan(anghor*pi/180); %yin=-R*cos(anghor*pi/180); y(1)=d; y(2)=x; y(3)=yin; y(4)=z; y(5)=R; y(6)=anghor;