clear;clc;
A=0.0001;B=0.0002;l=0.04;Cs=1.6e-5;L=3.6e-7;C=2e-7;U0=2000;R0=0.06;
g0=(2.5e9)*(U0/(2e+11*L))^0.1;
s=(2.5e8)*(U0/(2e+11*L))^0.1;
[x,y,z]=dsolve('Dx=z^2/Cs^2','Dy=z','Dz=(1/L)*(U0-z*(R0+(l/Cs)*(A*(1-sech(x/g0))+B*exp(-((x-g0)/s)^2)))-y/C)','x(0)=0','y(0)=0','z(0)=0','t') |
回复此楼
