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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
#include "stdafx.h"
- K4 {) ]0 `: Q8 a' m0 L) i* b
- K# G% X# t- g4 }' J#include <stdio.h>
#include <stdlib.h>
( x0 y2 Y' c; N4 P#include <string.h>

#include "mathlink.h"
! W* h3 G6 x: s: \/ c% Z6 `
double f(MLINK lp,double x);

; P% g5 i3 U  e. X6 h7 dint main()
5 B& ^! T) Y+ E  Z2 H7 y4 ]& ~4 \{
" E: b5 @% c4 P0 B0 b        int argc=4;
+ c5 ~0 R9 x! b& Q        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
& ]6 Z" f0 S1 v4 O) L        MLEnvironment env;//mathlink的环境变量
' Q* {, j8 J# E* w5 |        MLINK lp;
        env=MLInitialize(NULL);//初始化
        if(env==NULL)
        {
                printf("initialization error!\n");
: i- T" b8 i2 h' e; W6 t                return 1;
8 E  s0 Q6 [# z, t* N  g: V# ?) K        }
# ]& w  g# X% x3 f: p- T7 z        lp=MLOpen(argc, argv);
        if(lp==NULL)
4 R' O' i; m0 v0 |) d# c6 o        {
5 ?; V8 E. R5 D  [                printf("open error!\n");
( G1 |; C6 D: U  @2 b; B. P                return 1;
5 z( \, w% k$ S3 f* N        }
        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
+ z/ s  D( E, g5 `
. H; t2 o3 B$ p& G7 e3 l8 A% c' o* Y        MLClose(lp);
        MLDeinitialize(env);

        //printf("Hello World!\n");
' E# X1 s- {" i- V        return 0;
* |7 y# A1 h& X1 K}
! H$ n/ m- i& y" z
& E3 a0 L; s! o/ V" {. j7 C
double f(MLINK lp, double x)
1 p( g8 i# u3 }$ ]) h9 m1 ~7 V{
        double res;
        int pkt;
        static int done=0;
        if(!done)
        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
% B; B5 x" \- e1 ~9 e                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
8 W5 R* S9 n  F                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
( S- e" B! A5 W                MLNewPacket(lp);
                done=1;
% ]7 q+ \4 n/ s  y# Q: M* j" ]  d        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
+ I% g9 N! L, w5 O        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
3 [6 t& n+ |& K( J2 W- S( w3 @        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
        MLGetReal(lp, &res);
  }, h3 O7 i( Q: y4 U: v6 O+ d        return res;
}

# a- C3 F/ U1 D& w) h

遥望科西嘉

这么大个论坛没人了么？

遥望科西嘉

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