#include <math.h>
#include <stdio.h>

float r(float x)
{
	return exp(x);
}

float f(float x)
{
	return sin(x)*(1.0+exp(x));
}

void main ()
{
	void femlag(float [], float [], int, float (*)(float),
					float (*)(float), int, float []);
	int n,i,order;
	float pi,x[21],y[21],e[7],rho,d;

	printf("FEMLAG delivers:\n");
	for (n=10; n<=20; n += 10) {
		e[1]=e[4]=1.0;
		e[2]=e[3]=e[5]=e[6]=0.0;
		pi=3.14159265358979;
		for (i=0; i<=n; i++) x[i]=pi*i/n;
		printf("N=%2d\n",n);
		for (order=2; order<=6; order += 2) {
			femlag(x,y,n,r,f,order,e);
			rho=0.0;
			for (i=0; i<=n; i++) {
				d=fabs(y[i]-sin(x[i]));
				if (rho < d) rho=d;
			}
			printf("     ORDER=%1d    MAX.ERROR= %7.3e\n",order,rho);
		}
	}
}