function t=traycuad(a,b,c,d,e,alfa,fi);
delta=atan(a*sin(fi)./(d-a*cos(fi)));
g=sqrt(a^2+d^2-2*a*d*cos(fi));
mu=acos((b^2+g.^2-c^2)./(2*b*g));
x=a*cos(fi)+e*cos(mu-delta+alfa);
y=a*sin(fi)+e*sin(mu-delta+alfa);
t=[x;y];
A0=[0;0];
A1=[a*cos(fi(1));a*sin(fi(1))];
B0=[d,0];
B1=A1+d*[cos(mu(1)-delta(1));sin(mu(1)-delta(1))];
P=[x(1);y(1)];
figure
plot(A0(1),A0(2),'.');
hold on
plot(A0(1),A0(2),'o');
plot(B0(1),B0(2),'o');
plot(B0(1),B0(2),'.');
plot([A0(1) A1(1) B1(1) B0(1)],[A0(2) A1(2) B1(2) B0(2)]);
fill([A1(1) P(1) B1(1)],[A1(2) P(2) B1(2)],'c')
plot([A1(1) P(1)],[A1(2) P(2)]);
plot([A1(1) P(1) B1(1)],[A1(2) P(2) B1(2)],'r.');
plot(P(1),P(2),'r*')
plot(x,y)
axis equal
x2=a*cos(fi)+e*cos(-mu-delta+alfa);
y2=a*sin(fi)+e*sin(-mu-delta+alfa);
plot(x2,y2,':')
B12=A1+d*[cos(-mu(1)-delta(1));sin(-mu(1)-delta(1))];
P2=[x2(1);y2(1)];
plot([A0(1) A1(1) B12(1) B0(1)],[A0(2) A1(2) B12(2) B0(2)],':');
plot([A1(1) P2(1) B12(1)],[A1(2) P2(2) B12(2)],'g:')
plot([A1(1) P2(1)],[A1(2) P2(2)],'g:');
plot([A1(1) P2(1) B12(1)],[A1(2) P2(2) B12(2)],'g.');
plot(P2(1),P2(2),'g*')