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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
* u- m9 e+ z% V$ b3 J#include "stdafx.h"
! D8 Q$ a3 ?7 x2 W$ f+ r# w
3 R7 {. j3 b5 p#include <stdio.h>
6 \5 b: D' w: N+ v7 }8 ~3 y  A; Z#include <stdlib.h>
# [7 c1 B+ b1 `7 L" U' ^#include <string.h>

#include "mathlink.h"

4 E/ _. G: B( z! Fdouble f(MLINK lp,double x);

5 u- Y5 G5 r9 E$ mint main()
9 J$ D: }# s# m{
" N  |0 j8 c5 w( V; H        int argc=4;
6 W' f/ e- d+ Y, z9 ^6 d9 |        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
1 {: T' |, n) M0 ~! ^/ E6 p        MLEnvironment env;//mathlink的环境变量
. |9 |8 G$ w( o% E6 n: b4 H        MLINK lp;
0 U$ W3 ~% B4 p. N2 i        env=MLInitialize(NULL);//初始化
        if(env==NULL)
( N% E0 r9 l4 O5 Z        {
                printf("initialization error!\n");
                return 1;
        }
        lp=MLOpen(argc, argv);
        if(lp==NULL)
5 D* h+ B6 c! X+ f' q6 H        {
                printf("open error!\n");
                return 1;
! ]* E7 k& l: s        }
        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
2 q0 s/ T- j3 z0 q        MLDeinitialize(env);
/ h3 C$ ?9 K" ?
        //printf("Hello World!\n");
        return 0;
}
% {0 a% s& r  L* c7 H9 A

double f(MLINK lp, double x)
  u; v; h7 s+ i- x* y2 i( h# J{
        double res;
        int pkt;
+ j. ?; k7 S& r" }2 K        static int done=0;
        if(!done)
        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
4 r, M, c& L" t2 q                MLPutFunction(lp,"ToExpression",1);
5 R/ X4 C6 |$ P6 V                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
1 f2 c( d( ?5 H0 z                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
$ ?4 Y7 Z' O0 j# O. l                MLNewPacket(lp);
+ W( `( l9 \- V) \  g                done=1;
        }
        MLPutFunction(lp,"EvaluatePacket",1);
' ~$ ^( ~4 y" o9 P        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
* s, A8 k  b6 u2 F9 x/ d0 t        MLNewPacket(lp);
/ E8 v; z( g3 O  I        MLGetReal(lp, &res);
        return res;
) r) S" }$ u, [1 }: h( ]& y9 l}
/ ^; l" G5 U6 L) B2 P

遥望科西嘉

这么大个论坛没人了么？

遥望科西嘉

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