#include <stdlib.h>

double *vector(long size)
/* 配列v[]用 動的メモリ確保 */
{
	double *v;

	v=(double *)malloc((size_t) ((size+1)*sizeof(double))); // 1個余分に確保
	return v; // サイズがsで先頭アドレスがaなら,a〜a+s-1の範囲が使って良い領域です。
}

void free_vector(double *v)
/* メモリ領域開放 */
{
	free((char*) v);
}