Question

怎么用matlab解这个微分方程

初始条件a(0)＝0

Answer 1

求解常微分方程（组），一般可以用dsolve函数求其解析值，但有的微分方程只能通过ode函数求出其数值解。题主给出的微分方程，如用dsolve函数求解，则会出现错误，所以该方程可以用ode45函数得到其数值解。求解过程如下：

a0=0;

tspan=[0,5];

fun=@(t,a)2/(-3+2*(1-a)^(-1/3)+4*(1-a)^(-2/3));

[t,a]=ode45(fun,tspan,a0);

A=[t,a];

disp(' t a(t)')

disp(A)

plot(t,a);

xlabel('t'),ylabel('a(t)')

Answer 2

把等号右边的那个式子含有a的复杂的分母乘过去呀！然后再把等号左边那个分母dt乘过来。再两边积分。