编程计算子午圈曲率半径，卯酉圈曲率半径和纬线圈半径，我编的程序没有错误可是计算不出来，求指教~

#include<stdio.h>

#include<math.h>

#define a 6378245

#define b 6356863

void main()

{

int M,N,r,x

double y,z,v,w

double f

printf(" 子午圈曲率半径 卯酉圈曲率半径 纬线圈半径\n")

f=sqrt((a*a-b*b)/a*a)

y=f*f

for(x=0x<=60x++)

{

z=sin(x)*sin(x)

v=sqrt(1-y*z)

w=(1-y*z)

M=a*(1-y)/(w*v)

N=a/v

r=a*cos(x)/v

printf("纬度=%d,%d, %d, %d\n",x,M,N,r)

}

getchar()

}

计算模型

f(x)=sin(x)+x^2曲率计算

计算方法

利用符号运算

算法流程

（1）用符号表示表示f(x)

（2）求f(x)的一阶导数，二阶导数

（3）根据曲率半径定义r=|(1+y'^2)^(3/2)/y''|

计算代码。

clcclear all

syms x

f=sin(x)+x^2

f1=diff(f,x)

f2=diff(f1,x)

r=(1+f1^2)^(3/2)/f2

xx=0:0.1:30rr=zeros(size(xx))

MATLAB 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。

MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂(矩阵实验室)。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式，代表了当今国际科学计算软件的先进水平。

MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。

---