解差分方程C语言代码

#include<stdio.h>
void main(){
//x(k+1)=M*x(k);M=[x1,x2] ;x(k)=[y1]
//                             [x3,x4]       =[y2]
double x1 = 0.95, x2 = 0.03, x3 = 0.05, x4 = 0.97;
int step = 1;
double y1 = 350000, y2 = 650000,y3=0.0;
while (1){
printf("%d ", step);
y3 = x1*y1 + x2*y2;
y2 = x3*y1 + x4*y2;
y1 = y3;
step++;
printf("\ny1=%f\ny2=%f", y1, y2);
getchar();
}


}