数学建模之倾倒的啤酒杯
代码
clear;clc
a=10;k=0.3
w1=200;w2=230;w3=50;
r=15;s3=1;e=0.3;Pi=3.1415926;
for i=4:7
k=i/10;
f=@(x)(w1*3*(x^4-k^4)*x^3/(4*(1-k^3))+2*r*e*Pi*a*(1-
k^3)/3+w3*s3)/(w1*(x^3-k^3)/(1-k^3)+w2+w3)
fplot(f,[0,1]);
hold on;
end
ylabel("质心高度");
xlabel("酒的液面高度");
legend("k=0.4","k=0.5","k=0.6","k=0.7");