#include <stdlib.h>
#include "real.h"
#define rnd (real)rand()/RAND_MAX

void pseudo(int n,real a[],real b[], real x[])
{
int i; for(i=0;i<n;i++) x[i]=a[i]+rnd*(b[i]-a[i]);
}