#include "stdafx.h"
float v,p,m,i;
float a[2000],t,c,b,d,y1,y2,y3,y4; //全局变量
void main()
{
printf("请分别输入V0,P,m的值");
scanf("%d%d%d",&v,&p,&m);
printf("请输入dt的值");
scanf("%f",&t);
printf("请分别输入C,ρ,A的值");
scanf("%f%f%f",&c,&d,&b);
float f(float e ); //子函数声明
a[0]=v;
for(i=0;i<200/t;i++)
{
y1=f(a[i]); //调用子函数过程时a[i]不变
y2=a[i]+0.5*f(y1)*t;
y3=a[i]+0.5*f(y2)*t;
y4=a[i]+f(y3)*t;
a[i+1]=a[i]+(y1+2*y2+2*y3+y4)*t/6;
}
for(i=0;i<200/t;i++)
{
printf("%f\n",a[i]);
}
}
float f(float e) //子函数用来计算导数 ,返回导数值(传递参量)
{
f=0;
f=p/(m*e)-c*d*b*e*e/m;
return f;
}