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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
#include "stdafx.h"

4 P  ~8 r5 _5 b' s3 J7 u#include <stdio.h>
#include <stdlib.h>
6 y6 Q) e6 u4 ?* X5 P* q#include <string.h>

#include "mathlink.h"

double f(MLINK lp,double x);
4 a/ O' J% w$ K. u3 i( o) H4 Y" ~, E
" l+ C2 A# y4 @; R' b; p' o/ U% u4 n+ Uint main()
, D* H+ r- ~( l, Y{
        int argc=4;
* K; x2 R. g. u5 n; D9 R7 R        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
2 q. j  Q) r8 R# S- s        MLEnvironment env;//mathlink的环境变量
  H0 m. Z2 j$ X# S! n- p% ?7 W        MLINK lp;
        env=MLInitialize(NULL);//初始化
1 G5 R9 L. C: g* N        if(env==NULL)
$ C8 m7 h9 H, `  \7 g. }/ d) N        {
                printf("initialization error!\n");
: {9 J+ ]  |6 I; s# s) e- W                return 1;
" N/ z4 [! s# Y& v4 C$ m- v9 ?8 q        }
        lp=MLOpen(argc, argv);
3 _9 U; F- f8 Q) t( R        if(lp==NULL)
, u# D3 U" f% O$ F; S        {
- n7 X) l* S4 S6 l: M) |3 N                printf("open error!\n");
                return 1;
1 o7 A  r* P3 O" n        }
* Z# \; @- O' Q8 K. f5 y. Q        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
  S: `# W' z) r3 G, W" U- I! J- O        MLDeinitialize(env);
( H: i1 b6 q: d8 P3 \( L
        //printf("Hello World!\n");
        return 0;
}


double f(MLINK lp, double x)
$ E& P" s- K( r4 z( u{
        double res;
0 c  L) j- l* J4 X        int pkt;
        static int done=0;
, ]5 p- b/ L% Z% z( r7 M$ q        if(!done)
" w% q8 ~! ]: i3 b: X4 ]) p2 p        {
0 n2 B! `+ o$ R                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
) c$ d" \4 ^  ?0 s8 l             MLEndPacket(lp);
* q% _' y* |4 |6 D, h                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
7 g5 g) [0 {4 F, D$ A                done=1;
  V0 u( i+ ~! F' ~- ]$ x' j, a        }
        MLPutFunction(lp,"EvaluatePacket",1);
' |- \& c$ n) `9 F1 y5 l        MLPutFunction(lp,"N",1);
' x; Y  M. F. Y+ S3 r4 v" T/ Q        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
' ]; v. D, Z+ b8 B+ y8 Q        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
0 Z8 j* V4 g+ ^' f% \        MLNewPacket(lp);
' f1 j( j5 u$ g! U        MLGetReal(lp, &res);
( J. W' I, E$ q7 m3 P% Y% x  s        return res;
}

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？