xref: /plan9/sys/src/cmd/eqn/integral.c (revision 3e12c5d1bb89fc02707907988834ef147769ddaf)
1*3e12c5d1SDavid du Colombier #include "e.h"
2*3e12c5d1SDavid du Colombier #include "y.tab.h"
3*3e12c5d1SDavid du Colombier 
4*3e12c5d1SDavid du Colombier extern int Intps;
5*3e12c5d1SDavid du Colombier extern double Intht, Intbase, Int1h, Int1v, Int2h, Int2v;
6*3e12c5d1SDavid du Colombier 
integral(int p,int p1,int p2)7*3e12c5d1SDavid du Colombier void integral(int p, int p1, int p2)
8*3e12c5d1SDavid du Colombier {
9*3e12c5d1SDavid du Colombier 	if (p1 != 0)
10*3e12c5d1SDavid du Colombier 		printf(".ds %d \\h'%gm'\\v'%gm'\\*(%d\\v'%gm'\n", p1, -Int1h, Int1v, p1, -Int1v);
11*3e12c5d1SDavid du Colombier 	if (p2 != 0)
12*3e12c5d1SDavid du Colombier 		printf(".ds %d \\v'%gm'\\h'%gm'\\*(%d\\v'%gm'\n", p2, -Int2v, Int2h, p2, Int2v);
13*3e12c5d1SDavid du Colombier 	if (p1 != 0 && p2 != 0)
14*3e12c5d1SDavid du Colombier 		shift2(p, p1, p2);
15*3e12c5d1SDavid du Colombier 	else if (p1 != 0)
16*3e12c5d1SDavid du Colombier 		bshiftb(p, SUB, p1);
17*3e12c5d1SDavid du Colombier 	else if (p2 != 0)
18*3e12c5d1SDavid du Colombier 		bshiftb(p, SUP, p2);
19*3e12c5d1SDavid du Colombier 	dprintf(".\tintegral: S%d; h=%g b=%g\n", p, eht[p], ebase[p]);
20*3e12c5d1SDavid du Colombier 	lfont[p] = ROM;
21*3e12c5d1SDavid du Colombier }
22*3e12c5d1SDavid du Colombier 
setintegral(void)23*3e12c5d1SDavid du Colombier void setintegral(void)
24*3e12c5d1SDavid du Colombier {
25*3e12c5d1SDavid du Colombier 	yyval = salloc();
26*3e12c5d1SDavid du Colombier 	printf(".ds %d %s\n", yyval, lookup(deftbl, "int_def")->cval);
27*3e12c5d1SDavid du Colombier 	eht[yyval] = EM(Intht, ps+Intps);
28*3e12c5d1SDavid du Colombier 	ebase[yyval] = EM(Intbase, ps);
29*3e12c5d1SDavid du Colombier 	lfont[yyval] = rfont[yyval] = ROM;
30*3e12c5d1SDavid du Colombier }
31