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

遥望科西嘉

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

% U/ P; T8 d1 e5 H#include <stdio.h>
#include <stdlib.h>
2 C1 u' o# U3 y' m. X7 t#include <string.h>
1 @5 L1 w- f1 J0 N. I
#include "mathlink.h"
* Q' U2 ~; O. y9 i; J. [  S& U3 O
( r* S- n1 i1 D/ r. c% r2 r! edouble f(MLINK lp,double x);

& P  f( P. H9 q# Tint main()
! e) U3 t/ W% b: N{
2 f$ G/ A, Z7 H+ K, e+ n        int argc=4;
, s4 U/ W( w/ }8 e! q4 {, `: X/ x        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
0 y3 T. |: G( ?" p; t        MLEnvironment env;//mathlink的环境变量
5 a1 L6 n4 j4 ]1 y& q1 V7 ]        MLINK lp;
2 O  Z6 u/ o* a2 d& K# V* c! X        env=MLInitialize(NULL);//初始化
        if(env==NULL)
        {
* F2 S; z7 ^" X9 u! ?! V+ R                printf("initialization error!\n");
$ W5 v: S5 O, m                return 1;
        }
6 t( E+ @+ K2 v        lp=MLOpen(argc, argv);
4 q# N% y$ g8 K# L9 [9 Y4 A' r5 O        if(lp==NULL)
# k4 A2 Y8 X) Q' j2 g! ^8 f0 N$ H        {
  V5 G* p8 q  q2 H                printf("open error!\n");
                return 1;
( g  C1 W2 M, Q* s7 b" x        }
& _' `  B5 S, w$ F, S        for(int i=1;i<=10;i++)
5 s1 p) `5 H- {. n9 N2 H                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
        MLDeinitialize(env);

( }! J" Y0 t  O) D/ z5 B# I        //printf("Hello World!\n");
) u3 [2 x: i: I; ~* p7 a        return 0;
8 S, F2 r' l/ s}

( \& a; y( n9 F5 _& ^0 J# q
double f(MLINK lp, double x)
3 a0 E. q7 P9 d" P# d) k{
; I- |* h; q4 C8 T/ Q  c        double res;
. X. A( h# a; @        int pkt;
5 \# ?5 g, j* h$ z+ [1 r* E9 B        static int done=0;
        if(!done)
- q' B: h6 H" v$ H  X        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
* @1 I6 a! I4 s3 a. u+ x6 G                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
                done=1;
        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
  ^' i. }$ [' w& ^        MLPutFunction(lp,"F",1);
2 E$ U) g! h: W9 ?% B: q        MLPutReal(lp,x);
        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
+ V! T; f/ t; o+ @- x7 W; i( W        MLGetReal(lp, &res);
; h: g) B# ?$ z0 C! J        return res;
8 m& ~8 `; j! f% G}
  J, {$ H; r/ V) g) o8 f, O

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？