Fortran编写龙贝格公式

1.programRiceProblemintegerNgrid,NumPerTonrealTonNumdoubleprecisionSumNgrid=64NumPerTon=5d7Sum=0d0do

EXTERNALFSX0=1.5EPS=1.0E-06CALLDNEWT(X0,EPS,FS,L)IF(L.NE.0)THENWRITE(*,10)X0ENDIF10FORMAT(1X,'X=',E

运行一下看提示信息,如果单从给出的程序来看,类似15.0(x-30)的写法不和语法,运算符号是不能省略的

错误原因如ls至于“但当我把所有的i和x交换,其他不变,就没有问题,能够的到正确的结果.”是因为fortran的I-N规则,以字母IJKMLN开头的变量,如果没有定义,默认为整形,交换以后,虚参I是整

公式：=IF(A1

加个“=”,然后写需要的公式,或者插入函数

程序如下：（复制就能运行）t=-5:0.001:5;w=(1-t.^2/10).*exp(-t.^2/2).*cos(0.5.*t);%所要编写的公式plot(t,w)%该公式对应的函数图gridon

C++的行么?Romberg.h：文件//Romberg.h:interfacefortheRombergclass./////////////////////////////////////////

%输入一个数据后回车,接着输入第二个数据fprintf('请输入24个X的值(输入一个数据后回车)\n');for i=1:24    

A=0.027*B**(-1/7)再问：为什么B值我无论输入多大，A值都等于0.027，即B的负7分之1次方总是等于1，这是怎么回事了再答：这个好办，改为A=0.027*B**(-1.0/7.0)，变

IMPLICTNONEREAL::A(4,4),B(4,4),C(4,4)=0INTEGER::I,J,KREAD(*,*)A,B!输入方式也可改为读文件DOI=1,4DOJ=1,4DOK=1,4C(

直接for循环呗再问：？？？FORTRAN也能用for循环再答：可以的

比如a%b表示结构体变量a的b成员类似于C语言的a.b,结构体（struct）inpKey%nKeysinpKey是一个结构体,你可以搜索一下TypeinpKey找到它的定义nKeys是它的一个结构体

去年学的Fortran语言,许多算法都忘了.给你一个去年写的牛顿插值程序,需要改动的地方自己改,貌似循环计算（多层计算）那块处理得不错.real*8x(0:5),y(0:5),c(0:5),u1,u2

代码如下,供参考：      program test      im

programbaiduinteger::I,J,Nreal*8::Cr,Treal*8,dimension(:),allocatable::P,XN=3!变量X的个数Cr=5.0d0!常量Cr,自己

不知道你的三维网格需要建成什么形状?如果告诉tecplot需要读的文件的具体格式的话是挺好编写的.再问：三维网格形状是六面体，tecplot读取文件格式貌似是.F90谢谢再答：tecplot可以调用f

表示很难理解你的需求.再问：就是每十分钟有N个数，然后提取第6,12,18...时刻的数（即一个小时，两个小时，三个小时）。在每个时刻数第一行加上对应的时间和相应的步数，然后每个时刻的数加上1到N的编

你的问题可进行简化,fortran中可以直接用matmul函数进行矩阵的乘积,下面是我编的程序愿您满意!其中矩阵C为矩阵A的转置,有问题可以问我programmainimplicitnoneinteg

=IF(R7=0.9,R7