내것도 아닌데 오지랖때문에 ㅋㅋㅋ

일단 식은 이러해

초기조건이 h(0)=5

그래서 내가 한번 해봤는데

clc
clear all
Cd=0.5;
r=3;
pi=3.14;
A=r^2*pi;
g=9.81;
h=dsolve('Dh= (-Cd*A*sqrt(2*g))/(pi*(2*r*h^(0.5)-h^(1.5)))','h(0) = 5')

이러면 내 예상으로는 h에 대해서 t변수만 남기고 죄다 변수처리 될 줄 알았는데

h =
 RootOf(15*t*Cd*A*g^(1/2)+10*_Z^(3/2)*2^(1/2)*pi*r-3*_Z^(5/2)*2^(1/2)*pi-50*5^(1/2)*2^(1/2)*pi*r+75*5^(1/2)*2^(1/2)*pi)

이렇게 떠버리네

Z는 어디서 튀어나온걸까 뭐하는놈일까,, 접근 방법이 잘못된걸까?

누가 시간에따른 h의 식을 푸는 것좀 해줄 개럴없을까??

-------------------

추가

clc
clear all

Q = quad(@myfun,5,0);

%% myfunc.m%% 요건 그냥 m파일 따로 만들어놓은거
function t = myfun(h)
Cd=0.5;
r=3;
pi=3.14;
A=r^2*pi;
g=9.81;
t=(pi*(2*r*h^(0.5)-h^(1.5)))/(-Cd*A*(2*g)^0.5);
%%

요렇게하면 h를 5부터 0까지 적분해서 t가 나올줄 알았는데 애러가 뜨네,,, 으어,,,알려줘,, 울프 뭐시기 안되,, 무조건 맬랩이야...답변고마워 개럴