谁有使用FORTRAN90程序做高斯积分的例子？

有一组 fortran 的通用程序，可以在库中找到。

公共公共变量的数据类型声明语句和公共声明语句必须在所有使用它的主程序或子程序（模块）中重新写入，使用它的人必须编写它...

另外，common 用法是一种已经被淘汰的低效且不安全的写法，如果你已经使用了 module 的现代结构，你应该使用 module 方法，那么为什么要混入向后 common？ 事实上，你可以直接在主程序中编写 use global（下一行隐式 none），你根本不需要 common

这是 fortran95 编程中 p205 的一个示例问题，我也遇到过这个问题，研究过，我认为是：模块中定义的值不能以子程序参数的形式调用。 调用后，将发生此错误。

如果你在大型程序中遇到这样的问题，你可以使用调用 sub（） 的形式，如示例中所示。 如果遇到合理的解释，希望大家能一起讨论。

% Gauss-Legendre 数值积分。

fun=@(x)sin(x).*x.^3;%fun：积分表达式。

A、B：点的上限和下限。

a = 0;

b = pi;

tol=1e-8;%tol：积分精度，默认 1e-6% 计算正交节点。

syms x

p=sym2poly(diff((x^2-1)^(n+1),n+1))/2^n*factorial(n));

tk=roots(p);正交节点。

计算正交系数。

ak=zeros(n+1,1);

for i=1:n+1

xkt=tk;

xkt(i)=[

pn=poly(xkt);

fp=@(x)polyval(pn,x)/polyval(pn,tk(i));

ak(i)=quadl(fp,-1,1,tol);正交因子。

end%积分变量模具熔透替换，[a，b]变量猜测作弊为[-1,1]xk=（b-a） 2*tk+（b+a） 2;

测试了积分集群兆慢函数FUN的有效性。

fun=fcnchk(fun,'vectorize');

计算变量替换后积分函数的值。

fx=fun(xk)*(b-a)/2;

计算点值。

ql=sum(ak.*fx)

quadl（fun，0，pi） % 调用 MATLAB 的内部积分函数。

公共公共变量的数据类型声明语句和公共声明语句必须在所有使用它的主程序或子程序（模块）中重新写入，使用它的人必须编写它...

另外，常见的用法是一种已经淘汰的低效、不安全的编写方式，如果你已经用过模块的现代结构，就应该使用模块方法，那么何必费心混入后向常见的看其实可以直接在主程序中写 use global（下一行隐式无）， 而且你根本不需要普通

c 正确。 因为 2>3 的结果是合乎逻辑的。

没错t.和f.它们都不是真的。

逻辑常量应为true.和flase.（两边都有小数点）。

该计划的要求是什么。

理想。

对于没有定积分的函数，找一个现成的数值计算库，比如IMSL。

只是写在手机上，没有编译。

integer num,sum

num=0sum=0

11 do i=103,193,10do j=2,i-1

if(mod(i,j)==0) goto 11enddo

num=num+1

sum=sum+i

write(*,i

enddowrite(*,num

write(*,sumend