level 2
相距1647
楼主
想请教一下大家,这个程序跑完了为什么每次都显示0.2不是一个有效变量
a=0;
b=1;
R21[y_]=c11*y^7+c12*y^6+c13*y^5+c14*y^4+c15*y^3+c16*y^2+c17*y+c18 ; (*y<x*)
R22[y_]=c21*y^7+c22*y^6+c23*y^5+c24*y^4+c25*y^3+c26*y^2+c27*y+c28 ; (*y>x*)
d11[y_]=D[R21[y],y];
d12[y_]=D[R21[y],{y,2}];
d13[y_]=D[R21[y],{y,3}];
d14[y_]=D[R21[y],{y,4}];
d15[y_]=D[R21[y],{y,5}];
d16[y_]=D[R21[y],{y,6}];
d17[y_]=D[R21[y],{y,7}];
d21[y_]=D[R22[y],y];
d22[y_]=D[R22[y],{y,2}];
d23[y_]=D[R22[y],{y,3}];
d24[y_]=D[R22[y],{y,4}];
d25[y_]=D[R22[y],{y,5}];
d26[y_]=D[R22[y],{y,6}];
d27[y_]=D[R22[y],{y,7}];
solution=Solve[{R21[x]==R22[x],d11[x]==d21[x],d12[x]==d22[x],d13[x]==d23[x],d14[x]==d24[x],d15[x]==d25[x],d16[x]==d26[x],d27[x]-d17[x]==1,R21[a]-10*d11[a]==0,10*d17[0]-d16[a]+10*R21[0]+d11[0]==0,d12[a]+d15[a]==0,d13[a]-d14[a]==0,R22[b]+d21[b]==0,d24[b]==0,d25[b]==0,d26[b]+d27[b]==0},{c11,c12,c13,c14,c15,c16,c17,c18,c21,c22,c23,c24,c25,c26,c27,c28}];
solution=FullSimplify[solution];
u[x_]=3x^3-4 x^2;
u1[x_]=D[u[x],x];
u2[x_]=D[u[x],{x,2}];
u3[x_]=D[u[x],{x,3}];
yzh[y_]=u[0]*R21[0]+u1[0]*d11[0]+d12[0]*u2[0]+d13[0]*u3[0]++;
c11=solution[[1,1,2]];
c12=solution[[1,2,2]];
c13=solution[[1,3,2]];
c14=solution[[1,4,2]];
c15=solution[[1,5,2]];
c16=solution[[1,6,2]];
c17=solution[[1,7,2]];
c18=solution[[1,8,2]];
c21=solution[[1,9,2]];
c22=solution[[1,10,2]];
c23=solution[[1,11,2]];
c24=solution[[1,12,2]];
c25=solution[[1,13,2]];
c26=solution[[1,14,2]];
c27=solution[[1,15,2]];
c28=solution[[1,16,2]];
Simplify[yzh[x]]
R21[x_,y_]=Simplify[R21[y]]
R22[x_,y_]=Simplify[R22[y]]
R[x_,y_]=Piecewise[{{R22[y],0<=x<y<=1},{R21[y],0<=y<x<=1}}](*If[y>x,R22[x,y],R21[x,y]]*)
r1[x_,y_]=1/Gamma[0.9] Integrate[(y-s)-0.1 (D[R[x,s],{s,2}])s,{s,0,y},Assumptions->0<y<1]+(2y+6)(D[R[x,y],x])
r2[x_,y_]=1/Gamma[0.9] Integrate[(x-t)-0.1 (D[r1[x,t],{t,2}])t,{t,0,y},Assumptions->0<x<1]+(2*x+6)(D[r1[x,y],x])
r2[0.2,0.2]
2021年05月16日 08点05分
1
a=0;
b=1;
R21[y_]=c11*y^7+c12*y^6+c13*y^5+c14*y^4+c15*y^3+c16*y^2+c17*y+c18 ; (*y<x*)
R22[y_]=c21*y^7+c22*y^6+c23*y^5+c24*y^4+c25*y^3+c26*y^2+c27*y+c28 ; (*y>x*)
d11[y_]=D[R21[y],y];
d12[y_]=D[R21[y],{y,2}];
d13[y_]=D[R21[y],{y,3}];
d14[y_]=D[R21[y],{y,4}];
d15[y_]=D[R21[y],{y,5}];
d16[y_]=D[R21[y],{y,6}];
d17[y_]=D[R21[y],{y,7}];
d21[y_]=D[R22[y],y];
d22[y_]=D[R22[y],{y,2}];
d23[y_]=D[R22[y],{y,3}];
d24[y_]=D[R22[y],{y,4}];
d25[y_]=D[R22[y],{y,5}];
d26[y_]=D[R22[y],{y,6}];
d27[y_]=D[R22[y],{y,7}];
solution=Solve[{R21[x]==R22[x],d11[x]==d21[x],d12[x]==d22[x],d13[x]==d23[x],d14[x]==d24[x],d15[x]==d25[x],d16[x]==d26[x],d27[x]-d17[x]==1,R21[a]-10*d11[a]==0,10*d17[0]-d16[a]+10*R21[0]+d11[0]==0,d12[a]+d15[a]==0,d13[a]-d14[a]==0,R22[b]+d21[b]==0,d24[b]==0,d25[b]==0,d26[b]+d27[b]==0},{c11,c12,c13,c14,c15,c16,c17,c18,c21,c22,c23,c24,c25,c26,c27,c28}];
solution=FullSimplify[solution];
u[x_]=3x^3-4 x^2;
u1[x_]=D[u[x],x];
u2[x_]=D[u[x],{x,2}];
u3[x_]=D[u[x],{x,3}];
yzh[y_]=u[0]*R21[0]+u1[0]*d11[0]+d12[0]*u2[0]+d13[0]*u3[0]++;
c11=solution[[1,1,2]];
c12=solution[[1,2,2]];
c13=solution[[1,3,2]];
c14=solution[[1,4,2]];
c15=solution[[1,5,2]];
c16=solution[[1,6,2]];
c17=solution[[1,7,2]];
c18=solution[[1,8,2]];
c21=solution[[1,9,2]];
c22=solution[[1,10,2]];
c23=solution[[1,11,2]];
c24=solution[[1,12,2]];
c25=solution[[1,13,2]];
c26=solution[[1,14,2]];
c27=solution[[1,15,2]];
c28=solution[[1,16,2]];
Simplify[yzh[x]]
R21[x_,y_]=Simplify[R21[y]]
R22[x_,y_]=Simplify[R22[y]]
R[x_,y_]=Piecewise[{{R22[y],0<=x<y<=1},{R21[y],0<=y<x<=1}}](*If[y>x,R22[x,y],R21[x,y]]*)
r1[x_,y_]=1/Gamma[0.9] Integrate[(y-s)-0.1 (D[R[x,s],{s,2}])s,{s,0,y},Assumptions->0<y<1]+(2y+6)(D[R[x,y],x])
r2[x_,y_]=1/Gamma[0.9] Integrate[(x-t)-0.1 (D[r1[x,t],{t,2}])t,{t,0,y},Assumptions->0<x<1]+(2*x+6)(D[r1[x,y],x])
r2[0.2,0.2]