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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
1 j. H: H3 ]6 _0 Y#include "stdafx.h"
/ j( k, R2 G- U
& Y0 Z5 c7 A9 L#include <stdio.h>
#include <stdlib.h>
#include <string.h>
7 C  s7 M7 z- P  r! s  F
7 Y  s; w9 b8 f1 g  z  y#include "mathlink.h"
7 ^7 o* g- n3 w3 i6 a1 W$ v
double f(MLINK lp,double x);

8 R4 l9 M) @0 h8 m' Pint main()
, U5 F* J# y1 C4 Y% E{
        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
6 s! K7 {8 H) L' G1 r        MLEnvironment env;//mathlink的环境变量
9 V  m! n2 U- {" }9 K" ]        MLINK lp;
+ n( s  M$ b! a; ^0 m" V! A+ P4 ?        env=MLInitialize(NULL);//初始化
5 F( G: t# q) `8 v        if(env==NULL)
% K1 b3 V9 M# {$ l& g* O! [        {
                printf("initialization error!\n");
2 V" U) e, B  [5 l! ^                return 1;
        }
        lp=MLOpen(argc, argv);
        if(lp==NULL)
        {
                printf("open error!\n");
                return 1;
2 b; _. _) ]& i1 s9 n9 B: ?) o        }
1 ~' \) _$ A( F  C% x, H        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
6 r, B: Z4 j$ a! r, N- w8 ?% ^0 z        MLDeinitialize(env);

# P7 |$ h+ C7 ?1 l' g+ p' Q        //printf("Hello World!\n");
        return 0;
3 G. O0 I1 @7 A- d! l3 ?5 c}
# k/ i. t7 e8 K* p0 a

' q6 J" e) ?& Ddouble f(MLINK lp, double x)
* f2 p% _+ j) m& U4 @* p4 Y{
5 w+ d% \0 L. p5 h0 a' l        double res;
        int pkt;
, K6 X0 q( @4 n2 ~% f. R        static int done=0;
9 \* y; t0 @& _        if(!done)
        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
' q( d! u9 Z/ x0 j# S             MLEndPacket(lp);
' w0 i, x$ T1 g( X4 _+ k3 K                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
; p5 P2 S! B# {- w  O$ @                MLNewPacket(lp);
                done=1;
        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
" ~  @" d. v0 g5 A/ K0 s        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
3 j& f7 h& q- R( P) y        MLNewPacket(lp);
0 Z& `* N9 s/ x. T        MLGetReal(lp, &res);
0 z+ }4 k! w/ C        return res;
}
# k6 H. I5 b& p

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？