Mathematica或者Matlab解Gamma方程
来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/07/15 03:45:39
Mathematica或者Matlab解Gamma方程
![](http://img.wesiedu.com/upload/2/c4/2c421d42d8bb156996690737107f95dc.jpg)
这个值小于等于0.05 解出n=1到n=500的所有t的值,该怎么操作.会求一个n下的t也行!
k=0.5
![](http://img.wesiedu.com/upload/2/c4/2c421d42d8bb156996690737107f95dc.jpg)
这个值小于等于0.05 解出n=1到n=500的所有t的值,该怎么操作.会求一个n下的t也行!
k=0.5
![Mathematica或者Matlab解Gamma方程](/uploads/image/z/290765-29-5.jpg?t=Mathematica%E6%88%96%E8%80%85Matlab%E8%A7%A3Gamma%E6%96%B9%E7%A8%8B)
要求:matlab R2008b版本以前,之后的版本无法调用maple内核
因为你要的是范围,这个比较讨厌,所以我把你这个函数的图形画出来,你自己看着办,要求正好等于0.05的解可以对这个函数用fsolve或fzero求解,比如t=fzero(inline(ft-0.05),10)%10是迭代的初值,因为有时有两个解,见图,所以可能需要手动换初值才能得到想要的解
clear;clc;
n=11;
ft=sym(maple(['gamma(',num2str(n),'-1,0.5*t)/gamma(',num2str(n),'-1)']));
ezplot(ft,[-20,100]);ylim([-pi,pi]);set(gca,'ytick',[.05]);grid on;
![](http://img.wesiedu.com/upload/e/83/e8334c13fd0821911f853ec308cdee5f.jpg)
因为你要的是范围,这个比较讨厌,所以我把你这个函数的图形画出来,你自己看着办,要求正好等于0.05的解可以对这个函数用fsolve或fzero求解,比如t=fzero(inline(ft-0.05),10)%10是迭代的初值,因为有时有两个解,见图,所以可能需要手动换初值才能得到想要的解
clear;clc;
n=11;
ft=sym(maple(['gamma(',num2str(n),'-1,0.5*t)/gamma(',num2str(n),'-1)']));
ezplot(ft,[-20,100]);ylim([-pi,pi]);set(gca,'ytick',[.05]);grid on;
![](http://img.wesiedu.com/upload/e/83/e8334c13fd0821911f853ec308cdee5f.jpg)