www.mathworks.com feval(symengine,'simplify',YOUR_EXPRESSION_HERE,'IgnoreAnalyticConstraints')

a = 10; b = 10; N = 20; x = linspace(-a,a,N); y = linspace(-b,b,N); [X,Y] = meshgrid(x,y); f = @(n) (-1)^((n-1)/2)/n^3*sinh(n*pi*Y/(2*a))/cosh(n*pi*b/(2*a))*sin(n*pi*X/(2*a)); omega = 1; Z = omega*(X*Y-32*a^2/pi^3*(f(1)+f(3)+f(5))); %n=5で…

外圧を受ける円管 ラーメの式 clear; close all; a=1; b=2; p=1; r = linspace(a,b,20); sr = -b^2*p/(b^2-a^2)*(1-a^2./r.^2); ss = -b^2*p/(b^2-a^2)*(1+a^2./r.^2); plot(r,sr); title('\sigma_r'); figure(); plot(r,ss); title('\sigma_\theta'); 回転…