你好,你把你的ode45那一行改为:
flag=, n=, B_b=, I_p=, r_r= ;%这些数据在这里定义好 [t,y] = ode45(@(t,y)equation(t,y,flag,n,B_b,I_p,r_r), tspan, y0)
有帮助望采纳呢,非常感谢