#include <stdio.h>
#define is ==

void basis(
	int c,
	double t,
	int np,
	int x[],
	float n[]
);

void knot(
	int n,
	int c,
	int x[]
);

void rbsurf(
    float b[],
	int k,
	int l,
	int npts,
	int mpts,
	int p1,
	int p2,
	int *p_itest,
	int ibnum,
	float bold[],
	float niku[],
	float mjlw[],
	float rsumij[],
 	float q[]
);

float sumrbas(
	float b[],
	float nbasis[],
	float mbasis[],
	int npts,
	int mpts
);