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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
#include "stdafx.h"
/ w3 {; |  f/ g# r1 q& i  l
- ]. ?5 t9 F# I' d: \. A  u! ?- E#include <stdio.h>
( X  h: T; m/ ^8 l& @* }#include <stdlib.h>
#include <string.h>
# r& f/ T. ^2 L/ ?
5 E/ p- D) p; i* g; l#include "mathlink.h"
' p) G$ ]1 R% W- }# V& S; ?( ]
; m; O; B& x% l4 E5 b1 Mdouble f(MLINK lp,double x);

4 N, n( N9 g, E, W: c5 o$ Eint main()
: e3 B4 E( P2 c1 z5 ~1 o{
        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
9 ^6 }+ ~; r' v: O        env=MLInitialize(NULL);//初始化
3 X; {( o2 k$ F8 R( k% \% r- m        if(env==NULL)
) b% M! z, k( W4 d        {
                printf("initialization error!\n");
* @7 p, Q8 n: |$ B& p  \2 g                return 1;
        }
        lp=MLOpen(argc, argv);
, D' ?! U' {% o- E; }        if(lp==NULL)
        {
                printf("open error!\n");
  j* W6 D4 x, W                return 1;
/ q" Q- x% G$ {" \2 o        }
        for(int i=1;i<=10;i++)
: s* O5 q0 h+ I+ D  W& X" ?& ^                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
        MLDeinitialize(env);
, {( s# o2 v8 x8 ?7 ^
        //printf("Hello World!\n");
        return 0;
5 ^( L2 N' W5 N- P6 F/ M+ o4 W" g}


double f(MLINK lp, double x)
, V6 Z( f3 ^, f" K{
1 E5 J7 q9 W4 g/ P* {( c6 T        double res;
        int pkt;
        static int done=0;
        if(!done)
        {
                //char* expression="cos(x)^5";
% V9 m+ |5 i. G9 C8 @                MLPutFunction(lp,"EvaluatrPacket",1);
; x; c$ o* T0 t+ D* W/ V7 I                MLPutFunction(lp,"ToExpression",1);
5 Q8 U: j# m+ r" a: W" U                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
) F# m* h: O8 J$ x9 S                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
  Y( M7 K+ w$ a6 @; V; a# a2 t                MLNewPacket(lp);
                done=1;
, i3 J  ?: B; X* z9 \3 }        }
# u3 M5 Q- k- B: x        MLPutFunction(lp,"EvaluatePacket",1);
7 t! H7 `: o2 b5 p: f% q        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
        MLEndPacket(lp);
/ q$ N* @' B  i% _8 e        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
' Y  }5 _: t! ^( ~; M1 Y        MLGetReal(lp, &res);
) t( J/ g% N2 N0 X. }& N4 D        return res;
}

  z7 k, ~; y* s& }5 N: f8 L

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？