用ode求解二阶非线性微分方程

输入参数数目不足求大神帮助

 

clear all;
close all;
clc;
f0=2*10^4; 
R0=6*10^(-6); 
P0=1.013*10^5;
time=1/f0;
Pa=1.2*10^5; 
Ki=1.67; 
rho=1000; 
sigma=0.0725; 
mu=0.001; 
Pi=3.1415926;
c=1481; 
a=R0/10; 
w=2*Pi*f0; 
P1=1*10^5;
P2=1*10^5;
P3=1*10^5;
P4=1*10^5;
step=time/30000; 
ttime=1*time;
ts=0:step:ttime;
y0=[R0,0];
[t,y]=ode45(fun5,ts,y0); 
plot(t/time,y(:,1)/R0)
function f=fun5(t,y)
f0=2*10^4; 
R0=6*10^(-6); 
P0=1.013*10^5; 
time=1/f0;
Pa=1.2*10^5; 
Ki=1.67; 
rho=1000; 
sigma=0.0725; 
mu=0.001; 
Pi=3.1415926;
c=1481; 
a=R0/10; 
w=2*Pi*f0; 
f=[y(2);((P0+2*sigma/R0)*(R0/y(1))^(3*Ki)-2*sigma/y(1)-4*mu*y(2)/y(1)-(P0-Pa*sin(w*t)))/rho/y(1)-(3/2*(y(2)^2))/y(1)-3*(P0+2*sigma/R0)*Ki*R0^(3*Ki)*y(2)*y(1)^2/(rho*c*(y(1)^3-a^3)^(Ki+1))+(Pa*w*cos(w*t))/(rho*c)];
end