Easy MatLab Q - ODE 15S

Question

我正在尝试解决微分代数方程并继续得到以下错误：

脚本中的函数定义必须出现在文件的末尾。 将“pscalc”函数定义之后的所有语句移动到第一个本地函数之前 定义

%prompt = 'Enter the radius of the sphere (m): ';
%R =  input(prompt);

%prompt = 'Enter the area of the outlet of the sphere (m): ';
%A =  input(prompt);

g = 9.81;  % m2/s
R = 10;    %radius
A = 0.1; % outlet area m^2

y0 = [0; 19.99];
tspan = [0 3600];

function out=pscalc(t,y)
out = [y(1)*((pi/A)*y(2)-(R/A)-g*y(2))
      (A/(pi*((y(2))^2 - (2*R*y(2))))*((y(1))^2+ 2*g*y(2)))];
end

[t,y]=ode15s(@pscalc, tspan, y0);

不确定如何修复它。

任何帮助都是适用的。