我试图通过该程序找出拉格朗日的插值。我用数组解决了它,但是当使用动态内存分配时,程序给了我垃圾结果。
#include<stdio.h>
#include<conio.h>
#define SIZE 100
int main()
{
float *x,*y;
float value = 0,ask,temp;
int i,j,n;
printf("Enter size");
scanf("%d",&n);
x = (float*)malloc(n*sizeof(float));
y = (float*)malloc(n*sizeof(float));
for(i = 0; i < n;i++)
{
printf("x[%d]: ",i);
scanf("%f",(x+i));
printf("y[%d]: ",i);
scanf("%f",(y+i));
}
printf("Enter value to find");
scanf("%f",&ask); //cin >> ask;
for(i = 0; i < n;i++)
{
temp = 1;
for(j = 0; j < n; j++)
{
if(i != j)
{
temp = temp * (ask-(*(x+i))/(*(x+i)-*(x+j)));
}
}
value = value + temp * *(y+i);
}
printf("%f",value);
}
答案 0 :(得分:7)
您需要#include <stdlib.h>
,因为这是声明执行动态内存分配的函数的标头(malloc()
等)。
假设您使用C编译器而不是C ++编译器,您将发现malloc()
调用上的类型转换不是必需的。
x = (float*)malloc(n*sizeof(float)); /* (float *) is unnecessary */
问题是,如果没有stdlib.h
,编译器会假定malloc()
返回int
。类型转换可能允许代码在没有stdlib.h
的情况下进行编译,但结果是后续使用malloc()
d指针将具有未定义的行为,因为指针不一定能够在往返过程中存活(转换为int
然后回来。)
如果您使用的是C ++编译器,则需要(float *)
类型转换AND #include <stdlib.h>
以避免未定义的行为。