子程序
中相应的哑元S确设置成了real*8,这可能就是不一致的原因。(不知道你把uz_int设置成什么?)
再有
kxy_integral(k_end,k_start,uz_int)
应该是
call
kxy_integral(k_end,k_start,uz_int)
#include <stdio.h>#include <math.h>
void main()
{
int i=1
int n,m
for (i=1i<100000i++)
{
m=sqrt(100+i)
n=sqrt(268+i)
if(m*m==i+100)
if(n*n==i+268)
printf("%ld ",i)
}
}