mathematica c语言混合编程的问题 大神求救

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
! l, i! M6 i: Z% Z' S+ I#include "stdafx.h"

#include <stdio.h>
#include <stdlib.h>
" Y9 T; O  O) e+ U#include <string.h>
0 a0 R1 O/ @2 x  c
#include "mathlink.h"
5 x# }, n1 v+ b- U& N: k
double f(MLINK lp,double x);

: F1 }# r7 G( G) p8 D0 Fint main()
{
+ k: j) [+ N, |1 @        int argc=4;
' K9 }) e. x& _        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
        MLEnvironment env;//mathlink的环境变量
. v5 w. Z) K2 h8 x        MLINK lp;
        env=MLInitialize(NULL);//初始化
& z/ \" @  D2 q0 N( z# p        if(env==NULL)
        {
                printf("initialization error!\n");
                return 1;
        }
        lp=MLOpen(argc, argv);
% q( B; }5 v$ ~1 j4 B        if(lp==NULL)
& M( n" G; f6 C: M        {
                printf("open error!\n");
                return 1;
        }
/ u  C! _9 a9 @( o+ o+ H5 m        for(int i=1;i<=10;i++)
/ u( i( Q1 \1 |8 f                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
- X$ v* Z0 V  b5 h5 \& `; _
4 K& u6 A' O8 c5 j8 @$ {        MLClose(lp);
        MLDeinitialize(env);
1 o$ X3 Q% `7 I$ H" X! Y
1 h+ M7 a5 e; d( q7 n        //printf("Hello World!\n");
6 K& Y: D! H2 X2 K- v        return 0;
}


double f(MLINK lp, double x)
{
        double res;
+ t* X) a2 K! l1 ?- L        int pkt;
8 @5 h: z, z7 l2 h" J        static int done=0;
3 C1 l$ U/ c' t  j        if(!done)
        {
4 x% s9 O6 B6 Q+ k3 \' G3 Y. y                //char* expression="cos(x)^5";
4 e6 O5 o* O  j) R, S" K                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
; @, f9 V4 h6 e4 W. j             MLEndPacket(lp);
0 i9 k% l+ i: ~% k! m                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
                done=1;
        }
; S! g& [; P' P        MLPutFunction(lp,"EvaluatePacket",1);
9 O) y& M3 a, B8 o9 _+ J2 I6 B' l% f        MLPutFunction(lp,"N",1);
. a/ u, W* D8 |, C% N6 V6 L        MLPutFunction(lp,"F",1);
9 \; F1 ]/ i; j6 G$ K) J        MLPutReal(lp,x);
/ E. u( P1 o5 ]2 H4 y" z/ r  T        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
) C! ~7 _& i$ u; }* h0 c' |        MLGetReal(lp, &res);
# I6 z% m1 O- o) y        return res;
}

% _, b- C2 @6 e7 y- r1 c

遥望科西嘉

这么大个论坛没人了么？

遥望科西嘉

我觉得是不是因为在mathematica上面还要写个接口，提供c数据交换的接口