C / C++ / MFC

// plain old array return
double* MyFunc()
{
    double* pd = new double(SIZE OF ARRAY);
    // DO SOMETHING WITH ARRAY
    // THEN RETURN ARRAY
    return db;
}

void MyOtherFunc()
{
    double* pd = MyFunc();
    // DO SOMETHING WITH ARRAY
    // THE DELETE ARRAY
    delete dp;
}

void bessj0(double x[20],double j0[20])			
//Returns the Bessel function j0(x) for any real x. 
{
	double ax[20],z[20],xx[20],y[20],ans1[20],ans2[20];

	for (int i=0;i<20;i++)
	{
		if ((ax[i]=fabs(x[i]))<8.0)		//Direct rational function fit
		{
			y[i]=x[i]*x[i];
			ans1[i]=57568490574.0+y[i]*(-13362590354.0+y[i]*(651619640.7
				+y[i]*(-11214424.18+y[i]*(77392.33017+y[i]*(-184.9052456)))));
			ans2[i]=57568490411.0+y[i]*(1029532985.0+y[i]*(9494680.718
				+y[i]*(59272.64853+y[i]*(267.8532712+y[i]*1.0))));
			j0[i]=ans1[i]/ans2[i];
		}
		else
		{
			z[i]=8.0/ax[i];
			y[i]=z[i]*z[i];
			xx[i]=ax[i]-0.785398164;
			ans1[i]=1.0+y[i]*(-0.1098628627e-2+y[i]*(0.2734510407e-4
				+y[i]*(-0.2073370639e-5+y[i]*0.2093387211e-6)));
			ans2[i]=-0.1562499995e-1+y[i]*(0.1430488765e-3
				+y[i]*(-0.6911147651e-5+y[i]*(0.7621095161e-6
				-y[i]*0.934945152e-7)));
			j0[i]=sqrt(0.636619772/ax[i])*(cos(xx[i])*ans1[i]-z[i]*sin(xx[i])*ans2[i]);
		}
	}
	//return j0[i];
	//return 0;
}

bessj0(p,j0);

void Integrals(double Ip[20])	
{	
	double alpha[20],beta[20];
	for (int t=0;t<20;t++)
	{
		double j0[20],j1[20],H0[20],H1[20],p[20];

		void bessj0(double x[20],double j0[20]);	
		void bessj1(double x[20],double j1[20]);
		//int t = 1;

		alpha[t] = beta[t] = 2*t*PI;	// t = 1
		p[t] = 2*sqrt(pow(alpha[t],2.0)+pow(beta[t],2.0));
		
		bessj0(p,j0);
		bessj1(p,j1);


		H0[t] = sum0(p[t]);//Watch this!(***)
		H1[t] = sum1(p[t]);

		
		Ip[t] = ((PI*p[t])/2) * ((2*p[t]*j0[t]) + (((PI*p[t])/2) * ((H0[t]*j1[t])-(H1[t]*j0[t]))) - j1[t]);//WATCH THIS: WOULD IT BE PROPER TO USE 'j0[t]' in this line? the function call is :bessj0(p,j0), but I need to use the values of j0[t] from the function. Is this appropriate?

	//	return 0;//Ip[20];
	
	}
}

double sum0(double x)	//using the power series expansion, till the 11th power
{
	double H0;
	H0 = x - (pow(x,3)/9) + (pow(x,5)/(9*25)) - (pow(x,7)/(9*25*49))
		+ (pow(x,9)/(9*25*49*81)) - (pow(x,11)/(9*25*49*81*121));

	return H0;
}

// I forget but this may be passing by
// value, which is not what you want.
void Integrals(double Ip[20])
{
    ....
    for( i=0; i<20; ++i ) // SIZE KNOWN
        Ip[i] = what ever;
}

// This is closer to normal passing
void Integrals(double *Ip)
{
    ....
    for( i=0; i<20; ++i ) // SIZE KNOWN
        Ip[i] = what ever;
}

// This is safer and will not require a rewrite
// if array size (ever) changes.
void Integrals(double *Ip, int nSize)
{
    ....
    for( i=0; i<nSize; ++i ) // SIZE UNKNOWN
        Ip[i] = what ever;
}

#define WINVER       0x0501
#define _WIN32_WINNT 0x0501
#define _WIN32_IE    0x0600

int *ji=new int;<br />
	*ji=10;<br />
	::SendMessage((HWND)this->m_lHandle,WM_USER+90,0,(LPARAM)*ji);<br />
	CString s;<br />
	s.Format("%d",*ji);<br />
	MessageBox(s);

        int * i=(int*)lParam;<br />
         char Text[100];<br />
	  sprintf(Text,"%d ",*i);<br />
	 //char **Nam2=(char**)wParam;<br />
	  ::MessageBox(NULL,Text,Text,MB_OK);<br />
<br />
        *i=40;

int *ji = new int;
*ji = 10;
::SendMessage((HWND)this->m_lHandle,WM_USER+90,(LPARAM)ji);
...
int *i = (int *)lParam