A. 超等数学 黎曼重积分定理证明 3
7. 交换积分次序即可,∫[0->x]∫[0->y]g(t)dtdy的积分次序是
先把t从0到y积分,即0≤t≤y,然后在把y从0到x积分,0≤y≤x
交换次序后便是先对y积分,而由上可知 t≤y≤x ∴先把y从t->x,
然后在对t积分,而t的积分下线是0,上限是x,∴再把t从0->x
即有 ∫[0->x]∫[0->y]g(t)dtdy=∫[0->x]∫[t->x]g(t)dydt
=∫[0->x](x-t)g(t)dt
B. 用matlab中,在解微分方程组时,dydt=zeros(2,1) 这个式子表示什么意思
解微分方程组的时候后安排位置用的的,zeros(2,1)设置了一个两行一列的矩阵,后面的第一个方程Dy(1)就是这个矩阵的第一行,以此类推。如果方程组有5个方程就是zeros(5,1)了。
C. matlab中dydt函数的用法
y(1) = y
y(2) = y'
dydt 输出的是一个列向量
dydt(1)=y'
dydt(2)=y''
D. matlab 定解问题求解
#include<iostream.h>
#include<math.h>
#include<stdlib.h>
void main()
{ double u[16][16],x[16];
double h=0.0625,r=0.5,y;
int a=1,i,j;
y=r*h*h/a;
for(i=0;i<=16;i++) {x[i]=i*h;
u[i][0]=0; }
for(j=0;j<=8;j++)
{ u[0][j]=0;
u[16][j]=0; }
u[8][0]=u[8][0]+0.5*0.5*0.5*0.5*0.5*0.5*0.5*0.5*0.5*0.5;
for(j=0;j<=8;j++)
{
for(i=1;i<=15;i++)
{ u[i][j+1]=r*u[i+1][j]+(1-2*r)*u[i][j]+r*u[i-1][j];
}
}
for(j=1;j<=8;j++)
{for(i=1;i<=16;i++)
{cout<<u[i][j]<<" ";
}
cout<<"\n";
}
}