mathematica子函数共享

nyy

1. Clear["Global`*"];(*清除所有变量*)
   
2. deg=Pi/180;(*角度制下1°所对应的弧度*)
   
3. (*子函数,利用三边计算角的余弦值,角是c边所对的角*)
   
4. cs[a_,b_,c_]:=((a^2+b^2-c^2)/(2*a*b))
   
5. (*子函数,四面体体积公式，a,b,c分别是从一个顶点出发的三条棱，x,y,z分别是对棱*)
   
6. fun[a_,b_,c_,x_,y_,z_]:=Sqrt[Det[{{0,1,1,1,1},{1,0,a^2,b^2,c^2},{1,a^2,0,z^2,y^2},{1,b^2,z^2,0,x^2},{1,c^2,y^2,x^2,0}}]/288]
   
7. (*子函数,海伦公式,利用海伦公式计算三角形的面积*)
   
8. heron[a_,b_,c_]:=Module[{p=(a+b+c)/2},Sqrt[p*(p-a)*(p-b)*(p-c)]]
   
9. (*子函数,已知△ABC的a b c三边长度,求c这条边上的中线长度*)
   
10. zx[a_,b_,c_]:=Sqrt[(a^2+b^2-c^2/2)/2]
    
11. (*子函数,计算角平分线长度,求的是角C的平分线(也是c边上的角平分线)*)
    
12. ca[a_,b_,c_]:=Module[{p=(a+b+c)/2},2/(a+b)*Sqrt[a*b*p*(p-c)]]
    
13. (*miller rabin子函数*)
    
14. MR[n0_,a0_]:=Module[{n=n0,a=a0,s,m,t1,k},
    
15.     s=0;m=n-1;While[Mod[m,2]==0,m=m/2;s=s+1];
    
16.     t1=PowerMod[a,m,n];
    
17.     If[t1==1,Return[True]];
    
18.     k=0;While[k<s-1&&t1!=n-1,k=k+1;t1=Mod[t1^2,n]];
    
19.     If[t1==n-1,Return[True],Return[False]]
    
20. ]
    
21. Do[If[Not@MR[nn,k],Print[{k}]],{k,1,307}]
    
22. 
    

复制代码

有可能会经常用到的，就写成子函数，共享一下，也是备份一下

王守恩

1. (*子函数,四面体体积公式，a,b,c分别是从一个顶点出发的三条棱，x,y,z分别是对棱*)
   
2. fun[a_,b_,c_,x_,y_,z_]:=Sqrt[Det[{{0,1,1,1,1},{1,0,a^2,b^2,c^2},{1,a^2,0,z^2,y^2},{1,b^2,z^2,0,x^2},{1,c^2,y^2,x^2,0}}]/288]

复制代码

1,{a,b,c,x,y,z}={8,19,14,4,4,9},   四面体体积=?
2,四面体的6条边都是正整数，体积=33。来一个？

nyy

1. (*子函数,给定两个点,计算通过这两个点的直线的斜率.输入两点:{a,b},{c,d},返回:(d-b)/(c-a)*)
   
2. kk[pt1_,pt2_]:=Module[{dpt},dpt=pt2-pt1;dpt[[2]]/dpt[[1]]]
   

复制代码

nyy

nyy 发表于 2023-8-28 09:40

1. (*子函数,圆内接四边形面积公式(婆罗摩笈多公式)*)
   
2. plmjd[a_,b_,c_,d_]:=Module[{p=(a+b+c+d)/2},Sqrt[(p-a)*(p-b)*(p-c)*(p-d)]]
   

复制代码